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SUMMARY 

This report contains the detail information for the analytical and the experimental pro- 
grams described in the report "Study of Effects of Design Details on Structural Response 
to Acoustic Excitation, " NASA CR- 1959 (reference 1) and is intended as a supplement 
to that report. Since this report is a supplement the contents of reference 1 are contin- 
uously quoted, and the reader must refer to the original text for complete continuity. 

INTRODUCTION 


The analytical and experimental program described in reference 1 was concerned with 
the application of the finite element displacement method for the prediction of dis- 
placements and strain distributions for the normal mode vibration of flat stiffened panels. 
Panels with stiffeners in one direction (one dimensional panels) and panels with ortho- 
gonal stiffening (two-dimensional panels) were considered. Details of the stiffener and 
its attachment to the plate were taken into account. In the case of the one-dimensional 
panel, elastic edge conditions on the boundary along the panel width were considered 
analytically. The material presented in reference 1 is essentially a summary of the 
analytical approach and the experimental results, and this report provides essential 
documentation for persons interested in more detailed information. 

The stiffness and consistent mass matrix for both the stiffener model and the plate ele- 
ment developed in reference 1 is discussed and presented in detail in the appendix. 

An outline of the computer programs is presented along with pertinent details of the 
analysis. Detailed experimental data from all the specimens described in reference 1 
is discussed and presented in tabulated form. 



ANALYTICAL PROGRAM 


One aspect of the analytical program described in reference 1 considered the 
derivation of a beam element to represent a thin-walled open-section stiffener as is| 
usually encountered in aircraft structure. This element was used to model the stiffeners 
for the panel configurations under consideration. For the one-dimensional panel con- 
figuration, the stiffener warping coordinate was taken to be zero since warping is an 
odd function along the stiffener length and is zero at the panel centerline for the assumed 
fundamental mode in the direction of the stiffener (see equations 20a and 20b, reference 
1) . For the two-dimensional panel analysis the stiffener warping coordinate was taken 
equal to the panel twist to insure compatibility of slope between the stiffener element 
and the plate element. The two-dimensional panel configuration required consideration 
of a coordinate transformation to describe stiffeners parallel to the x-axis and the y-axis, 
and since it was necessary to consider rotation of the stiffener about a general point on 
the stiffener profile line (the 'attach' point), it was also necessary to develop trans- 
formations for the elastic forces from the shear center to the attach point and for the 
inertia forces from the centroid to the attach point (equations 10 through 13, reference 
1) . All of these transformations have been carried out and are summarized as a com- 
posite stiffness and mass matrix for the stiffener element. The adjective 'composite' 
is used to denote the use of logic numbers to compute the stiffness or the mass matrix 
for a stiffener parallel to the x-axis or the y-axis. The composite stiffness and mass 
matrices are presented in Appendix A. The notation and the sign conventions are 
as described in reference 1 . 

The rectangular plate bending element described in reference 1 was based upon the 16 
degree-of-freedom plate bending element described by Bogner, Fox, and Schmidt 
(reference 2). The modification introduced in reference 1 was to consider an internal 
mode for the element in the form of clamped -clamped beam functions described by the 
coordinate W q (equation 23, ref. 1). The introduction of the coordinate W q resulted 
in the definition of modifying terms for the basic stiffness and mass matrices. These 
modifying terms (equations 26 and 27, ref. 1) are presented in Appendix B. The 
notation and sign convention are as described in reference 1 . 

COMPUTER PROGRAMS 

The basic computer program flow chart for the one-dimensional panel analysis is pre- 
sented in figure 1. The program computes, sequentially, the element properties for a 



bay of structure and assembles the element in the free-free stiffness and mass 
matrices by application of displacement compatibility and equilibrium conditions 
at each element node. The desired elastic supports are assembled in the free-free 
stiffness and mass matrices introducing the elastic constraints. If a lumped mass is 
desired, the data is introduced as a lumped support with zero stiffness. The stiffness 
and mass matrices are non-dimensionalized, and kinematic constraints are applied 
at either end of the structure as desired. The kinematic constraint is of the form of 
a clamped support at either or both ends of the structure and is realized computa- 
tionally by deleting the row-column terms in the stiffness and mass matrices corres- 
ponding to the constrained coordinates and appropriately reordering the stiffness and 
mass matrices. The eigenvalue problem is formulated and the eigenvalues and eigen- 
vectors are obtained using standard routines (reference 3) based upon the Jacobi's 
method (reference 4). The non-dimensionalizing parameters and the eigenvalues 
and corresponding eigenvectors are printed. If mode shapes, shear, and bending 
moment distributions are desired, the values are computed, normalized to the maximum 
value, and printed. Values for the displacement, shear, and bending moment for 
points interior to an element are computed using equations 16, 17, and 18 of 
reference 1 . 

The basic flow chart for the two-dimensional panel computer program is presented in 
figure 2. The program computes the system stiffness characteristics, non-dimensionalizes 
the matrix, removes the constrained coordinates (c lamped-edges) , and reassembles the 
stiffness matrix. The plate stiffness is computed first and then the rib stiffness is com- 
puted. Each element is introduced into the free-free system by applying displacement 
compatibility and equilibrium conditions at each grid point. The consistent mass 
matrix is assembled identically to the stiffness matrix, the eigenvalue problem is 
formulated, and the eigenvalues and eigenvectors are obtained as previously described. 
Figure 3 illustrates the stiffened plate in plan view showing the bay (plate element) 
and the rib nomenclature. The plate stiffness and mass matrices are assembled in the 
sequence indicated by the plate bay number, and the rib stiffness and mass matrices 
are assembled in the sequence indicated by the rib number and rib segment number as 
indicated in figure 3. The assembly of the rib and the plate elements at the inter- 
section of two orthogonal ribs is illustrated in figure 4. The positive coordinate 
directions are as indicated. Computer program listings, flow charts, and descriptions 
for the one-dimensional and two-dimensional panel arrays are presented in Appendix C. 
The necessary information for data input and program output is presented in Appendix D. 
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EXPERIMENTAL DATA 


One of the objectives of the program described in reference 1 was to provide data 
for comparison with the analytical results. The technique used to determine mode 
shapes and strain distributions is described in reference 1 . Comparison of theory 
and experiment is given in reference 1 for frequencies, mode shapes, and strain 
(bending moment) distribution (in the case of one-dimensional panels). 

Each specimen was mounted in the test frame and cork particles were sprinkled on the 
specimen. The specimen was excited by discrete frequency sinusoidal excitation using 
a specially designed speaker enclosure as described in reference 1 . Frequency sweeps 
were conducted for four speaker phase conditions (ref. 1) . The predominant modes as 
indicated by the Chladni patterns formed by the cork particles were photographed. 
These patterns, for the indicated specimen, frequency, and speaker phase condition, 
are presented in Appendix E. 

For the one-dimensional panel specimens, mode shapes were determined using two (2) 
accelerometers. One accelerometer was fixed in position for a reference value and 
the other accelerometer was stepped in position along the centerline of the panel . 

At each position of the stepped accelerometer, amplitude and phase of both accelero- 
meters (as observed on an oscilloscope) was recorded. The accelerometer positions 
for the one-dimensional specimens is given in figure 5. The accelerometers were 
calibrated to give identical output for a given input, but no attempt was made to 
force an absolute output since only relative acceleration (displacement) was desired. 
Data reduction was accomplished by determining the accelerometer output in millivolts 
for both accelerometers, dividing the value at a position by the reference value for 
the positions, and then normalizing the set of data to the largest value of the ratio 
(not necessarily the reference position) . The normalized acceleration data (mode 
shape) for the specimen, the indicated position (figure 5), frequency, and speaker 
phase conditions are tabulated in Tables 1 through 8. A minus sign as a value indicates 
a 180° phase shift with respect to the reference value and an asterisk denotes a 90° 
phase shift. To compare the normalized experimental values with the computed 
normalized mode shapes, the ratio of the experimental value and the computed value 
at a point on the structure was determined, all experimental values for the mode 
were multiplied by this ratio, and the data plotted. The plotted data comparison 
for the one-dimensional panel specimens is given in figures 22 through 38 of reference 
1 . 



Strain measurements in the direction of the panel length (perpendicular to the ribs) 
was accomplished by placing fifteen strain gages along the panel as indicated in 
figure 5. Detailed location of the strain gages is indicated in figure 18 of ref- 
erence 1 . Data reduction for the strain measurements was accomplished as described 
for the acceleration data. The strain gage system was calibrated so that one millivolt 
of output corresponded to 417 microinches per inch of strain. Strain gage output in 
millivolts and phase with the indicated reference is tabulated in Tables 9 through 16 
for the indicated specimen, strain gage, frequency, and speaker phase condition. 

The plotted comparison between the experimental and calculated values for strain 
are given in figures 22 through 38 of reference 1 . 

As described in reference 1, experimental determination of mode shapes for the two- 
dimensional panel specimens was more difficult than for the one-dimensional specimens. 
A detailed experimental mode investigation was possible only for the machined panel 
specimen. For specimens SP 1 1- T and SP 11-2, only the basic phase relationship 
between adjacent panel bays could be determined. Acceleration measurements for 
the machined panel specimen were taken at the locations illustrated in figure 6. 
Tabulated values for accelerometer output in millivolts with phase relative to the 
reference are given in Table 17 for the indicated position, frequency, and, 
speaker phase condition. For the machined panel specimen the strain gage locations 
are indicated in figure 7 with the exact location indicated by the (x, y) coordinate 
position given in Table 18. 

For specimens SP 1 1 — 1 and SP 11-2 the tabulated values for acceleration in millivolts 
for the center of each panel bay is given in Table 19. Strain gage location 
and nomenclature for specimens SP II -1 and SP 11-2 are given in figures 8 and 9, 
respectively, with the exact location tabulated. Strain measurements for the indi- 
cated specimen, frequency, and speaker phase condition are given in Tables 20 
through 22. The strain gage calibration was suchthat417 microinches per inch equaled 
one millivolt of strain gage output. 

Damping was measured for selected strain gages by determining the logarithmic 
decrement from the photograph of the decaying strain signal (ref. 5) . The specimen 
was excited in a given mode with the selected strain gage signal displayed on an 
oscilloscope. The excitation was suddenly stopped and the decaying strain signal 
photographed with a camera mounted on the oscilloscope. The logarithmic decre- 
ment and the damping ratio were determined from the photograph. For the indicated 
specimen, strain gage, and frequency, values of the damping ratio (percent of 
critical damping) are given in Table 23. 


5 




FIGURE 1. FLOWCHART: PROGRAM BMPROP(MAIN)/ONE-DIMENSIONAL 
PANEL ARRAYS (CONTINUED) 




A 


PRINTED' 
YES / STIFFNESS 

MASS MATRIX. 


NO 


l 


WRITE NUMBER OF FINAL 
COORDINATES, NON- 
DIMENSIONALIZING 
PARAMETERS, STIFFNESS 


WF 

NONDIMENS 
PARA A 

LITE 

ilONALIZING 

PETERS 

AND MASS MATRICES 






STEP 8 


3 


COMPUTE AND 
WRITE EIGENVALUES 
AND EIGENVECTORS 
IN DESCENDING 
ORDER 


STEP 9 


ARE 

YES M0DE SHAPES, 

SHEAR, AND BENDING 
•M0MENT DISTRIBUTION- 
DESIRED? 


NO 


COMPUTE AND WRITE NORMALIZED 
DISPLACEMENTS, SLOPES, SHEAR 
AND BENDING MOMENT DISTRIBUTIONS 
ALONG CENTERLINE OF THE PANEL 


STEP 10 



FIGURE 1. FLOW CHART: PROGRAM BMPROP(MAIN)/ONE-DIMENSIONAL 
PANEL ARRAYS (CONCLUDED) 
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FIGURE 2. FLOWCHART: PROGRAM PLTVIB(MAIN)AWO-DIMENSIONAL 
PANEL ARRAYS (CONTINUED) 


8 




i SEGMENT 2, PLTSTF | 
^SEGMENT 3, ASSYPJ 


CALCULATE 
PLATE STIFFNESS 
FOR PANEL BAY 


PLTSTF 


STEP 5 


/ ASSEMBLE \ 

'PARTITIONED PLATE' 
STIFFNESS MATRIX 
IN SUPER MATRIX, R, 

\ ” ASSYP / 


STEP 6 


J SEGMENT 5 i 
| RIB, ASSYR J 


COMPUTE PARTITIONED 
RIB STIFFNESS MATRIX 

FOR^m element; 

RIB 


ASSEMBLE RIB 
STIFFNESS MATRIX 
IN SU£ER_MA]_RIX, R 

ASSYR 


HAVE ALL RIBS 

BEEN CONSIDERED? 


COMPUTE STIFFNESS 
PARAMETER, TK 


HAVE ALL BAYS ' 
BEEN CONSIDERED?,. 


STEP 7 


STEP 8 


STEP 9 


STEP 10 


T " SEGMENT 6 ~”~1 

NONDIM, ORDER, FILL, DELETE] 

1 

/nondimensionalize free-free\ 
f _ _ STIFFNESS [MATRIX, R \ STEP 1 1 

\ NONDIM 7 

/ APPLY EDGE CONSTRAINTSV CTCn 
/ AND REORDER MATRIX, R ) STEP 12 

\ ORDER, DELETE "/ 

~ V 

/ FILL CONSTRAINED STIFFNESS \ 

/ MATRIX, SF, FROM REORDERED \ STEP 13 
\ SUPER MATRIX, R / 

\ FILL / 


SET IOP = 2 
GO TO STEP 4 


FIGURE 2. FLOWCHART: PROGRAM PLTVIB(MAIN)/TWO-DIMENSIONAL 
PANEL ARRAYS (CONTINUED) 


9 




10 




STEP 17 


STEP 18 


STEP 19 


STEP 20 


STEP 21 


STEP 22 


FIGURE 2. FLOWCHART: PROGRAM PLTVIB(MAIN)AWO-DIMENSIONAL 
PANEL ARRAYS (CONTINUED) 


11 





FIGURE 2. FLOW CHART: PROGRAM PLTVIB(MAIN)AW0-DIMENSI0NAL 
PANEL ARRAYS (CONCLUDED) 
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FIGURE 6. ACCELEROMETER LOCATION FOR 

NINE-BAY MACHINED PANEL SPECIMENS 
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TABLE 1 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-1 

SPI-1 

SPI-1 

SPI-1 

SPI-1 

SPI-1 

Frequency 

82 Hz 

106 Hz 

119 Hz 

126 Hz 

108 Hz 

127 Hz 

Speaker Cond . 

A 

A 

A 

B 

A 

B 

Position 1 

0 

0 

0.20* 

0 

- 

- 

2 

0.75 

0.70 

0.25* 

-0.50 

- 

- 

3 

0.25 

0 

0.15* 

0 

- 

- 

4 

-0.05 

0 

0.40 

-0.24 

0.25 

-0.18 

5 

0 

0.95 

1.00 

-1 .00 

1 .00 

-1 .00 

6 

-0.10 

0 

0.43 

0.10* 

0.55 

-0.40 

7 

0.25 

0 

0.63 

0.19* 

0.20 

-0.30 

8 

0.75 

0.61 

1.00 

0.27* 

0.21* 

-0.17 

9 

0.25 

0 

0.60 

0.29* 

0.25 

-0.15 

10 

-0.10 

0 

0.42 

0.35* 

0.65 

-0.20 

11 

0 

1.00 

0.90 

0.81 

1.00 

0.43 

12 

-0.10 

0 

0.35 

0.35 

0.35 

0 

13 

0 

0 

0.45 

0.31 

- 

- 

14 

1 .00 

0.94 

0.50 

0.50 

- 


Position 15 

0 

0 

0.05 

0 

- 

— 


*90° Phase shift to reference 
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TABLE 2 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-2-1 

SPI-2-1 

SPI-2-1 

SPI-2-1 

SPI-2-1 

SPI-2-1 

SPI-2-1 

Frequency 

88 Hz 

92 Hz 

103 Hz 

113 Hz 

120 Hz 

91 Hz 

109 Hz 

Speaker Cond . 

A 

A 

B 

A 

A 

A 

A 

Position 1 

0.01 

-0.24 

0.09 

0.05 

0.05 


- 

2 

0.04 

-0.24 

0.21 

0.01 

-0.11 

- 

- 

3 

0.03 

-0.43 

-0.93 

0.03 

-0.27 

- 

- 

4 

0.03 

-0.46 

-0.75 

0.01 

-0.26 

- 

. - 

5 

0.04 

-0.37 

-0.54 

0.02 

-0.13 

- 

- 

6 

0.03 

-0.24 

0.22 

0.08 

-0.11 

- 

- 

7 

-0.02 

-0.30 

0.22 

0.14 

-0.10 

0.10 

0 

8 

-0.08 

0.63 

0.23 

0.26 

-0.10 

0.26 

0.57* 

9 

-0.40 

0.22 

0.23 

0.42 

-0.24 

0.50 

0.59 

10 

-0.67 

0.21 

0.23 

0.47 

-0.53 

0.60 

0.69 

11 

-0.50 

0.28 

0.25 

0.35 

-0.44 

0.69 

0.69 

12 

-0.13 

0.81 

0.24 

0.25 

-0.20 

0.40 

0.42* 

13 

-0.01 

-0.35 

0.20 

0.25 

0.19* 

0.33 

0.11* 

14 

0.01 

-0.01 

0.10 

0.26 

0.28 

0.54 

0.14 

15 

0.18 

0.07 

.0.05 

0.28 

0.36 

0.82 

0.36 

16 

0.20 

0.15 

0 

0.26 

0.36 

1.00 

0.46 

17 

0.18 

0.12 

0.05 

0.21 

0.30 

0.74 

0.36 

18 

0.11 

-0.10 

0.10 

0.24 

0.22 

0.54 

0.18 

19 

-0.01 

-0.35 

0.19 

0.26 

0.24 

0.33 

0.11* 

20 

-0.16 

0.96 

0.20 

0.33 

0.80 

0.40 

0.46* 

21 

-0.70 

0.37 

0.24 

0.63 

1.00 

0.60 

0.91 

22 

-1 .00 

0.36 

0.21 

1.00 

0.85 

0.60 

1 .00 

23 

-0.60 

0.39 

0.34 

0.79 

0.48 

0.50 

0.91 

24 

-0.14 

1.00 

0.19 

0.40 

0.48 

0.28 

0.50 

25 

-0.02 

-0.30 

0.26 

0.12 

0.31 

0.14 

0 

26 

0.05 

-0.07 

0.16 

0.10 

0.29 

- 

- 

27 

0.10 

-0.03 

-0.56 

0.07 

0.40 

- 

- 

28 

0.106 

-0.04 

-0.75 

0.07 

0.49 

- 

- 

29 

0.110 

-0.11 

-1.00 

0.08 

0.54 

- 

- 

30 

0.09 

-0.11 

0.34 

0.08 

0.58 

- 

- 

Position 31 

0.02 

-0.08 

0.20 

0.09 

0.14 

- 

- 


*90° Phase shift to reference. 
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TABLE 3 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-2-1D 

SPI-2-1D 

SPI-2-1D 

SPI-2-1D 

SPI-2-1D 

SPI-2-1D 

SPI-2-1D 

Frequency 

80 Hz 

88 Hz 

111 Hz 

128 Hz 

94 Hz 

108 Hz 

134 Hz 

Speaker 
Cond . 

A 

A 

A 

B 

A 

B 

B 

Position 1 

0 

0 

0.09 

-0.01 

- 

_ 

- 

2 

-0.01 

0.32 

0.24 

0.40 

- 

- 

- 

3 

-0.09 

0.77 

0.28 

0.58 

- 

- ' 

- 

4 

-0.09 

1.00 

0.29 

0.70 

- 

- 

- 

5 

-0.08 

0.82 

0.25 

0.48 

- 

- 

- 

6 

-0.01 

0.27 

0.21 

0.25 

- 

- 

- 

7 

-0.06 

0.30* 

0.29 

0.25 

-0.06 

0.05* 

0.22* 

8 

-0.14 

-0.30 

0.44 

0.38 

-0.14 

0.12* 

0.44* 

9 

-0.33 

-0.68 

0.73 

0.60 

-0.28 

0.28* 

0.78* 

10 

-0.40 

-0.30 

1.00 

0.70 

-0.40 

-0.44 

0.96* 

11 

-0.35 

-0.93 

0.77 

0.60 

-0.40 

-0.44 

0.85* 

12 

-0.16 

-0.24 

0.44 

0.38 

-0.20 

0.12* 

0.47* 

13 

-0.08 

-0.25 

0.43 

0.15 

0.07 

0.25 

0.25* 

14 

0.06 

0.17 

0.52 

0.01* 

0.40 

0.53 

0.31 

15 

0.13 

0.34 

0.63 

0.15* 

0.80 

1.00 

0.53 

16 

0.18 

0.27 

0.69 

0.15* 

1.00 

1 .00 

0.53 

17 

0.06 

0.15 

0.63 

0.21* 

0.80 

0.89 

0.47 

18 

0.06 

0.09 

0.49 

0.15* 

0.58 

0.50 

0.25 

19 

-0.06 

0.09 

0.44 

0.18* 

0.16 

0.40 

0.22* 

20 

-0.20 

0.16 

0.61 

0.45 

-0.40 

0.40 

0.44* 

21 

-0.62 

0.30 

0.94 

0.84 

-0.46 

0.60 

0.85* 

22 

-0.98 

0.43 

1.00 

0.84 

-0.46 

0.90 

1 .00* 

23 

-0.62 

0.35 

0.88 

0.72 

-0.40 

0.80 

0.78* 

24 

-0.20 

0.26 

0.55 

0.45 

-0.28 

0.50 

0.41* 

25 

-0.09 

- 

0.39 

0.30 

-0.10 

0.24 

0.28* 

26 

-0.02 

-0.32 

0.28 

0.42 

- 

- 

- 

27 

-0.47 

-0.54 

0.38 

0.95 

- 

- 

- 

28 

-0.59 

-0.63 

0.39 

1.00 

- 

- 

- 

29 

-0.59 

-0.58 

0.36 

0.92 

- 

- 

- 

30 

-0.16 

-0.12 

0.32 

0.60 ' 

- 

- 

- 

Position 31 

-0.02 

0.09 

0.15 

0.20 

- 

- 

- 


*90° Phase shift to reference. 
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TABLE 4 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-2-2 

SPI-2-2 

SPI-2-2 

SPI-2-2 

SPI-2-2 

SPI-2-2 

Frequency 

57 Hz 

61 Hz 

68 Hz 

51 Hz 

61 Hz 

67 Hz 

Speaker Cond. 

A 

A 

A 

A 

A 

A 

Position 1 

0.03 

0.05 

0.05 

- 

— 

- 

2 

0.08 

0.13 

0.10 

- 

- 

- 

3 

0.16 

0.15 

0.15 

- 

- 

- 

4 

0.16 

0.17 

0.18 

- 

- 

- 

5 

0.12 

0.15 

0.17 

- 

- 

- 

6 

0.04* 

0.09 

-0.10 

0.04* 

0.07 

-0.11 

7 

-0.18 

-0.23 

-0.37 

-0.11 

-0.17 

-0.38 

8 

-0.36 

-0.43 

-0.79 

-0.22 

-0.37 

-0.81 

9 

-0.40 

-0.63 

-1.00 

-0.36 

-0.54 

-1.00 

10 

-0.40 

-0.61 

-1.00 

-0.36 

-0.52 

-1.00 

11 

-0.24 

-0.30 

-0.47 

-0.22 

-0.29 

-0.56 

12 

-0.04 

0.07 

-0.13 

0.11 

-0.07 

-0.14 

13 

0.27 

0.29 

0.30 

0.49 

0.29 

0.31 

14 

0.63 

0.67 

0.67 

0.69 

0.67 

0.67 

15 

0.74 

0.92 

0.87 

0.95 

0.83 

0.93 

16 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

17 

0.91 

0.92 

0.93 

1.00 

0.83 

0.93 

18 

0.70 

0.69 

0.73 

0.78 

0.60 

0.67 

19 

0.40 

0.36 

0.33 

0.39 

0.40 

0.33 

20 

0.06 

0.07 

0.05 

0.09 

0.08* 

-0.10 

21 

-0.20 

-0.29 

-0.37 

-0.15 

-0.30 

-0.47 

22 

-0.36 

-0.43 

-0.63 

-0.28 

-0.50 

-0.87 

23 

-0.40 

-0.54 

-0.73 

-0.28 

-0.60 

-1.00 

24 

-0.36 

-0.46 

-0.60 

-0.19 

-0.54 

-0.87 

25 

-0.16 

-0.21 

-0.33 

-0.12 

-0.27 

-0.47 

26 

-0.06 

-0.09 

-0.11 

0.02 

-0.10 

-0.14 

27 

0.07 

0.17 

0.08 

- 

- 

- 

28 

0.18 

0.18 

0.13 

- 

- 

- 

29 

0.18 

0.17 

0.17 

- 

- 

- 

30 

0.08 

0.13 

0.12 

- 

- 

- 

Position 31 

0.05 

0.08 

0.05 


- 

- 


*90° Phase shift to reference. 
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TABLE 5 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-2-2D 

SPI-2-2D 

SPI-2-2D 

SPI-2-2D 

SPI-2-2D 

Frequency 

55 Hz 

59 Hz 

68 Hz 

96 Hz 

64 Hz 

Speaker Cond . 

A 

C 

D 

B 

A 

Position 1 

0.09 T 

0 

0.05 

0.04 

- 

2 

0.14 

0 

0.20 

-0.04 

- 

3 

0.19 

0 

0.35 

-0.08 

- 

4 

0.19 

0 

0.36 

-0.07 

- 

5 

0.13 

0 

0.16 

-0.04 

- 

6 

0.07 T 

0 

0.06 

0.14 

0 

7 

-0.13 

-0.21 

-0.14 

0.44 

0 

8 

-0.19 

-0.29 

-0.36 

0.77 

-0.06 

9 

-0.25 

-0.35 

-0.40 

1.00 

-0.11 

10 

-0.23 

-0.29 

-0.36 

0.77 

-0.17 

11 

-0.15 

-0.23 

-0.20 

0.44 

-0.17 

12 

0.06 T 

0 

-0.06 

0.22 

-0.11 

13 

-0.32 

0.38 

0.22 

-0.11 

0.07 

14 

-0.50 

0.44 

0.89 

-0.16 

0.27 

15 

-1.00 

0.50 

1 .00 

-0.18 

0.56 

16 

-1.00 

1.00 

1 .00 

0.16* 

0.89 

17 

-0.63 

0.50 

1.00 

0.22 

1.00 

18 

-0.45 

0.38 

0.70 

0.22 

0.82 

19 

-0.25 

0.21 

0.22 

0.18 

0.45 

20 

0.11 T 

0.12 

-0.14 

0.13* 

0.18 

21 

0.15 T 

-0.18 

-0.28 

0.16* 

0.11 

22 

0.20 T 

-0.29 

-0.60 

0.18* 

0.06* 

23 

0.23 T 

-0.29 

-0.80 

0.20* 

0.08* 

24 

0.20 T 

-0.29 

-0.70 

0.20* 

0.08* 

25 

0.10 T 

-0.15 

-0.24 

0.16* 

0.08* 

26 

0.06 T 

0 

-0.10 

0.06* 

0 

27 

0.10 T 

0.12 

0.10 

0.06 

- 

28 

0.13 T 

0.12 

0.20 

0.11 

- 

29 

0.13 T 

0.12 

0.20 

0.11 

- 

30 

0.10T 

0 

0.10 

0.07 

- 

Position 31 

0.06 T 

0 

0.04 

0.03 

- 


*90° Phase shift to reference. 

T response at twice frequency of reference. 
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TABLE 6 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-3-1 

SPI-3-1 

SPI-3-1 

SPI-3-1 

SPI-3-1 

SPI-3-2 

SPI-3-2 

Frequency 

80 Hz 

88 Hz 

107 Hz 

1 17 Hz 

101 Hz 

62 Hz 

68 Hz 

Speaker Cond. 

A 

B 

A 

B 

A 

A 

A 

Position 1 

0.00 

- 

- 

- 

- 

- 

- 

2 

0.10 

- 

- 

- 

- 

- 

- 

3 

0.20 

- 

- 

- 

- 

- 

- 

4 

0.30 

0.60 

1.00 

0.77 

- 

0.00 

- 

5 

0.38 

- 


- 

- 

- 

- 

6 

0.20 

- 

- 

- 

- 

- 

- 

7 

0.00 

- 

- 

- 

- 

- 

- 

8 

-0.40 

- 

- 

- 

- 

- 

- 

9 

-0.83 

- 

- 

- 

- 

0.20 

0.00 

10 

-1 .00 

-0.60 

0.28 

0.42 

0.40 

- 

- 

11 

-0.80 

- 

- 

- 

- 

- 

0.00 

12 

-0.35 

- 

- 

- 

- 

- 

- 

13 

-0.06 

- 

- 

- 

- 

- 

0.25 

14 

0.10 

- 

- 

- 

- 

- 

- 

15 

0.30 

- 

- 

- 

- 

- 

- 

16 

0.39 

-0.15 

-0.18 

0.00 

1 .00 

1 .00 

1.00 

17 

0.24 

- 

- 

- 

- 

- 

- 

18 

0.12 

- 

- 

- 

- 

- 

- 

19 

-0.15 

- 

- 

- 

- 

- 

0.20 

20 

-0.40 

- 

- 

- 

- 

- 

- 

21 

-0.70 

- 

- 

- 

- 

- 

0.00 

22 

-0.80 

1 .00 

0.30 

0.50 

0.28 

- 


23 

-0.68 

- 

- 

- 

- 

0.58 

0.00 

24 

-0.40 

- 

- 

- 

- 

- 

- 

25 

-0.03 

- 

- 

- 

- 

- 

- 

26 

0.00 

- 

- 

- 

- 

- 

- 

27 

0.10 

- 

- 

- 

- 

- 

- 

28 

0.18 

-0.32 

0.82 

1 .00 

- 

0.00 

- 

29 

0.10 

- 

- 

- 

- 

- 

- 

30 

0.00 

- 

- 

- 

- 

- 

- 

Position 31 

0.00 

- 

- 

- > 

— 

- 

- 
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TABLE 7 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-3-1D 

SP1-3-1D 

SPI-3-1D 

Frequency 

105 Hz 

105 Hz 

109 Hz 

Speaker 
Cond . 

A 

D 

B 

Position 1 

0.09 

0.15 

0.01 

2 

0.15 

0.18 

0.08 

3 

0.56 T 

0.28 

0.20 

4 

0.20 T 

0.25 

0.35 

5 

0.22 T 

0.23 

0.27 

6 

0.25 T 

0.20 

0.15 

7 

0.35 T 

0.43 

-0.04 

8 

0.74 

0.75 

-0.30 

9 

0.87 

1 .00 

1 .00 

10 

0.94 

1 .00 

0.59 

11 

1.00 

1 .00 

0.63 

12 

0.81 

0.85 

-0.43 

13 

0.56 

0.55 

-0.07 

14 

0.22 

0.30 

0.71 

15 

0.17 T 

0.15 

0.25 

16 

0.17T 

0.15 

0.25 

17 

0.17 T 

0.15 

0.25 

18 

0.22 

0.25 

0.11 

19 

0.47 

0.55 

-0.10 

20 

0.70 

0.70 

-0.33 

21 

0.88 

0.80 

- 

22 

0.92 

0.75 

- 

23 

0.87 

0.80 

0.20 

24 

0.67 

0.75 

0.11 

25 

0.15 

0.20 

0.20 

27 

0.10 T 

0.20 

0.32 

28 

0.10 T 

0.25 

0.38 

29 

0.10 T 

0.25 

0.28 

30 

0.10 T 

0.25 

0.13 

Position 31 

0.10 

0.15 

0.03 


SPI-3-1D SPI-3-1D SPI-3-1D SPI-3-1D 


127 Hz 

134 Hz 

106 Hz 

139 Hz 

A 

B 

A 

B 

0.13 

0.12 



0.36 

0.61 

- 

- 

0.31 

0.81 

- 

- 

0.27 

0.91 

- 

- 

0.30 

0.76 

- 

- 

0.36 

0.44 


- 

0.16 

0.30 

0.26 

0.09 

0.26 

0.44 

0.63 

0.25 

0.36 

0.36 

1 .00 

0.62 

0.38 

0.96 

1 .00 

0.50 

0.38 

0.96 

1 .00 

0.49 

0.31 

0.57 

0.69 

0.35 

0.27 

0.52 

0.47 

0.12 

-0.36 

-0.09 

0.40 

0.08 

-0.24 

-0.21 

0.60 

0.09 

-0.24 

-0.18 

0.74 

0.12 

-0.24 

-0.12 

0.60 

0.13 

-0.36 

-0.12 

0.33 

0.08 

0.50 

-0.27 

0.33 

-0.21 

0.37 

-0.51 

0.20 

-0.44 

0.67 

-0.91 

0.74 

-0.80 

0.74 

-1.00 

0.87 

-1 .00 

0.61 

-0.88 

0.67 

-0.78 

0.49 

-0.52 

0.33 

-0.35 

0.25 

-0.43 

- 

- 

1.00 

-0.64 

- 

- 

0.30 

-0.76 

- 

- 

0.43 

-0.67 

- 

- 

0.25 

-0.55 

- 

- 

0.10* ' 

-0.15 

- 

- 


T response at twice frequency of reference. 
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TABLE 8 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS: 
ONE-DIMENSIONAL PANELS 


Panel 

SPI-3-2D 

SPI-3-2D 

SPI-3-2D 

SPI-3-2D 

SPI-3-2D 

Frequency 

63 Hz 

74 Hz 

107 Hz 

65 Hz 

71 Hz 

Speaker Cond . 

A 

A 

B 

A 

D 

Position 1 

0.13 T 

0.10T 

-0.08 

- 

- 

2 

0.10 T 

0.16T 

-0.39 



3 

0.10 T 

0.16 T 

-0.69 

- 

- 

4 

0.10 T 

0.10 T 

-0.54 

- 

- 

5 

0.15 T 

-0.13 

-0.23 

- 

- 

6 

0.45 T 

-0.31 

0.19 

- 

- 

7 

0.50 * 

-0.57 

0.50 

0.07 

0.19T 

8 

1 .00 

-0.72 

0.85 

0.07 T 

-0.38 

9 

0.75 

-0.61 

1.00 

0.10 T 

-0.50 

10 

0.75 

-0.36 

1 .00 

-0.10 

-0.50 

11 

0.35 

-0.14 

0.42 

-0.10 

-0.30 

12 

0.12 

0.30 

0.23 

0.07 

0.09 T 

13 

-0.10 

0.68 

0.31 * 

0.25 

0.33 

14 

-0.20 

0.85 

0.39* 

0.72 

0.67 

15 

-0.25 

1.00 

0.39* 

1.00 

0.50 

16 

-0.25 

1.00 

0.62 

1.00 

1.00 

17 

-0.30 

0.82 

0.58 

1 .00 

0.50 

18 

-0.22 

0.33 

0.39 

0.80 

0.40 

19 

0.12 T 

0.27 

0.27 

0.42 

0.25 

20 

0.13 

-0.10 

0.31 

0.17 

0.08 T 

21 

0.36 

-0.31 

0.39* 

-0.17 

-0.29 

22 

0.55 

-0.63 

0.58 * 

-0.33 

-0.50 

23 

0.63 

-0.63 

0.69 * 

-0.33 

-0.63 

24 

0.55 

-0.56 

0.69 

-0.25 

-0.48 

25 

0.30 

-0.31 

0.39 * 

0.12 T 

-0.25 

26 

0.13 

-0.10 

0.08 * 

- 

- 

27 

0.05 T 

-0.06 

0.12 * 

- 

- 

28 

0.08 

-0.08 

0.23 * 

- 

- 

29 

0.10 

-0.10 

0.31 * 

- 

- 

30 

0.08 

0.06 

0.15 * 

- 

- 

Position 31 

0.03 

0.06 

0 • 

- 

- 


*90° Phase shift to reference. 

T response at twice frequency of reference. 
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TABLE 9 


NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, u in/in 
SPECIMEN SPI-1 


Configuration 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

3 Bay 

3 Bay 

3 Bay 

Frequency, Hz 

82 

106 

119 

126 

108 

115 

127 

Strain Gage 1 

0.00 

0.00 

0.00 

0.00 

- 

- 

- 

2 

1.00 

0.26 

-0.50 

0.40 

- 

- 

- 

3 

-0.87 

-0.40 

0.61 

-1.00 

- 

- 

- 

4 

-0.93 

-0.40 

0.56 

-0.73 

-0.56 

-0.43 

-0.20 

5 

0.00 

0.40 

0.28* 

0.73 

1 .00 

0.60 

0.38 

6 

-0.80 

-0.58 

0.84 

-0.77 

-0.83 

-0.50 

-0.60 

7 

-0.67 

-0.45 

0.56 

-0.77 

-0.83 

-0.50 

-0.56 

8 

0.67 

-0.40 

-1.00 

0.25 

0.00 

-1 .00 

0.40 

9 

-0.67 

-0.32 

0.70 

0.50 

-0.56 

-0.63 

0.56 

10 

-0.80 

-0.32 

0.70 

0.40 

-0.67 

-0.68 

0.46 

11 

0.00 

-0.84 

-0.42 

-1.00 

0.95 

0.50 

-1 .00 

12 

-0.67 

-0.47 

0.20 

0.94 

-0.67 

-0.38 

0.60 

13 

-0.67 

-0.40 

0.20 

0.94 

- 

- 

- 

14 

0.67 

1 .00 

-0.28 

-0.67 

- 

- 

- 

Strain Gage 15 

0.60 

-0.13 

0.00 

0.11T 

- 

- 

- 

e, p. in/in 

63 

158 

150 

250 

75 

168 

208 


*90° Phase shift to reference. 

T response at twice reference frequency. 
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TABLE 10 


NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, p in/in 
SPECIMEN SPI-2-1 


Configuration 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

3 Bay 

3 Bay 

Frequency, Hz 

88 

92 

103 

113 

120 

91 

109 

Strain Gage 1 

0.00 

0.00 

0.00 

0.00 

0.00 

- 

- 

2 

0.80 

0.28 

-0.40 

1 .00* 

0.70 

- 

- 

3 

-0.27 

-1 .00 

0.45 

-0.89 

-1 .00 

- 

- 

4 

-0.13 

-0.83 

0.30 

-0.72 

-1 .00 

0.54 

-0.36 

5 

-1.00 

1 .00 

1 .00 

0.54 

0.60 

-0.54 

0.36 

6 

-0.27 

-0.95 

-0.70 

0.79 

-0.80 

-0.77 

-0.29 

7 

-0.27 

-0.83 

-0.63 

-0.54 

-0.80 

-0.77 

-0.29 

8 

1.00 

0.27 

0.00 

0.79 

0.25 

0.77 

0.36 

9 

-0.40 

-0.45 

0.18 

-0.64 

0.20 

-0.93 

-0.72 

10 

-0.40 

-0.56 

0.20 

-0.64 

0.35 

-1.00 

-1 .00 

11 

-0.40 

0.25 

-0.40 

0.47 

-0.75 

-0.77 

0.89 

12 

-0.27 

-0.89 

0.15 

-0.54 

1.00 

0.62 

-0.72 

13 

-0.27 

-0.56 

-0.20 

-0.47 

0.75 

- 

- 

14 

0.54 

0.55 

0.60 

0.82 

-0.75 

- 

- 

Strain Gage 15 

0.00 

0.00 

0.00 

0.00 

0.00 

- 

- 

e, p in/in 

312 

75 

83 

117 

83 

54 

100 


*90° Phase shift to reference. 
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TABLE 11 


NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, u in/in 
SPECIMEN SPI-2-1D 


Configuration 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

3 Bay 

3 Bay 

3 Bay 

Frequency, H 

z 

80 

88 

111 

128 

94 

108 

134 

Strain Gage 

1 

-0.32 

0.00 

0.33 

0.00 

- 

- 

- 


2 

0.80 

0.73 

=0.67 

0.55 

- 

- 

- 


3 

-0.60 

0.00 

-0.33 

-0.55 

- 

- 

- 


4 

-0.52 

0.00 

-0.33 

-0.55 

0.50 

-0.56 

0.84 


5 

0.40 

-1 .00 

0.33 

0.91 

-0.57 

0.63 

-1.00 


6 

-0.32 

0.33 

-0.60 

-0.46 

0.33 

-0.63 

0.67 


7 

-0.28 

0.33T 

-0.74 

-0.46 

0.33 

-0.63 

0.56 


8 

0.40 

0.67 

1 .00 

0.00 

0.50* 

1.00 

0.56 


9 

-0.44 

-0.67 

-0.60 

0.55 

0.33 

-0.63 

-0.45 


10 

-0.60 

-0.67 

-0.47 

0.55 

0.67 

-0.63 

-0.67 


11 

1.00 

0.87 

0.40 

-1 .00 

-1 .00 

0.31 

1 .00 


12 

-0.60 

0.33 

-0.40 

0.82 

0.83 

-0.38 

-0.56 


13 

-0.56 

0.47 

-0.47 

0.82 


- 

- 


14 

0.28 

-0.80 

-0.67 

-0.73 

- 

- 

- 

Strain Gage 

15 

0.28 

0.47 

0.27 

0.00 

- 

- 

- 

e, i-l in/in 


105 

125 

312 

92 

125 

67 

75 


*90° Phase shift to reference. 

T Response at twice reference frequency. 
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TABLE 12 


NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, u in/in 
SPECIMEN SPI-2-2 


Configuration 

5 Bay 

5 Bay 

5 Bay 

3 Bay 

3 Bay 

3 Bay 

Frequency, Hz 

57 

61 

68 

61 

51 

67 

Strain Gage 1 

-0.13 

-0.27 

0.0 

- 

- 

- 

2 

0.20 

0.31 

0.17 

- 

- 

- 

3 

0.16 

0.31 

0.27 

- 

- 

- 

4 

0.27 

0.50 

0.37 

0.70 

0.33 

0.85 

5 

-0.55 

-0.69 

-0.50 

-O.89 

-0.33 

-1.0 

6 

-0.36 

-0.19 

0.17 

-0.30 

-0.40 

0.50 

7 

-0.65 

-0.54 

-0.20 

-0.50 

-0.80 

-0.25 

8 

1.00 

1 .00 

1 .00 

1 .0 

1.00 

1 .00 

9 

-0.69 

-0.58 

-0.27 

-0.60 

-0.83 

-0.35 

10 

-0.44 

-0.27 

0.17 

0.25* 

-0.47 

0.40 

11 

-0.44 

-0.54 

-0.37 

-0.65 

-0.33 

-0.75 

12 

0.25 

0.35 

0.27 

0.60 

0.40 

0.75 

13 

0.18 

0.35 

0.20 

- 

- 

- 

14 

0.20 

0.23 

0.13 

- 

- 

- 

Strain Gage 15 

-0.09 

0.15 

0.0 

- 

- 

- 

e, n in/in 

230 

108 

125 

83 

125 

167 


*90° Phase shift to reference. 



TABLE 13 


NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, p in/in 
SPECIMEN SPI-2-2-D 


Configuration 

5 Bay 

5 Bay 

5 Bay 

3 Bay 

Frequency, Hz 

55 

68 

96 

64 

Strain Gage 1 

0.5 

0.0 

0.0 

- 

2 

0.7 

0.0 

0.0 

- 

3 

-0.5 

0.0 

0.40 

- 

4 

-0.5 

0.0 

0.55 

0.0 

5 

0.7 

-0.83 

-0.75 

-0.6 

6 

-0.5 

-1 .0 

0.40 

-1.0 

7 

-0.5 

-1.0 

0.40 

-1 .0 

8 

1.0 

1.0 

0.50 

-0.8 

9 

-0.6 

-0.83 

-0.75 

-1.0 

10 

-0.8 

-0.83 

-0.75 

-1.0 

11 

1.0 

-0.83 

1.00 

-1.0 

12 

-0.5 

0.0 

-0.75 

0.5 

13 

-0.5 

0.0 

-0.50 

- 

14 

0.0 

0.0 

0.0 

- 

Strain Gage 15 

0.0 

0.0 

0.0 

- 

e, p in/in 

45 

25 

83 

42 
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TABLE 14 


NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e= NORMALIZING STRAIN, u in/in 
SPECIMEN SPI-3-1 


Configuration 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

3 Bay 

3 Bay 

Frequency, 1 

Hz 

80 

88 

92-93 

107 

117 

90 

101 

Strain Gage 

1 

0.13 

0.0 

-0.38 

-0.25 

-0.19 

- 

- 


2 

0.36* 

0.35 

0.94 

1.00 

0.75 

- 

- 


3 

0.36* 

0.0 

-1.00 

-0.45 

-0.93 

- 

- 


4 

-0.39 

0.25 

-0.94 

-0.42 

-0.93 

-0.83 

-0.72 


5 

0.71 

-0.75 

0.56 

0.50 

1.00 

1.00 

1.00 


6 

0.21* 

0.50 

-0.31 

0.33 

-0.55 

-0.63 

-0.43 


7 

0.27* 

0.40 

-0.31 

0.50 

-0.44 

-0.21 

0.43 


8 

-0.71 

0.50 

0.0 

-1.00 

0.15 

-0.42 

-0.86 


9 

0.21* 

-0.75 

-0.47 

0.67 

0.37 

0.0 

0.50 


10 

-0.64 

-1.00 

-0.63 

0.50 

0.37 

-0.42 

-0.36 


11 

1.00 

1.00 

1.00 

0.33 

-0.74 

1 .00 

0.79 


12 

-0.50 

-0.50 

-0.78 

-0.67 

0.74 

-0.63 

-0.64 


13 

0.27* 

0.75 

-0.78 

-0.67 

0.74 

- • 

- 


14 

-0.72 

-1 .00 

0.47 

0.83 

-0.55 

- 

- 

Strain Gage 

15 

0.30 

0.50 

0.0 

0.0 

0.0 

- 

- 

e, (i in /in 


232 

83 

133 

250 

225 

100 

290 


*90° Phase shift to reference. 
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TABLE 15 


NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, H In/in 
SPECIMEN SPI-3-2 


Configuration 

5 Bay 

5 Bay 

3 Bay 

3 Bay 

3 Bay 

Frequency, Hz 

54 

62 

62 

68 

101 

Strain Gage 1 

0.0 

0.0 

- 

- 

- 

2 

0.0 

0.0 

- 

- 

- 

3 

0.0 

0.0 

- 

- 

- 

4 

0.0 

0.13 

0.0 

0.18 

0.61 

5 

0.0 

-0.50 

0.0 

-0.21 

-0.89 

6 

-0.62 

-0.30 

-0.42 

-0.36 

0.33 

7 

-1.00 

-0.85 

-0.84 

-0.89 

-0.89 

8 

1.00 

0.85 

1.00 

1.00 

1.00 

9 

-0.93 

-0.50 

-0.88 

-0.89 

-1 .00 

10 

-0.54 

-0.40 

-0.50 

-0.36 

0.28 

11 

0.0 

-1.00 

-0.21 

-0.27 

-0.56 

12 

0.0 

0.50 

0.0 

0.18 

0.45 

13 

0.0 

0.65 

- 

■ - 

- 

14 

0.0 

0.30 

- 

- 

- 

Strain Gage 15 

0.0 

0.0 

- 

- 

- 

e, |a in/in 

54 

167 

100 

233 

75 
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TABLE 16 

NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, p in/in 
SPECIMENS SPI-3-1D AND SPI-3-2D 

SPI-3-1D SPI-3-2D 


Configuration 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

5 Bay 

3 Bay 

3 Bay 

Frequency, Hz 

105 

115 

127 

74 

107 

65 

71 

Strain Gage 1 

0.18* 

0.17* 

0.25 

0.33 

- 

- 

- 

2 

0.42 

0.83 

-0.63 

-0.47 

- 

- 

- 

3 

-0.58 

-0.63 

0.42* 

-0.60 

-0.60 

- 

- 

4 

-0.79 

-1.00 

0.42* 

-0.67 

-0.80 

0.50 

0.50 

5 

1.00 

0.83 

0.42 

1 .00 

0.67 

-0.25 

-1 .00 

6 

-0.66 

-0.42 

-0.63 

0.67 

-0.87 

-0.90 

1.00T 

7 

-0.53 

-0.33 

-0.63 

0.53 

-0.67 

0.0 

0.60 

8 

0.50 

-0.17 

1 .00 

-1 .00 

0.67 

1.00 

0.90 

9 

-0.58 

-0.63 

-0.84 

0.53 

-0.67 

-0.65 

-0.50 

10 

-0.63 

0.58 

-0.75 

0.67 

-0.80 

-0.90 

-0.50 

11 

0.79 

-0.75 

0.57 

0.87 

1.00 

0.40 

0.40T 

12 

-0.66 

0.83 

-0.50 

-0.80 

-0.53 

0.0 

0.40 

13 

-0.53 

0.79 

0.42* 

-0.67 

-0.47 

- 

— 

14 

0.58 

-0.63 

-0.75 

-0.47 

-0.80 

- 


Strain Gage 15 

-0.36 

0.33 

0.42 

0.47* 

0.40 

- 

- 

e, p in/in 

158 

200 

100 

63 

63 

83 

83 


*90° Phase shift to reference. 

T Response at twice frequency of reference. 
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TABLE 17 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS 
NINE-BAY MACHINED PANEL SPECIMEN 
(See Figure 6) 


Frequency, Hz 

88 

94 

148 

175 

188 

Position 1 

0.008 

0.004 

0.014 

0.021 

0.028 

2 

0.005 

0.006 

0.030 

0.050 

0.065 

'3 

0.015 

0.020 

0.052 

0.067 

0.102 

4 

0.028 

0.038 

0.074 

0.100 

0.185 

5 

0.045 

0.064 

0.100 

0.133 

0.296 

6 

0.055 

0.150 

0.104 

0.150 

0.296 

7 

0.070 

0.180 

0.111 

0.150 

0.352 

8 

0.088 

c 0.200 

0.111 

0.133 

0.407 

9 

0.125 

0.170 

0.100 

0.217 

0.444 

10 

0.163 

0.280 

0.100 

0.183 

0.537 

11 

0.188 

0.310 

0.082 

0.217 

0.537 

12 

0.219 

0.270 

0.067 

0.250 

0.556 

13 

0.263 

0.370 

0.044 

0.233 

0.593 

14 

0.288 

0.410 

0.022* 

0.333 

0.593 

15 

0.300 

0.420 

-0.030 

0.383 

0.593 

16 

0.300 

0.420 

-0.044 

0.433 

0.593 

17 

0.275 

0.390 

-0.044 

0.433 

0.593 

18 

0.238 

0.360 

-0.044 

0.433 

0.574 

19 

0.200 

0.320 

-0.041 

0.433 

0.556 

20 

0.163 

0.270 

-0.032 

0.399 

0.482 

21 

0.125 

0.230 

-0.044 

0.367 

0.463 

22 

0.100 

0.170 

-0.026 

0.367 

0.407 

23 

0.075 

0.120 

-0.024 

0.317 

0.389 

24 

0.063 

0.092 

-0.022 

0.283 

0.352 

25 

0.050 

0.068 

-0.020 

0.233 

0.278 

26 

0.025 

0.040 

-0.018 

0.158 

0.204 

27 

0.013 

0.024 

-0.013 

0.108 

0.120 

28 

0.005 

0.008 

-0.007 

0.049 

0.065 

Position 29 

0.008 

0,004 

- 

0.037 

0.028 


*90° Phase shift to reference. 
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TABLE 17 
(CONTINUED) 


Frequency, Hz 

88 

94 

148 

175 

188 

Position 30 

0.008 

0.060* 

0.011 

0.167 

-0.092 

31 

0.050 

0.080 

0.556 

-0.108 

-0.278 

32 

0.088 

0.096 

0.815 

-0.108 

-0.130 

33 

0.125 

0.128 

0.593 

-0.733 

-0.074 

34 

0.150 

0.144 

0.593 

-0.733 

0.093* 

35 

0.163 

0.160 

0.815 

-0.750 

0.148 

36 

0.163 

0.160 

0.704 

0.517 

0.315 

37 

0.150 

c 0.160 

0.269 

0.383 

0.482 

38 

0.150 

0.220 

0.133 

0.267 

0.482 

39 

0.225 

0.290 

0.222 

0.517* 

0.556 

40 

0.375 

0.440 

0.315 

-0.800 

0.407 

41 

0.500 

0.600 

0.296 

-0.450 

0.259 

42 

0.725 

0.740 

0.052* 

-0.450 

0.148 

43 

0.925 

0.860 

-0.185 

-0.617 

-0.667 

44 

1.000 

1 .000 

-0.241 

-0.750 

-1 .000 

45 

0.950 

0.900 

-0.370 

0.750* 

-0.482 

46 

0.800 

0.900 

-0.704 

0.300 

-0.259 

47 

0.550 

0.780 

-1 .000 

0.500 

0.185* 

48 

0.300 

0.560 

-0.324 

0.708 

0.482 

49 

0.175 

0.300 

-0.185 

1 .000 

0.704 

50 

0.158 

0.260 

-0.048 

0.917 

0.482 

51 

0.113 

0.200 

0.030 

0.583 

0.556 

52 

0.138 

0.200 

0.048 

0.250 

0.519 

53 

0.125 

0.200 

0.048 

0.083* 

0.250 

54 

0.113 

0.160 

0.052 

-0.100 

0.093 

55 

0.088 

0.100 

0.052 

-0.217 

-0.185 

56 

0.063 

0.080 

0.052 

-0.333 

-0.463 

57 

0.038 

0.080* 

0.059 

-0.333 

-0.648 

Position 58 

0.008 

0.020* 

0.044 

-0.125 

-0.222 


*90° Phase shift to reference. 


TABLE 17 
(CONTINUED) 


Frequency, Hz 

88 

94 

148 

175 

188 

Position 59 

0.010 

0.004 

0.013 

0.021 

0.028 

60 

0.005 

0.006 

0.037 

0.067 

0.074 

61 

0.008 

0.018 

0.059 

0.083 

0.111 

62 

0.023 

0.040 

0.085 

0.100 

0.185 

63 

0.033 

0.060 

0.104 

0.133 

0.241 

64 

0.050 

0.100 

0.119 

0.167 

0.296 

65 

0.065 

0.140 

0.126 

0.167 

0.333 

66 

0.088 

c 0.180 

0.133 

0.150 

0.370 

67 

0.125 

0.200 

0.133 

0.250 

0.444 

68 

0.163 

0.300 

0.133 

0.200 

0.482 

69 

0.200 

0.340 

0.133 

0.217 

0.519 

70 

0.238 

0.380 

0.126 

0.225 

0.519 

71 

0.275 

0.430 

0.096 

0.225 

0.574 

72 

0.300 

0.440 

0.082 

0.225 

0.611 

73 

0.313 

0.460 

0.044 

0.267 

0.630 

74 

0.300 

0.460 

0.022 

0.283 

0.630 

75 

0.275 

0.420 

-0.030 

0.299 

0.630 

76 

0.250 

0.400 

-0.030 

0.317 

0.593 

77 

0.188 

0.350 

-0.030 

0.333 

0.574 

78 

0.163 

0.300 

-0.028 

0.333 

0.555 

79 

0.125 

0.220 

-0.037 

0.267 

0.463 

80 

0.100 

0.180 

-0.026 

0.367 

0.444 

81 

0.075 

' 0.120 

-0.022 

0.333 

0.407 

82 

0.063 

0.086 

-0.022 

0.317 

0.370 

83 

0.050 

0.064 

-0.020 

0.300 

0.333 

84 

0.025 

0.040 

-0.016 

0.200 

0.250 

85 

0.013 

0.020 

-0.011 

0.138 

0.120 

86 

0.005 

0.008 

-0.007 

0.071 

0.111 

Position 87 

0.008 

0.004 

- 

0.042 

0.046 
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TABLE 17 
(CONTINUED) 


Frequency, Hz 

88 

94 

148 

175 

188 

Position 88 

0.003 

0.014 

0.007 

0.033 

0.080 

89 

0.013 

0.038 

0.019 

0.054 

0.135 

90 

0.045 

0.080 

0.032 

0.071 

0.185 

91 

0.063 

0.110 

0.048 

0.104 

0.259 

92 

0.088 

0.140 

0.061 

0.150 

0.324 

93 

0.108 

0.180 

0.089 

0.200 

0.370 

9 

0.125 

0.190 

0.100 

0.217 

0.444 

94 

0.138 

0.200 

0.119 

0.250 

0.444 

95 

0.145 

0.220 

0.122 

0.267 

0.463 

38 

0.150 

0.220 

0.133 

0.267 

0.482 

96 

0.150 

0.220 

0.141 

0.267 

0.482 

97 

0.150 

0.200 

0.141 

0.250 

0.444 

67 

0.125 

0.200 

0.133 

0.250 

0.444 

98 

0.100 

0.190 

0.111 

0.200 

0.370 

99 

0.095 

0.150 

0.104 

0.183 

0.296 

100 

0.063 

0.110 

0.082 

0.150 

0.222 

101 

0.1575 

0.070 

0.048 

0.100 

0.185 

102 

0.025 

0.036 

0.032 

0.058 

0.148 

103 

0.003 

0.016 

0.017 

0.042 

0.074 

104 

0.069 

0.112* 

-0.082 

-0.050 

-0.185 

105 

0.250 

0.030 

-0.163 

-0.200 

-0.426 

106 

0.438 

0.480 

-0.178 

-0.250 

-0.370 

107 

0.500 

0.560 

-0.178 

-0.283 

-0.241 

108 

0.425 

0.500 

-0.148 

0.217* 

0.148* 

109 

0.313 

0.380 

-0.089 

0.200 

0.630 

15 

0.300 

0.420 

-0.030 

0.433 

0.593 

110 

0.500 

0.520 

-0.010 

-0.383 

- 

111 

1.000 

0.780 

-0.222 

-0.667 

-0.741 

44 

1.000 

1.000 

-0.241 

-0.750 

-1 .000 

112 

0.925 

0.780 

-0.222 

-0.683 

-0.704 

113 

0.500 

0.520 

-0.089 

-0.150 

0.093 

73 

0.625 

0.460 

0.044 

0.283 

0.630 

114 

0.375 

0.400 

0.704 

0.400 

0.556 

115 

0.50 

0.520 

0.370 

0.217 

0.185 

116 

0.625 

0.620 

0.444 

0.133* 

-0.185 

117 

0.525 

0.520 

0.444 

0.167* 

-0.370 

118 

0.325 

0.520 

0.741 

-0.167 

-0.407 

Position 1 19 

0.125 

0.100 

0.185 

0.042* 

-0.167 


*90° Phase shift to reference. 
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TABLE 17 
(CONTINUED) 


Frequency, Hz 

88 

94 

148 

175 

188 

Position 120 

0.003 

0.016 

-0.007 

0.050 

0.074 

121 

0.013 

0.042 

-0.009 

0.104 

0.111 

122 

0.045 

0.080 

-0.015 

0.129 

0.167 

123 

0.063 

0.120 

-0.022 

0.183 

0.222 

124 

0.088 

0.160 

-0.011 

0.275 

0.296 

125 

0.108 

0.190 

-0.037 

0.308 

0.333 

21 

0.125 

0.230 

-0.044 

0.317 

0.407 

126 

0.138 

0.250 

-0.044 

0.333 

0.444 

127 

0.145 

0.260 

-0.044 

0.350 

0.462 

50 

0.158 

0.260 

-0.048 

0.367 

0.482 

128 

0.158 

0.260 

-0.044 

0.350 

0.462 

129 

0.150 

0.250 

-0.037 

0.267 

0.444 

79 

0.125 

0.220 

-0.037 

0.267 

0.407 

130 

0.108 

0.200 

-0.030 

0.183* 

0.370 

131 

0.083 

0.160 

-0.026 

0.167* 

0.333 

132 

0.063 

0.130 

-0.015 

0.117* 

0.259 

133 

0.045 

0.090 

-0.013 

0.100 

0.185 

134 

0.025 

0.048 

-0.005 

0.067 

0.158 

Position 135 

0.003 

0.020 

-0.007 

0.042 

0.074 


*90° Phase shift to reference. 
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TABLE 18 


STRAIN GAGE LOCATIONS 
NINE-BAY PANEL SPECIMENS 


Machines Panel Specimen 
(See Figure 7) 


SG 

(x, y) 

SG 

(x, y) 

XI 

( 3.25,18.0) 

Y1 

(13.0,12.35) 

X2 

( 6.40,18.0) 

Y2 

(13.0,16.25) 

X3 

( 9.75,18.0) 

Y3 

(13.0,18.3 ) 

X4 

(10.25,18.0) 

Y4 

(13.0,19.85) 

X5 

(13.00,18.0) 

Y5 

(13.0,21.85) 

X6 

(14.25,18.0) 

Y6 

(13.0,23.75) 

X7 

(15.75,18.0) 

Y7 

(13.0,24.25) 

X8 

(16.25,18.0) 

Y8 

(13.0,27.10) 

X9 

X10 

SG 

(19.60.18.0) 

(22.80.18.0) 

SPII-1 

(See Figure 8) 

(x,y) 

SG 

(x, y) 

XI 

( 6.0 ,18.5) 

Y1 

(13.5, 7.50) 

X2 

( 8.0 ,18.5) 

Y2 

(13.5,10.50) 

X3 

( 9.0 ,18.5) 

Y3 

(13.5,11.50) 

X4 

(10.75,18.5) 

Y4 

(13.5,14.90) 

X5 

(13.5 ,18.5) 

Y5 

(13.5,18.50) 

X6 

(15.75,18.5) 

Y6 

(13.5,21.90) 

X7 

(18.0 ,18.5) 

Y7 

(13.5,25.50) 

X8 

(19.0 ,18.5) 

Y8 

(13.5,26.50) 

X9 

(21.0 ,18.5) 

Y9 

(13.5,30.00) 

SG 

SPIt-2 

(See Figure 9) 

(*, y) 

SG 

(x, y) 

XI 

( 6.7 ,18.5) 

Y1 

(13.5, 8.5) 

X2 

( 8.2 ,18.5) 

Y2 

(13.5,10.75) 

X3 

( 9.7 ,18.5) 

Y3 

(13.5,13.00) 

X4 

(10.7 ,18.5) 

Y4 

(13.5,14.10) 

X5 

(13.5 ,18.5) 

Y5 

(13.5,18.5 ) 

X6 

(16.3 ,18.5) 

Y6 

(13.5,23.0 ) 

X7 

(17.2 ,18.5) 

Y 7 

(13.5,24.0 ) 

X8 

(18.6 ,18.5) 

Y8 

(13.5,26.25) 

X9 

(20.1 ,18.5) 

Y9 

(13.5,28.50) 
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TABLE 19 


NORMALIZED ACCELERATION (DISPLACEMENT) MEASUREMENTS 
NINE-BAY PANEL SPECIMENS SPII-1 AND SPII-2 


Frequency, Hz 

90 

97 

Specimen SPII-1 
101 107 112 

134 144 

168 

Bay Number 1 

- 

0.20 

0.23 

-0.23 

-0.38 

-0.40 

0.70 

2 

-0.25 

0.33* 

0.30* 

-0.27 

-0.62 

1.00 -0.27 

0.50 

3 

- 

0.20 

0.30 

-0.19 

-0.38 

0.27 

1 .00 

4 

-1.00 

-0.27 

0.33 

-0.23 

-0.38 

0.25 0.67 

0.65 

5 

0.19 

1.00 

1.00 

1 .00 

1.00 

0.70 

0.85 

6 

-0.25 

-0.20 - 

-0.27 

-0.31 

-0.54 

-0.65 0.27 

0.45 

7 

- 

0.23 

0.30 

0.23 

0.38 

-1 .00 

-1 .00 

8 

-0.85 

-0.23 

0.30 

0.27 

0.31 

-0.30 -0.33 

-1.00 

Bay Number 9 

- 

0.20 

0.27 

0.15 

0.38 

0.40 

-0.70 

Frequency, Hz 

74 

77 

Specimen SPII 
82 

-2 

no 

112 

126 

Bay Number 1 

- 

- 


- 

- 

0.15T 

- 

2 

- 

- 


-0.57 

-1.00 

-1.00 

0.86 

3 

- 

- 


- 

- 

- 

- 

4 

- 

- 


-0.57 

-0.19 

- 

- 

5 

1.00 

1 .00 

1.00 

0.27 

0.35 

1 .00 

6 

- 

- 


-0.71 

-0.27 

- 

- 

7 

- 

- 


- 

- 

- 

- 

8 

- 

0.20 

-0.74 

-0.27 

- 

-0.37 

Bay Number 9 

- 

- 


- 

- 

- 

- 


Note: Measurements taken at center of each panel bay. 
*90° Phase shift to reference. 
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TABLE 20 

NORMALlZED_HALF AMPLITUDE MODAL STRAIN DISTRIBUTION 
e ^NORMALIZING STRAIN, n in/in 
NINE BAY MACHINED PANEL SPECIMEN 

(See Table 18 and Figure 7) 


Stra in 


Frequency 

, Hz/S pea 

ker Phase Condition 


Gage 

88/A 

94/A 

123/A 

148/A 

170/A 

188/A 

XI 

-0.79 

-0.80 

0.09 

0.76 

0.00 

-0.40 

X2 

0.48 

0.40 

-0.09 

0.00 

0.00 

0.32 

X3 

-0.79 

-0.80 

0.09 

0.51 

0.00 

-0.40 

X4 

-1 .00 

-1 .00 

-0.86 

-0.51 

-0.61 

-1.00 

X5 

0.79 

0.80 

0.86 

0.51 

0.72 

1.00 

X6 

0.31 

0.36 

0.57 

0.25 

0.40 

0.00 

X7 

0.79 

-0.80 

-0.86 

-0.51 

1.00 

-1.00 

X8 

0.00 

-0.60 

0.57 

-0.25 

-0.40 

-0.40 

X9 

-0.63 

0.50 

-0.57 

0.25 

0.28 

0.28 

X10 

-0.79 

-0.60 

1.00 

-0.51 

0.00 

-0.61 

Y1 

-0.48 

-0.50 

-0.57 

-1.00 

0.80 

-0.61 

Y2 

0.22 

0.10 

0.29 

0.51 

-0.40 

0.61 

Y3 

0.16 

0.30 

0.46 

0.25 

0.61 

0.20 

Y4 

0.00 

0.00 

0.29 

-0.25 

0.40 

-0.20 

Y5 

0.00 

-0.10 

0.00 

-0.25 

0.00 

-0.20 

Y6 

-0.16 

-0.40 

-0.43 

0.25 

-0.61 

0.00 

Y 7 

0.00 

-0.10 

0.00 

1.00 

-0.21 

0.00 

Y8 

0.00 

0.10 

0.00 

-0.51 

0.00 

0.00 

e 

132 

208 

146 

83 

104 

104 



TABLE 21 

NORMALIZED_HALF AMPLITUDE MODAL STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, n in/in 
NINE BAY PANEL SPECIMEN SPII-1 

(See Table 18 and Figure 8) 


Strain 


Frequency 

, Hz/S pea 

ker Phase Condition 


Gage 

90/A 

97/A 

112/6 

134/6 

144/6 

168/6 

XI 

0.14 

0.16 

0.00 

-0.51 

0.84 

0.79 

X2 

0.12 

-0.16 

-0.17 

-0.35 

0.66 

0.33 

X3 

-0.24 

-0.37 

-0.34 

0.35 

-1 .00 

-1 .00 

X4 

-0.26 

-0.95 

-0.86 

-0.76 

-0.42 

-1 .00 

X5 

0.22 

0.95 

0.74 

0.76 

0.00 

0.67 

X6 

-0.38 

-1 .00 

-1 .00 

-1 .00 

-0.34 

-0.92 

X7 

-0.76 

-0.39 

-0.43 

-0.60 

0.00 

-0.46 

X8 

0.28 

-0.47 

-0.17 

-0.25 

0.00 

-0.33 

X9 

0.86 

-0.37 

0.17 

0.35 

0.00 

0.33 

Y1 

0.34 

0.27 

0.00 

0.00 

-0.42 

0.21 

Y2 

0.65 

-0.40 

0.20 

0.00 

-0.58 

0.46 

Y3 

-1 .00 

0.63 

-0.34 

-0.25 

0.76 

-0.40 

Y4 

-0.17 

-0.68 

-0.83 

-0.60 

0.34 

-0.52 

Y5 

0.12 

0.58 

0.46 

0.40 

0.00 

0.46 

Y6 

-0.09 

-0.53 

-0.66 

-0.40 

0.00 

-0.33 

Y7 

-0.14 

0.00 

-0.34 

-0.35 

-0.58 

0.67 

Y8 

0.14 

0.00 

0.34 

-0.25 

0.42 

-0.46 

Y9 

0.14 

0.00 

0.17 

-0.25 

0.42 

-0.46 

e 

242 

158 

125 

83 

50 

63 
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TABLE 22 

NORMALIZED HALF AMPLITUDE STRAIN DISTRIBUTION 
e = NORMALIZING STRAIN, \x in/in 
NINE BAY PANEL SPECIMEN SPII-2 

' (See Table 18 and Figure 9) 


Strain 



Frequency, 

H z/Speaker Phase 

Condition 



Gage 

74/A 

77/B 

77/C 

81/D 

82/A 

110/A 

11 2/D 

126/B 

XI 

0.00 

0.00 

0.00 

0.00 

-0.25 

0.50 

0.00 

0.00 

X2 

-0.50 

-0.52 

-0.50 

0.00 

-0.39 

0.00 

0.00 

0.00 

X3 

-1.00 

-0.76 

-0.79 

-0.72 

-0.75 

-0.69 

0.00 

-0.46 

X4 

0.50 

0.39 

0.50 

0.72T 

0.45 

0.40 

0.64 

0.46 

X5 

1 .00 

1 .00 

1 .00 

0.72 

1.00 

1 .Q0 

1.00 

0.46 

X6 

0.50 

0.52 

0.50 

0.72T 

0.50 

0.48 

0.64 

0.33 

X7 

-1.00 

-1.00 

-0.79 

-1 .00 

-1.00 

-0.60 

-0.52 

-0.33 

X8 

0.57 

-0.64 

-0.50 

.0.00 

-0.39 

0.00 

0.00 

0.00 

X9 

0.00 

0.00 

0.00 

0.00 

-0.42 

-0.50 

-0.52 

0.00 

Y1 

0.00 

0.00 

0.00 

0.00 

-0.39 

-0.36 

-0.52 

1.00 

Y2 

0.00 

0.00 

0.00 

0.00 

0.25 

0.36 

0.52 

-0.86 

Y3 

-0.72 

-0.76 

-0.60 

0.00 

-0.60 

0.00 

0.64 

0.46 

Y4 

0.50 

0.00 

0.50. 

0.00 

0.35 

0.50 

0.64 

-0.67 

Y5 

0.50 

0.76 

0.60 

0.72 

0.50 

0.50 

0.64 

0.46 

Y6 

0.36 

0.00 

0.40 

0.00 

0.20 

-0.50 

-0.64 

0.67 

Y7 

-1.00 

-0.76 

-0.60 

0.00 

-0.55 

0.50 

-0.64 

-0.67 

Y8 

0.00 

0.00 

-0.50 

0.00 

0.20 

0.50 

0.64 

0.33 

Y9 

0.00 

0.00 

-0.69 

0.00 

-0.25 

-0.69 

-0.76 

-0.46 

e 

58 

33 

42 

29 

84 

42 

33 

63 



TABLE 23 


Specimen 

SPI-1 

SPI-1 

SPI-1 

SPI-1 

SPI-1 

SPI-2-1 

SPI-2-1 

SPI-2-1 

SPI-2-1 

SPI-2-1 D 
SP I —2—1 D 

SPI-2-2 

SPI-2-2 

SPI-2-2 

SPI-2-2 

SP I —2—1 D 

SPI-3-1 

SPI-3-1 

SPI-3-1 

SPI-3-1 

SPI-3-1 

SPI-3-2 
SP 1-3— 2 

SPI-3-2D 

SPI-3-2D 


DAMPING RATIOS (PERCENT OF CRITICAL DAMPING) 
ONE-DIMENSIONAL PANELS 


Number of 
Bays 

Frequency 

Hz 

Strain 
Gage No. 

Damping Ratio 
(percent) 

5 

82 

X10 

1.1 

5 

119 

Xll 

1 .3 

5 

126 

X10 

1 .3 

3 

108 

X5/X11 

3. 0/1. 6 

3 

127 

Xll 

2.0 

5 

88 

X8 

1 .6 

5 

103 

X5 

2.7 

5 

120 

X5 

2.0 

3 

91 

X8 

3.0 

5 

88 

X5 

1.5 

5 

111 

X8 

2.0 

5 

57 

X8 

2.0 

5 

61 

X8 

1 .5 

3 

51 

X8 

2.7 

3 

61 

X8 

1.5 

5 

68 

X8 

1.5 

5 

80 

X5 

1.0 

5 

88 

X5 

2.9 

5 

107 

X8 

2.0 

5 

117 

X5 

1.0 

3 

101 

X8 

2.0 

5 

62 

Xll 

1.8 

3 

68 

X8 

1.6 

5 

74 

X8 

1.6 

3 

71 

X8 

1.5 
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APPENDIX A 


STIFFNESS AND CONSISTENT MASS MATRIX FOR 
A THIN-WALLED OPEN-SECTION BEAM 


The stiffness and consistent mass matrices presented here are in the form of a composite 
array with logic numbers (&x,&y) defining terms relating to the orientation of the beam ele- 
ment. Subscripts 1 and 2 in the loading and displacement vectors refer to ends 1 and 2 on 
the element. End 2 is always in the positive direction of the element axis from end 1 . Sub- 
scripts x and y in the expressions in the matrices refer to the x-axis and the y-axis, respec- 
tively. The cross-section nomenclature and coordinate directions are defined in figures 3 
and 4 of reference 1 . 


For stiffener elements conforming to the edge rotations of the plate element described 
in Appendix B, the loading and displacement coordinates at station i of the element are 
defined as 


- 

p 


r 

shear in the z direction 

M 


bending moment about the x-axis 

X 

« 

. = - 


M 

y 


bending moment about the y-axis 

M 

x yJ 

i 

twisting moment 

» - 

• 

d 


• > 
displacement in the z direction 

8 

X 

► = « 

rotation about the x-axis 

► 

9 

y 


rotation about the y-axis 

9 

x yJ 

i 

twist 

• - 


The composite beam stiffness matrix has the form 

CP }, =[K 1] ]{d} 1 +[K 12 ]{d} 2 
CP } 2 = [K 12 ] T {d} 1 + [K 22 ]{d } 2 
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Appendix A 


where for ribs parallel to the x-axis (a =1; a =0) and for rib parallel to the y 

(a =0;a =1) X y 

x y 


[K..] = a [K..] +a [K..] 

|J X Ij X y IJ y 


[K n ] 

I I X 


[K 12 ] x 


[k 99 ] 

ll x 


[K n ] 

11 y 


128 

-12R 

-68 L 

-6r L 

yy 

y 

yy x 

y x 


12Ti 

6R L 

y x 

48 L 2 

6T 2 
4r L 2 

(symmetric) 


yy x 

y 

4Y 3 . 

-128 

12R 

- 68 L 

-6r L 

yy 

y 

yy x 

y 

12R 

y 

- 12t i 

6R L 
y x 

2 

6T 2 
o .2 

68 L v 

- 6R L 

28 L 

2r L 

yy x 

y x 

yy x 

y x 

6r L 

- 6t 

2r L 2 

2Y 4 

L y x 

l 

y 

4 J 

' 128 

-12R 

68 L 

6r L 

yy 

y 

yy x 

y x 

(symmetric) 

,2t 1 

- 6R L 
y x 

48 L 2 
yy x 

" 6T 2 
4r L 2 

y 

4Y 3 - 

128 

69 L 

1 2 R 

6r L ‘ 

XX 

xx y 

X 

x y 


48 L 2 

6R L 

4r L 2 

(symmetric) 

xx y 

x y 
,2T . 

x y 
- 6t 2 
4Y 3 


axis 
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Appendix A 


[ - K 12- , y 


[K„] = 

22 y 


’-123 

63 L 

-12R 

6r L 

XX 

xx y 

X 

x y 

-6g L 

2p L 2 

- 6R L 

2r L 2 

xx y 

xx y 

x y 

x y 

-12R 

6R L 

- ,2T , 

-6t 

X 

x y 

2 

- 6r L 

L x y 

2r L 2 

X 

6T 2 

1 

>- 

CM 

" 123 

- 63 L 

12R 

-6r L " 

XX 

xx y 

X 

x y 


43 L 2 

- 6R L 

4r L 2 

(symmetric) 

xx y 

X 

12t i 

x y 

% 
4V 3 . 

EI../L 3 

R 

= S 3 -S3 

- r 

'J 

X 

X XX zxz 

X 

ER ./L 3 

R 

= S 3 -S3 

+ r 

ei 

y 

y yy z y z 

y 


Y] =Er/L J + GJ/10L 
Y 2 = Er/L 2 + GJ/60 


Y 3 = Er/L + GJL/30 
Y 4 = Er/L - GJL/60 


for ribs parallel fo fbe x-axis 


T. =V, + S 2 p - 2S S 3 +S 2 P -2(Sr -Sr) 
I I z zz y z yz y yy z z y y 


t 9 = Yo " (S r - S r )L 
z z z z y y 


for ribs parallel to the y-axis 


T, =Y,+S 2 P -2SS3 +S 2 P -2 (Sr - Sr) 

I I z zz X Z XZ X XX XX z z 

= y 9 - (S r -Sr)L 


X X 


z z 


The composite beam mass matrix has the form 


+[M 12 ](3) 2 
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Appendix A 


1 + tM 22 ](d} : 


where [M..] = a [K..] + a [K..] 
ij x ij x y ij y 


[M„] = M 
1 1 x 


[M 10 ] = M 
12 x 


[M_ 0 ] = M 

22 x 


[M n ] = M 

11 y 


13/35 

13e /35 
y 

A 

-1 1 L /210 

X 

0 


13L I /35 

x p 

-lie L /210 
y * 

1 1 L 3 I*/210 
x p 

(symmetric) 


L 2 /105 

X 

0 

L 4 I*/105 
x p J 

9/70 

9e /70 
y 

13L /420 

X 

0 

9e /70 

y 

9L 2 I /70 

x p 

13e L /420 

y x 

-13L 3 I*/420 
x p 

-13L /420 

X 

-I3e L /420 
y * 

- L 2 /140 

X 

0 

0 

- 

-T3L 3 I*/420 
x p 

0 

1 

0 
2; 

vL 

* Cl. 

« 

^ X 

-J 

1 

13/35 

13e /35 

y 

A 

UL /210 

X 

0 


13L I /35 
x P 

lie L /210 

y x 

-11L 3 I*/210 
x p 

(symmetric) 

13/35 

UL /210 
y 

L 2 /105 

X 

-13e /35 

X 

0 

L 4 I*/105 
x p 

0 


L 2 /105 

y 

-lie L /210 
x y 

0 

(symmetric) 


13L 2 I /35 
y p 

-IIL 3 I*/210 

y p 

A 


LI*/ 1 05 

y p 
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Appendix A 



9/70 

-13L /420 
y 

- 9e /7 0 

X 

0 


13L /420 

- L 2 /140 

-13e L /420 

0 

[M j9 1 - M 

12 y 

y 

- 9 e /70 

y 

13e L /420 

x y 

9L 2 I /70 

13L 3 I*/420 


X 

x y 

y P 

y p 


0 

0 

-13L 3 I*/420 

- L 4 I*/140 


- 


y P 

y p 


13/35 

—i —i 

i 

-13e /35 

X 

lie L /210 

0 

0 

[M 0 J = M 
22 y 

For ribs parallel to th 

(symmetric) 
le x-axis 

y 

x y 

13L 2 I /35 
y P 

11L 3 I*/210 

y p 

L 4 I*/105. 

y p 


M = pA L e = C - S e = C - S 

xx yyy zzz 


I* = (I + I )/(A L 2 ) 

p yy zz x x 

2 

I = r +1* 

P P 


and for ribs parallel to the y-axis 


M = pA L 

y y 

2 2 2 2 
r =(e +e7L 

X z y 



I* = (I + I )/(A L 2 ) 

p xx zz y y 


I = r 2 + I* 

P P 


e = C -S 
zzz 
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APPENDIX B 


STIFFNESS AND CONSISTENT MASS MATRIX FOR INCLUDING 
THE FUNDAMENTAL INTERIOR MODE FOR A RECTANGULAR PLATE BENDING ELEMENT 


The stiffness and consistent mass matrix presented here are based upon equations (26) 
and (27) as described in the main text. The basic stiffness matrix is a sixteen -degree-of- 
freedom element as referenced in the main text, and the modifying interior mode functions 
are taken as the eigenfunctions for a clamped -clamped beam. The positive coordinate direc- 
tions for this element are given in figure 9. 


At each comer of the element, the loading and displacement coordinates for the i 
comer are 


.th 


p 


' 

shear in the z direction 

M 

X 


bending moment about x-axis 

M 

y 

■ =• 

bending moment about y-axis 

M 

xyj 

i 

twisting moment 

r 

d 


displacement in the z direction 

8 


rotation about x-axis 

, x 

> = - 


8 

y 


rotation about y-axis 

9 xy . 

i 

twist 


(i=l, 2,3,4) 


From equation (26), the stiffness matrix has the form 


P. 


= -Q_ 

ab 


[K iL ! ♦ [K ] j K ci 

- k t- -j-k 

Cl I 



d. 


_ J_ 


W 

J 

o 


D = ■ 


ihr 

12(1 -V 2 ) 


where P o is the generalized force in the W Q direction. The matrix [K..] has the form 


11 

K 12 

K 13 

K 14 


K 22 

^23 

"t 

CM 

1^ 



^33 

K 34 

(symmetric) 


*44. 
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Appendix B 


where 



"k/1 6 

kb/64 

-ka/64 

8k|ab 



kb 2 /256 

-kab/256 

2k^ab 2 


(symmetric) 


ka 2 /256 

-2k^a 2 b 

. 2.2 
k2a b 


"k/1 6 

kb/64 

ka/64 

-8k|ab 


kb/64 

kb 2 /256 

kab/256 

-2k^ab 2 

[k I2 ] = 

-ka/64 

-kab/256 

-ka 2 /256 

2k^a 2 b 


8k^ab 

2k^ab 2 

2k^a 2 b 

, 2 u 2 

-k2a b 


"k/1 6 

-kb/64 

-ka/64 

-8k^ab 


kb/64 

-kb 2 /256 

-kab/256 

-2k^ab 2 

i — i 

7SI 

CO 
1 1 

II 

-ka/64 

kab/256 

ka 2 /256 

2k^a 2 b 


8k^ab 

-2k^ab 2 

-2k^a 2 b 

. 2,2 
-kgO b 


"k/1 6 

-kb/64 

ka/64 

8k|ab 


kb/64 

-kb 2 /256 

kab/256 

2k ]a b 2 

II 

1 — 1 
i^r 

i i 

-ka/64 

kab/256 

-ka 2 /256 

-2k^a 2 b 


8k^ab 

-2k^ab 2 

2k^a 2 b 

. 2,2 
k2a b 

n 

i — i 

CM 

CM 

1^ 
l l 

'k/1 6 
(symmetric) 

kb/64 

kb 2 /256 

ka/64 

kab/256 

ka 2 /256 

-8k^ab 

-2k ] ab 2 

-2k^a 2 b 

, 2,2 
< 2 ° b 
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-kb/64 

-ka/64 

-8kjab 

-kb 2 /256 

-kab/256 

-2k^ab 2 

-kab/256 

-ka 2 /256 

-2k^a 2 b 

2k^ab 2 

2k|a 2 b 

, 2 k 2 

><2 0 b 

-kb/64 

ka/64 

8k-|0b 

-kb 2 /256 

kab/256 

2k^a 2 b 

-kab/256 

k a 2 /256 

2k^ab 2 

2k^ab 2 

-2k^a 2 b 

, 2 k 2 

-k2a b 

-kb/64 

-ka/64 

-8k^ab 

kb 2 /256 . 

kab/256 

2k^ab 2 


ka 2 /256 

2k^a 2 b 



. 2 k 2 

k 0 a b 


-kb/64 

ka/64 

8k^ab 

kb 2 /256 

-kab/256 

-2k^ab 2 

kab/256 

-ka 2 /256 

-2k^a 2 b 

2k^ab 2 

-2k-ja 2 b 

, 2 k 2 

-k2<3 b 

-kb/64 

ka/64 

8k^ab 

kb 2 /256 

-kab/256 

-2k^ab 2 


ka 2 /256 

2k^a 2 b 


, 2,2 
K 20 b 
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_ _o 
C 


[(b/o) 2 + (aA) 2 +2cf 1 (C 3| -2J 2 ] 


11 


k, = k/2048 - C Q C 
C Q = 1/(1. 5881 5) 2 
C u =a 1 /(P 1 L) 


C 31 a l P l L 


11 


k 2 = k/4096 - C q C^ 


a ] = 0.98250222 


pjL =4.7300408 


The coupling stiffness matrix, (K .} , has the form 

C I 




K 


C J 

c2 

C c3 

c4 


where 


C K d! 


t K c 3 ) 


-k/4 


'-k/4 

-kb/16 


-kb/16 

ka/16 

’ » tK c2 }=' 

-ka/1 6 

-32k, ab 

k 1 


32k,ab 

1 i 

-k/4 


'-k/4 

kb/16 


kb/16 


• ; tK .}=• 


ka/16 

C*t 

-ka/1 6 

32k, ab 
. * 


-32k,ab 
* > 


From equation (27), the consistent mass matrix has the form 


P. 

i 

~F~ 

- pbab 

o 



[M ] + [M ] | M c . 

■ m t :* _ r? 

Cl I o 



d. 

__J_ 


W 

o 
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The matrix [M..] has the form 


[M.j] 


M 11 M 12 

^13 

M u‘ 

m 22 

m 23 

M 24 

symmetric 

^33 

^34 




*44. 


where 


[M,,] 


[m )2 ] 


[m, 3 ] 


— 

r 2 b 

-r 2 a 

-» 

cr 

1 


r 3 b2 

- r 3 ab 

r 5 a b 2 

(symmetric) 

2 

r 3 a 

-r 5 a 2 b 

V^ 2 . 

’ r l 

r 2 b 

r 2 a 

- r 4 ab 

r 2 b 

r 3 b2 

r 3 ab 

-r 5 ab 3 

-r 2 a 

-r 3 ab 

" r 3 a2 

vf 2 b 

_a 

i 

r 5 ah 2 

r 5 a 2 b 

' 2 u 2 
-r 6 a b J 

"'i 

-r 2 b 

-r 2 a 

-r 4 ab 

r 2 b 

-r 3 b2 

-r 3 ab 

-r 5 ab 2 

-r 2 a 

r 3 ab 

r 3 a 2 

r 5 a 2 b 

_Q 

— _1 

-r 5 ab 2 

-r 5 a 2 b 

-r 6 a 2 b 2 _ 


57 



Appendix B 


[M, 4 ] 


[M 22 ] 


[M 23 ] 




[M33] 


' r l 

-r 2 b 

r 2 a 

r 4 ob 

r 2 b 

_r 3 b2 

r 3 ab 

r 5 a b 2 

-r 2 a 

r 3 ab 

2 

” r 3 a 

-r 5 a2 b 

_r 4 ab 

-r^b 2 

2, 

r 5 a b 

r 6 a 2 b 2 

" r l 

r 2 b 

r 2 a 

-r 4 ab 


r 3 b2 

r 3 ab 

-r 5 a b 2 

(symmetric) 

2 

r 3 a 

2, 

-r 5 a b 

2, 2 
r 6 a b 

" r l 

-r 2 b 

-r 2 a 

“ r 4 ab 

r 2 b 

-r 3 b 2 

-r 3 ab 

-r^b 2 

r 2 a 

-r 3 ab 

2 

-r 3 a 

-r 5 a 2 b 

-r 4 ab 

r 5 a b 2 

2, 

r 5 a b 

r 6° 2b2 

" r l 

-r 2 b 

r 2 a 

r 4 ab 

r 2 b 

-r 3 b 2 

r 3 ab 

r 5 ab 2 

r 2 a 

-r 3 ab 

2 

r 3 a 

r 5 a 2 b 

1 

Q 
c r 

r 5 ab 2 

2, 

-r 5 a b 

-r 6 a2 b 

_r i 

“ r 2 b 

" r 2 a 

" r 4 ab 


r 3 b2 

r 3 ab 

r 5 a b 2 

(symmetric) 

2 

r 3 a 

r 5a 2 b 

r 6 a 2 b 2 
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[M34] 


&*?- 


r l 

-r 2 b 

r 2 a 

r 4 ab 

-r 2 b 

r 3 b2 

-r 3 ab 

" r 5 ab2 

' r 2 a 

r 3 ab 

2 

_r 3 a 

" r 5 ° 2b 

-r 4 ° b 

r 5 ab2 

-r 5 a 2b 

-r 6 a 2 h 

r l 

-r 2 b 

r 2 a 

r 4 ab 


r 3 b2 

~ r 3 ab 

-r 5 ab2 

(symmetric) 

2 

r 3 a 

r 5 a 2b 

2 , 2 
r,a b 


r i = c o( c o/ ,6 - 2C n) 
r 2 ° C o( C c/ 64 “ C 1 1 (C 21 + C 1 \ /4] ) 
r 3 = C o ( C </' 28 " C 1 l C 2 lV 2 
r 4 = c 0 ( c o/ 256 - c 2i - c n/’ 6 ) 
r 5 = C c(V 256 - C 21- C 21 C n/ 4 V 4 
f 6 = C o( C </ 512 - C 2 l )/ 8 

C 21 =C l / C 31 


The coupling matrix, [M has the form 

0 I 


r m 


cl 


M « 
c2 -| 


C I 


M 


c3 


M 


c 4 J 


where 
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CM ,} = 


CM c3 } = 


h 


~ r 7 

_D 

00 

L. 

{m c2 }- 


r 8 a 


-r 8° 

_r 9 ab 


r ? ab 


-r 7 


_r 7 

r 8 b 

Cm c4 }=' 

r 8 b 

r 8 a 

-r 8 a 

. r 9 ab 


_Q 

w 

1 


r 7 = 2C o (C o/ ,6 - c n> 

r 8 = ~ C ll C 2p 

r 9 = W 56 ' C 21> 
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APPENDIX C 


COMPUTER PROGRAM DESCRIPTIONS, FLOW CHARTS, 
AND CARD IMAGE LISTINGS 


This appendix contains the computer program descriptions, flow charts, and listings for 
the one-dimensional panel and the two-dimensional panel modal analyses. The program 
descriptions list the purpose of the program or subprogram, the subprograms required, 
definition of primary variables, restrictions, accuracy (when applicable), and compiled 
size (octal) of the program for the NASA Langley CDC 6600 computer. 
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ONE-DIMENSIONAL PANEL ANALYSIS: MAIN PROGRAM AND SUBPROGRAMS 
IN ORDER OF APPEARANCE W 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

INPUT DATA: 
VARIABLES: 


RESTRICTIONS: 
SIZE: 005152 g 


Appendix C 


PROGRAM BMPROP (MAIN) 

(See Figure 1 for Flow Chart) 

Main program for computing natural frequencies, normal mode 
shapes, modal shear distribution, and modal bending moment dis- 
tribution along the centerline of a one-dimensional panel array 
undergoing cylindrical bending (beam analogy). A fundamental 
clamped -clamped mode across the panel width is assumed. Elastic 
supports are introduced as lumped spring -mass constants and lumped 
masses are introduced as a support with zero stiffness. Clamped or 
elastic supports can be introduced at either end of the structure. 

ELEM, ASSY, LOC, NONDIM, ORDER, DELETE, FREQ, NROOT, 
EIGEN, SSBM 

See Appendix D for description - 

NCASE, NDATA, NBAY, NSUP, IOUT, IBL, IBR, NINT, BW, 

PR, NEL(I), El(l), WB(I), BL(I), NCP(I), SL(I), SC(I), SR(I), 

RL(I), RC(I), PR(I) 

NDEL(N) - Constraint vector for an N degree-of-freedom system. 

If NDEL(J) = 0, the J* coordinate is unconstrained; if NDEL(J) = 1 
the J^ 1 coordinate is constrained (deleted from the equations of 
motion). 

R(N) - Consistent mass matrix - dimension N(N + l)/2 for N degree 
of-freedom system. 

S(N) - Stiffness matrix - dimension N(N + l)/2 for N degree-of- 
freedom system. 

SI (4), S2(4), S3(4) - Dummy arrays for assembling the element 
stiffness in matrix S(N) and the element mass in matrix R(N). 

DUM3(2293) - Dummy array for eigenvalue calculation. 

For the declared size of the arrays 
2(SNEL(I)+ 1) < 38, NBAY < 5. 
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: - Pr og ram bmpkqp ( iuPi it. output #tapes= input #tApe6=qutput). 

□ lMENSIOlxl BL(5) #N£L(5) ,NDeL( 38) *R(741) fS(74l) 

dimension sl(6) #sc(g) ,sru) »RL(6) #RC(f.) #rr(g) 

dimension si(4) tS2(4) »S3(4) »ncp(6> 

; D IMENS I ON £1(5) #Wa(G) .RLMT (5) #QUM 3(22 q3) 

.common r#s 

common- e i > wb # blmt , bvj • pr - ' : ^ • 

CQMMONt TL»SL#SC#SR#RL»RC»RRfNEL#NCP 

COMMON |BL» Si »S2»S:^rNDEL>DUM3 

IN=5 ; 

1 0 — 6 } ; 

. ICASE=1. 

REAO( I Nf 120) NCASE 

100 HEAD < IN $120) NDATa#NBaY»NSUP» I OUT# lBL f IBR#N iNT 

READ ( IN » 125) BW#PR ; 

READ(lNfl30) (NEL(I) ,EI(I) »WB(I) »BL(I) ,I=i»NBAY) 

■ IF(NSUP) Il0#110#l05 ; 

. . 105 DO 110 1=1# NSUP 

READ (IN >130) NCP(I) ^La) >SC(I) #SR(I) 

read ( in » i<+o ) rh d #rc( d #rR(I) 

no continue . . 

W« I TE ( 10 » 145 ) 

WRITE (I0» 150) ndata 
Vvr I TEE ( 10# 155) NBAy# NSUP 
WRITE (io #160) BW # PR 
WRITE (10# 165) 

WRlTE(I0#l7Q) ■ 

DO 175 I=1#NBAY 

WrIT£( I0#18Q) I#NEL(I) > EI(I)#WB(I)#BL(I) 

1 175 Continue 

Do 161 1 = 1 # NSUP 

write ( io # i82 ) i#ncp(d 

Wr ITE ( 10# 183) SL(I),SC(I)#SR(I) 

~ write (10# 184) RL ( I ) # RC ( I ) » RR ( I ) "1 

181 continue : - ; 

tl=o.o j 

Tk~ 0 » 0 

Tm=0.0 “ 

NcT=0 

Do 200 1 = 1 » NBAY • 

BlMT ( I ) =3L ( I ) /FLOAT ( NEL ( I )) 

! TL=TL+6L(l) 

1 NCT=NCT+NEL(I) 

200 continue 

NgP=NCT+1 
: NC0=2*N6P 

NuP=NC0*(NC0+l)/2- ; 

! Do 240 I0P=1 » 2 

Do 210 1 = 1 » NUP 
R ( 1 ) =0 . 0 

210 Continue ! 

NuT=0 : : ; 

Do 220 J=1»NBAY } 

CaLL ELEM(U» l0P#Si»S2#S3) 

■ Nu=NEL(J) t 

NUT=NUT + NU — ~~ • 

Do 220 1=1 #NU 
ICN=2* ( I+NUT-NU) -1 

PROGRAM BMPROP: CARD IMAGE L i STI NG 1/3 I 


IcS=ICN+2 " ' 

call assy (Mco# icn, ics.r»si*s2»S3»2> 
220 Continue' ~ T 

IF ( IOP-1 ) 250 >230 #240 ' 

230 Oo 235 I=1#NUP 
S(I)=R(I) 

. 235 Continue 

240 Contin ue 

Ip(NSUP) 250 * 250 t 245 
245 DC 250 I=1»NSUP 

o=ncp( i ) ; 

K=J+1 

jj=J+( J*J-J)/2 
JKrJ+(K*K-K)/2 
KK^K+(K*K-i<)/2 
S(UJ)=S( JJ)+SL( I) 

S( JK)=5< JK")+SC (I ) : ! 

; S(KK)=S<KK)+SR(H 

R ( J J ) =R(JJ)+RL(I)/306.0 
R ( JR ) =R ( JK ) +RC ( I ) /3B6 . 0 
R(KK)=R(KK)+RR( D/386.0 
250 Continue 

C NOnDTmEnSIONALIzE st IFFNESS MATRIX 
Oo 255 I=2»NC0#2 

! il=I+( 1*1-17/2 “ 

TK=TK+S(II) 

: rt^=rw+R(iTT 

255 continue 

i rr-'tL7FL-6Airwcn 

Tk=TK/FLOAT(NCT) 

” |> i= TM/FLOAT{NCT) ; 

Call nondiM(s»nco,tl»tk) 

CaLl nonuim<r» NCU, IL»TMJ 
;c APpLy CONSTRAINTS 
• i\iL)L=0 

Oo 261 J=l»NCO 
No ELI J)=CJ 
261 Continue 

. m I UL ) 26472647 262 - ~~ 

i 262 N[)EL(1)=1 

: NoEL“(2)'-I 

i Do 263 I=1'NSUP 

NcFTTl^NOP ( I) -2 

263 Continue 

264 If(IUR) 270*270*265 

265 IC=NCO-i 
~ NdLLUC) = 1 

NdEL(NCO)=1 

Call order (s» ndel, ngp,ndlj 
Call ORDER (R»NDEL»NGPf NDL) 

270 Nco=NCO-NOL 

NuP=NC 0*(NC0+1)/2 
1FUUUT-1) 3 15f275^275 
i 275 WRITE ( 10 r 280 ) 
i KTTTTE1T01 265) NCO , T K » TL 

WRlTE(I0r290) (S(I),I=1»NUP) 

f wrtte( iOf'300) ' 

WrITE(IO»305) NCO,TM,TL 
“ WrI |H 10 » 290 ) TRTTTTTTTTTTuPl 

iPROGRAM BMPROP; CARD IMAGE L 1ST I NG 2 /3 i 
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c compute eigenvalues and eigenvectors 

, 315 WRITER 10*320) TL>TK,TM . ' ■ ■ v ■ , 

i call freo (S#,r » tl * tk r tm»ndaTa » nco ) 

IEU0U1-1) r00»10Q»26n ; •' • 

260 call ssom(nco»nbay»nsup» ibl,nint»ndata) 

, I CAS E = IC ASE-kl 

. IF<NCASE-ICA§E) 205#100»100 
; 205 continue . H 

120 FoKMAK 15,61^) 

i?5 Format ( 3X>2Ej2. 5? 

. 130 FoRMAT( I3f3£i€.5) 

140 Format ( 3x,3Ei2. 5) 

l 145 FORMAT (1H1»7X»47HFREE VIBRATION OF A one DIMENSIONAL panel ARRAY) 
?' 150 FcRMAT(/#2bX»9H0ArA CaSE»I4) 

1rj5^Fo«r^T(7»4X#i5l : <NUMbER OF bAy'S=# I3»19X, 19HNUMBER Op SUPPORTS:, 13) 

; 160 FORMAT (/»4X» 12HPANEL ^IUTh=#E 12.5»7X» ibHPOISSON'S RATIO: » El2 . 5 ) 
16F'FoUMAiT/»5X73HBAY,5x»9HNUMljER OR » 3X , 7HBENQ IflG * 6X » jOHWElGHT PER* 
18X * 3HBAY ) 

170 format ( 4x , ohnumber * 3X , hhelEments » 4 X , qhRIGidity,5X, 

19HUNIT AREA#aX»faH L ENGTH#/) 

180 F0RMAT(/»5X»l3»7X,l3»5X*Ei2.5#lX»El2.5*3X»Ei2.5) 

> 182 FqRMAT (/» 4XrllHSUpP0RT NO. » 13»25X, 17HINPUT COQRDInATE=»I3) 

183 FORMAT ( 4X , 4HKZZ="i £l2 • 5 > 2X , 8HKZTHETA: , E 12 . 5 » 2X r 7 HKtHETA= , Ei2 . 5 ) 

104 FoRMAT(4Xf4HIZZ=.El2.5»2X,8HlZTHETA=,El2.5*2X»7HlTHETA=»Ei2.5) 

280 Format ( lHlf4X»16HsTlFFNESS MATRIX) 

265 FqRMAT (/f5X»4HNC0=» 14 »2X>3HTK=fEl2»5 > 2X>3HTL=»E12.5»/) 

290 Format ( 8 E 12 . 5 ) 

300 format (ihi» 4 x»hhwass matrix) 

305~FoRMAT (/r5X*4HNC0-» I4,2X»3HTN=»E12.5»2Xr3HTL=»E12.5»/) ' ' T 

310 Format ( 6 E 11 . 4 ) 

; 320 FORMAT (3E 12 . 5 ) 


End 
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PURPOSE: 


SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
SIZE: 0001 52 g 
REFERENCE: 


SUBROUTINE ELEM (I, IOP,Sl , S2,S3) 

To compute the stiffness (IOP = 1) and the consistent mass matrix 

th 

(IOP = 2) for a panel element in the I bay of the structure. The 
partitioned stiffness or mass matrix is assigned to the arrays S1,S2, 
S3 where SI is the direct stiffness or mass term at end 1 of the ele- 
ment, S2 is the coupling term between ends 1 and 2, and S3 is the 
direct term at end 2 of the element. 


None 

I - Panel Bay Number 

IOP - Logic number: IOP = 1, stiffness matrix for element is cal- 
culated; IOP = 2, consistent mass matrix is calculated. 

SI, S2 , S3 - See 'Purpose' above , 

El(l) - Bending rigidity of element 

WB(I) - Weight percent length of element 

BL(I) - Length of element 

None 


Reference 1, pp. 9- 11. 
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FLOWCHART: SUBPROGRAM ELEM (I, IOP, SI, S2, S3) 
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' ; * ' « ' I 

Subroutine elem d kop,si,S2»S3) 

iO I MENS I ON Sl(l) »S 2 < 1 ) »S3(1>»DUM1 (1482) 

pi MENS I ON El (5) »WB<5) »OL(5i 

PoMMON DUMl»El»W 8 -, BL»BW»PR 

A=bL(I) 

ftS-A*A • ; 

|gVdOP-l) 10, 10 » 15 

: 105|=El(I)/(A*AS) - 

f jfi ( 1)=12.0*0 ~ ~~ 

IjSl ( 2 ) -6 « 0*D*A . 

Si ( 3) =4 » 0*D*AS 

5 2 d ) =- 12 « 0*0 

S 2 ( 2 )=- 6 . 0 *D*A 

S 2 t 3 ) = 6 • 0*0*A 
S 2 l 4 ) =2 . 0*D*AS 

5 3 d ) = 12 • 0*0 
SyT 2 )=- 6 . 0 *D*A 
S 3 i 3) =4 • 0*0* AS 
RETURN 

15 0=WB( I ) * A/ 386 • 0 
Si ll)=13. 0*0/35.0 
Si(2)=11.0*0*A/210»0 
S 1 T 3 T = d*as/ioS.o 

S2(l)=9. 0 * 0 / 70.0 
S2<2)=l3.0*0*A/420»0 
S 2 (3)=-13.0*0*A/420.0 
S2t4)=-U*A/140,0 
S3d)=13. 0*0/35.0 

S 3 "r 2 )=- 11 . 0 *D*A/ 210.0 ' ■ 

S 3 ( 3) =D*AS/105 . 0 

RE'TDRR i 

end • 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 0001 43 g 


SUBROUTINE ASSY (N, ICN, ICS, A, El , E2, E3, NLM) 

To assemble (add) the partitioned element stiffness or mass matrices 
(E1,E2,E3) in the appropriate location of the global stiffness or 
mass matrix A . 


LOC 

N - Size of matrix A (See Subroutine LOC). 

ICN - Initial coordinate for diagonal location of El in A. 

ICS - Initial coordinate for column location of E2 in A for row ICN. 

NLM - The number of coordinates at each grid point for the ele- 
ment described by El, E2, and E3. For the one-dimensional 
analysis NLM = 2. 


None 

Not Applicable 
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STEP 4.1 


STEP 4.2 


STEP 4.3 


STEP 4.4 


STEP 4.5 


STEP 4.6 


FLOWCHART: SUBPROGRAM ASSY (N, ICN, ICS, A, El, E2, E3, NLM) 
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FL0W CHART: SUBPR0GRAM ASSY (N, ICN, ICS, A, El, E2, E3, NLM) 
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! Appendix C 


; OT7<UUTTFir~ ASi>Y(N, ICN, ICS, A,E1 e L2r L3#NLM) 
Dimension A(i) » eicd »E 2(i) »E3(i) 
l^O 20 K = 1»NLM 7 1 

Dp 2U L = 1 , NLM . ; 

call loc(k#l,kls,nlm,d 

call loC(k*l>klg,nlm,q) 

InK=ICN+K-1 

InL=ICN+L-1 

Call loc<ink»inl»nkl»n»i> ' 

IsK=ICS+K-l ) 

ISL=ICS+L-1 

call loc(isk,isl, ikl»n»i>. 

CaLl L0C(1NK#ISL»IKS»N»1) 

IF(K-L) 15*10*10 , 

10 A(NKL)=A(NKU+EKkLS) : 

A(IKL)=A( IKL)+E3(KLS) 

15 A(1KS)=A(IKS)+E2(KLG) j 

20 continue | 

return t ‘ 1 

End i 


■SUBPROGRAM ASSY: CARD IMAGE LISTING _ j 
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SUBROUTINE LOC (I, J, IR, N,MS) 

PURPOSE: To calculate a single subscript, IR, for a double subscripted square 

array. A, of size NXN for a row-column location (I, J). MS is a 
logic number for calculating the subscript IR. 

MS = 0, General Storage, array is assumed to be dimensioned as 
A(R) R = NXN 

MS = 1, Symmetric Storage, array is assumed to be dimensioned as 
A(R) R = N (N + 1 )/2 

MS = 2, Diagonal Storage, array is assumed to be dimensioned as 
A(R) R = N 


SUBPROGRAMS 

REQUIRED: None 

RESTRICTIONS: None 

ACCURACY: Not Applicable 

SIZE: 000057g 

REFERENCES: Reference 3 
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FL0W CHART: SUBPROGRAM L0C(l, J, IR, N, MS) 


Appendix C 
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Appendix C - 


. c 


subroutine loC(i»j,ir,n»ms> ; 

c 

c 


This subroutine compresses a two dimensional 
array into a one dimensional array 

c 

c 


Ms=o is for general storage 
M£=i is for symmetric storage 

c 

: C 


(|s=2 IS FOR DIAGONAL STORAGE . 

Ix = l 
1 Jx=J 

i 

1 

10 

IFIMS-D 10 r 20 * 30 
IrX=N*( jx-d+ix 

1 • 

20 

Go TO 36 

IF(IX-JX) 22 » 24 » 24 

! 

i 

22 

IrX=IX+(JX*JX-JX)/2 
GO TO 36 


. 24 

IRX=JX+( lX*IX“IX)/2 
Go TO 36 

i 

.30 

If<X=0 

IF(IX-JX) 36 » 32 » 36 


32 

36 

IRX=IX 

IR=IRX 

RETURN 
! End 


[SUBPROGRAM LOC: CARD IMAGE LI 1ST 1 NG 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 0001 33 g 


SUBROUTINE NONDIM (A,N,TL,TIC) 

To nondimensionalize the free-free stiffness and mass matrix of the 
structure. 


LOC 

A(N) - Stiffness on mass matrix to be nondimensionalized 
N - Size of array A (MS = 1 in LOC) 

TL - Arbitrary length parameter (taken as the average element 
length of the structural idealization) 

TK - Arbitrary stiffness parameter (taken as the average rotational 
stiffness of the structural idealization) 


None 

Not Applicable 
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FL0W CHART: SUBPR0GRAM N0NDIM (A, N, TL, TK) 
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fr • ■ • . ' ' 

| Appendix C 


:-l : : : : •- 

' *, • SUBROUTINE MONO IM(A»N»TL»TK) I 

1 ) DIMENSION A ( 1 ) 

1 L>0 3 I=2»M*2 

* Ii = I-l _ . j 

" 12=1+1 

• 3 13 = 1+2 , 

\ IU=U+(Il*Il-Il)/2 

■\ ni=ii+(i*i-i)/2 1 

A( IIl)=TL*TL*A(IIi) 

A ( 1 1 1 ) = TL*A (III) 1 

• — Tf(I-N) 2 » 3 » 3 

j - 2 Il2=U-t-(l2*l2-I2)/2 

Il3=Il+(l3*l3-I3)/2 

1X2=1+ ( 12*12-12) /2 

A( H2)=TL + TL*AT1 12) ; 

A( Il3)=TL*A( 113) 
a7TI2 ) = TL*A(H2) 

3 continue 

do 4 1=1 #N 
Do 4 J= I • N 

CaLl L0C(I»J»'IJ»N,1) '■ 

A(IJ)=A(IJ)/TK 

txoNTTNue: 

return 

r e n° : ' 


SUBPRO'GRAM'NO KD I M:"~CARb IMAGE LISTING 



Appendix C 


PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 

ACCURACY: 
SIZE: 000244g 


SUBROUTINE ORDER (A, NDEL, NGP, NDL) 

To remove (set to zero) specified (constrained) coordinates in array 
A, reorder the array A, and calculate the new size of array A. 


DELETE, LOC 

A(l) - Stiffness or mass matrix of structural idealization 
(symmetric storage mode: LOC) 

NDEL(I) - Array of logic numbers: See BMPROP 

NGP - Number of grid points of the structure 

NDL - Number of coordinates removed by this subprogram 

NC(38) - Array of coordinate numbers for which NDEL(I) = 1 

NGP is assumed to be equal to or less than 19. 

NC(I) must be dimensioned the same as NDEL(I) in program 
BMPROP 

Not Applicable 


80 



Appendix C 


^ BEGIN J 

* _ 

INITIALIZE DATA 


SELECT GRID 
P0INT, 1C 


I 


COMPUTE SUBSCRIPTS FOR 
COORDINATES AT GRID POINT 1C 
TRANSLATION SUBSCRIPT, J 
ROTATION SUBSCRIPT, II 


C0NTRAINT 
NDEL(I) > 0 


NDL = NDL + 1 CUMULATIVE C0UNT 

0N DELETED C00RDINATES 

NC(NDL) = 1 C00RDINATE N0. DELETED 


STEP 7.1 
STEP 7.2 

STEP 7.3 



STEP 7.4 


STEP 7.5 



STEP 7.6 


STEP 7.7 


STEP 7.8 


G0 T0 STEP 7.2 


BEGIN REM0VAL 
0F C0NSTRAINED 
C00RDINATES 


STEP 7.9 


FL0W CHART: SUBPR0GRAM 0RDER (A, NDEL, NGP, NDL) 
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FLOWCHART: SUBPROGRAM 0RDER (A, NDEL, NGP, NDL) 
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Appendix C 


SUBROUTINE DELETE (A,N,J) 

PURPOSE: To delete (set to zero) all elements in the J^ 1 row column of array A. 

SUBPROGRAMS 

REQUIRED: LOC 

VARIABLES: A(N) - Single subscripted array (symmetric storage mode: LOC) 

N - Size of array A 

J - Row-column number of elements of A to be set to zero 
RESTRICTIONS: None 

ACCURACY: Not Applicable 

SIZE: 000034g 
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Appendix C j ‘ SUBROUTINE ORDER ( A » NDEL » N6P » NOL ) 

— V ; D I MENS ION A ( 1 ) r NDEL ( 1 ) , NC ( 3a ) 1 

~ NoLzlT ~ i 

NcO=2*NGP 
DO 1 1=1 »NCO 
i I NO ( I ) =0 

■ i continue 

i r Do 5 ic=i » ngp i 

■ f. I=2*IC-1 - 

■ j j Il=2*IC ' 

i IF (NDEL ( I ) ) 3»3»2 

2;NpL=NDL+l . 

iNC ( NOL ) = 1 

! 3 |lF(NDEL(Ii) ) 5 i5»4 - . . : : 

i 4 ^UL=NDL+1 ' ' ~ 

Np ( NOL ) = 1 1 ; 

: 5 continue 

! Do 11 K = 1»NDL ■' 

M=NC(K) 

ML=NCO-K ■ 

N=ML+1 ; 

; call uelete(a>n>M) j 

Do 9 I=M#ML 
I 1=1+1 
DO 9 J=1 • N 
Call loC(i» j,ij»n,u 
IF(I-J) 7 r 6 # 7 

6 Ul=J+l 

TTaEl^LOC ( 1 1 # Jit lJi»N»i j ^ ! 

60 TO 8 I 

7 CaLL LOCTIl,J»IJl,N,l) : 

a a ( I j ) = a ( i ji ) : 

■"9 -Continue 

DO 11 1=1 » ndl 

. iF(Trc(i)-M) rrTTTTTo : ; 

10 Ncd)=NC(I)-i 

11 Continue ' 1 ;■ 

M=NC0-NDL+1 

DO 12 K = M f NCO ; | . 

call delete (a*nco,K) 1 

12 ConTTNUE" 

RETURN 

•• EmD • 1 


subroutine delete (A,n,j> 
dimension Ad) 

do i k=i >N 

C all lo C(k>j>kj»n»d 

A ( K J ) = 0 • 0 

i continue 

return 

end 


SUBPROGRAMS ORDER AND D E LETE: C AR D I MAG E LIST I NG ! 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 000256 g 


SUBROUTINE FREQ (SS / RR # TL / TK / TM / NCASE, N) 

To convert the stiffness, SS, and mass, RR, matrices (symmetric 
storage mode: LOC) to the stiffness, S, and mass, R, matrices 
(general storage mode: LOC); call for eigenvalue calculations, 
and write the eigenvalues (frequency) and eigenvectors. 

LOC, NROOT (NROOT requires EIGEN) 

SS(NUP) - Dimensionaless stiffness matrix 
RR(NUP) - Dimensionless mass matrix 

TL - Length parameter used in nondimensionalization 

TK - Stiffness parameter used in nondimensionalization 

TM - Mass parameter used in nondimensionalization 
NCASE - Data case identification number 
N - Number of coordinates 

For the declared size of the arrays, N £ 19 

See NROOT and EIGEN 
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FLOWCHART: SUBPROGRAM FREQ (SS, 


RR, TL, TK, TM, NCASE, N) 


36 




! Appendix C 


subroutine freq(Ss*rr » tl#tk»tm»ncase»n) 

! Dll MENS I ON SS(1) *RR( \ ) ,EVL(3b) »£VC( 115ft) 

! dimension stuse) f R(ii56 ) #dumk342) 

; dimension idumkii ) ,dum2(i7) : 

common evl*evc»dumi,idumi,dum2»s*r 

logo . 

NuP=N* (N+l )/2 

WrITE ( IO» l2Q ) , 

WrITE( 10 » 125) NCAsE 

. 200 DO £05 1 = 1 >N 

Do 205 J=I»N 

CaL^. LpC ( I * J > IJ»N»0) ' 

CaML LOC ( I » J * I JS » N » 1 ) 

S(IJ)=SS(IJS) 

K(IJ)=RR(US) 

205 continue 

call nroot <n»s*r,evl,evC) 

Dp 220 I = 1*N t 

IF(EVL(I) ) 210*215*215 

210 EyL(I)=ADS(EVL(I) ) ' 

215 EvL(i)=0.l59i55*SQRT(TK*E\/L(I)/TM) 

WRITE ( 10*225 ) EVL(I) ' 

WRITE ( 10,230) 

IX=N*(I-1)+1 

. Iy=N*I 

write ( io * 235 ) (Evcu) , j=ix*iy) 

220 continue 

1 RETURN - __ 

: 120 FoRMAT(1H1*27X*24HFREE VIBRATION OF A ONE-* 

123HDIMENS10NAL PANEL ARRAY) 

125 FORMAT (/»44X»9HDAT A CASE»I4#/) 

225 FqRMAT,(/*5X>10HFReQUENCY=: > Ei2»5*1X*3HhZ« ) 
230 Format </»5X»iohmode shape*/) 

: 235 Format (5X»6El2. 5) . 

236 FoRMAT(5El5.8) 

end 



■SUBPROGRAM FREQ: CARD IMAGE LISTING , 


I ■ , ‘ ' i 
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Appendix C 


PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 0003 12 g 
REFERENCES: 


SUBROUTINE NROOT (M,A,B,XL,X) 

To compute the eigenvalues and eigenvectors of a real symmetric 
matrix of the form B-inverse times A. 


EIGEN 

M - Order of square matrices A, B, and X 
A - Input matrix, MXM, (stiffness matrix) 

B - Input matrix, MXM, (mass matrix) 

XL - Output vector of length M containing eigenvalues of B-inverse 
times A 

X - Output matrix, MXM, containing eigenvectors column wise 
See EIGEN 
See EIGEN 


References 3 and 4 
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SUBROUTINE NROOT (M,A»B»XL»X) 

Dimension A(D»b(i) »xl( 1) »X(i) 

i ^ppenaix v* 

c 

CoD R 0TE~E TGL’NV A L U e 5 ' and"H1GeN VECTORS" oF 'B 

K=1 



L=M*(J-i) 


00 ioo r=i 7 j 

1 1.4 *•! 



K=K+1 % 

■ 

100 B ( K ) =13 (L } 1 

C“ 

j 

The matrix b is a real symmetric matrix 7 n 

My = 0 i ' 

i c 

" " CaL L EI6EN(B»X-»M»MV) ' “■ 

form reciprocals of square root of eigenvalues. The results! 

~ Z ~ 

" arl rklmmluplilu By the assoCIAieo eigenvectors. 

LrO 


W~TrO~J=T7W 

L~L+U 


“TlO - XLCj)=1.0/bQR"f (At3S(0(L) ) ) 

K = 0 

uo 115 j=i»m ' ' 

Do 115 I= 1 »M 

K _- K + 1 - ; 

115 B(K)=X(K)*XL(J) 

Z FoRFF (B**(-i72) )PrIME*A*(B** (“1/2TJ 

DO 120 I = 1 r M 
|N)? - 0 

Do 120 J=1»M 

; rTr=M^rr^rj 

L-M* ( J-l ) +1 
X(L)=0>0 * 

Do 12(1 K=1*M 

. =hi + i • ! ' 

N2=N2+1 

“”120 X(L)=X(U+b(Nl)*A(N2) 

L-0 

Uo 150 J=1#M 
DC 130 I=1»J 

|> j 1 = j L-M — — 

N2=M*(J-1) 

_ r ^ 1 _ I 

A ( L ) -0 . 0 

~~ ^W TJO^K-l »~M 1 

Ni=N1+M 
R2=N2+1 

130 A(L)=A(L)+X(Nl)*B(N2) 

“C CppOTF: EIGENVALUES AND EiGeNVECIOKS OF MaTr-IX A. 

Call eigen (a»x*m,mv) 

L^O 


DO 140 1 = 1 * M 

" l-l+ i : " 

140 XL(I)=A(L) 

~C . CoMPUlE NQKMAUZtD EIGENVECTORS : 

; Do lbO I = 1 r M 

DO 150 0=1 * M 

/ — Ri-i-iM ; 

L ' ~~ [SUBPROGRAM NRO OTi CARD IMAGE LISTING 1/2 1 


i 





Appendix C 

I 

) 

> ' J-D + I 

; a ( d -0, 6 . 

• ( Do i50 Ki=i»M 

' Ni=Nl+M , 

N?„=N2 + 1 \ 

I 150 A(L)=A(LjUB(Nl)*X(N2) 

LrO \ 

K=0 \ 

DO 180 U=\ »M 

Sum v= o.o i 

'* Do 170 . 1 = 1 i M 

•i LzL-t-1 

■ r , 170 SuMV=SUMV+A(L)*A(L) 

■ 175 SuMV=SQRT (SUMV) 

Do 180 1 = 1 * M 

v K-K+l 

180 X(K)=A(K)/SUMV 

RETURN 


end 










Appendix C 


PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 

ACCURACY: 

SIZE: 000522_ 

O 

REFERENCES: 


SUBROUTINE EIGEN (A,R,N,MV) 

To compute eigenvalues and eigenvectors of a real symmetric 
matrix by the Jacobi method. 


None 

A - Original matrix (symmetric storage mode: LOC), destroyed 
in computation. Resultant eigenvalues are developed in 
diagonal of matrix A in descending order. 

R - Resultant matrix of eigenvectors (stored column wise, in 
same sequence as eigenvectors). 

N - Order of matrices A and R 

MV - Input option 

0 Compute eigenvalues and eigenvectors 

1 Compute eigenvalues only (R need not be dimensioned but 
must still appear in calling sequence) 

The original matrix A must be real symmetric (storage mode I: 
LOC). Matrix A cannot be in the same location as matrix R. 

At each step of the diagonalization a norm is calculated and the 
diagonal ization continued with the magnitude of the off-diagonal 
term is sufficiently small to insure convergence. 


References 3 and 4 
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Append ix C "SuBrODTTWe EIGEiMTa # R » N * M \T) 

— - dimension .A(D »R(i) 

~c ‘ generate iuenity matrix ' 

5 KaNGE=1.0E-6 

IfIMV-1) 10 » 25 » 10 

10 I G--N : _ i 

DO, 20 J=1 »N \ 

iQ^IO-t-N - 

Do ; 20 1 = 1, N 
Ij-IQ+I 

. , ~KTT^')=0.0 

J iFtk-P) 20*15*20 . ■ 

' 15 R ( 1 0 ) = 1 . U 

20 Continue 

~C Co^hUTE .'INITIAL AND FINAL NORMS (ANORM AND ANORMX) 

25 AnO||m= 0.0 
Do 35 I=l7R 
Do 35 J=I ,N 
IfTI-J) 30 05 * 30 
30 Ia=I+(J*J-J) /2 

ANORM=ANO'i<"M + A TTA ) *A ( 1 a ) 

. 35 continue 

If (ANORM) 165*165,40 
40 AN0RM=1.414*SQRT (ANORM) 

- AnRMX=ANOrM*RANOL/RLOAT (N j 

! c initialize indicators and compute threshold, thr 

• InD=o 

j ThR=AN0RM 

• 45" lhR= I HR/FLOAT (N) 1 

50 L=1 

55 M=L+i 

c compute sin and Cos 

- 60 ->i u =lM*|V|-M)/2 

! Lq=(L*L-L)/2 

qq-C+MO 

62 IF ( ABS ( A ( LM ) ) -THR ) 130*65,65 

65 InO~T 

" LL=L+LG 

_ FigrT'i + fviQ 

X=0.5*(A<LL)-A(MM) ) 

66 Y=-A (LMT/SORH A (LM) *A (LM)+X*X J 
IF ( X ) 70*75,75 

70 Y=-Y 

75 SlNx=Y/SQRT(2.0*(l.0+(SQRT(l.0-Y*Y) ) ) ) 
bjN X 2=SiNX*SlNX 
78 CoSX=SQRT(1.0“SINx2) 

C0SX2=C0SX*C0SX 
SlNCS=SINX*COSX 
~C ROTATE l and m COLUMNS 

IlQ=N* (L“ l) 

I(v]Q=N* (M— i ) 

■ DO 125 1 = 1 * N 

“ Iq=U*I-1)/2 

IF(I-L) 60*115*80 
80 IfTTHM) 85*115*90 
85 Tm=1+MG 

, ; oo ro 95 ; 

_____ SUBPROGRAM EIGEN; CARD IMAGE LISTING 1/2 


90 IiM— M+1Q j Appendix C 

' 95 IjF t I -L ) 1Q0,105,1Q5 

100 1L=I+LQ 

6b TO 110 • 

105 IL=L+IG 

11 0 Xr A ( XL ) *COSx~A ( 1M ) I K'X / 

AtlM)=A(lL)*SlNX+A(IM)*COsX 

: A(.IL)=X ! ; ; 

115 IF (MV-I) 120*125*120 ' ] 

IPO II R=ILQ+I ; j 

1MR=IMQ+I 

X=R( ILR )*CQSX-R( IMR)-*SINX 

R(IHR)=R(ILR)*S1Nx+R(IMR)*C0SX I 

R( lllR )=X ; . 

125 continue 

X=2|o*A(LM)*SlNCS 

Y=/^LL)*C0SX2+A(Mm) *SINX2-X 

' X =)&A(LL)»SlNX2 + A(MM)»C0Sx2 

A(Lm7=(A(LL)-A(MM) )*SINCS+A(LM)*(C0SX2-SInX2) 

A(LL)=Y " ; 

; A(MM)=X ' , 

:_c te sts for completion _ 

; C TEST FOR M = LAST COLUMN 

• i 130 IF(M-N) 135 >140 *135 , 

135 M=M+1 ■ ! 

! Go TO 60 

c . test for l = second from last column 1 

: 140 IF(L-(N-1)) 145*150*145 j 

145 LzL 41 i 

60 TO 55 ! 

150 IF(INO-I) 160*155,160 

155 1nD=U ; ; 

Go TO 50 

c Compare threshold with final norm 

160 iF(THR-ANRMx) 165,165,45 

_C S ORT EIGENVALUES ANp EIGENVECTORS ' 

165 IO=-N 

DO 185 1 = 1 * N 

lQ = 10+N 

Ll=I-K I*I-I)/2 

Jq=N*(I-2) 

DO 185 J= 1 * N 

Jq=JO+N • 

: MiM- J4 ( J* J- J ) /2 

IF(A(LL)-A(MM)) 170,185*185 

170 XrA(LL) ' ' 

A ( LL ) - A ( MM ) ~ ----- 

A (MM) rX ; : 

IF(MV-I) 175*185*175 

175 Dp 180 K-l , N • . r ’ ] 

ILR=1Q+K , 

imr=jo+k 

XrR(ILR) 

R( ILR)=R( IMR) 

" 180 R ( IMR)=X ; 7 : 

• 185 Continue . • 

return 

end ■ ; ' 

" "SUBPROGRA M EIGEN : CARD IMAG E LISTING 2/2 i , 


Appendix C 


SUBROUTINE SSBM (NCO, NBAY, NSUP, IBL, NINT, NDATA) 

PURPOSE: To compute the displacement, slope, shear, and bending moment 

distribution along the centerline of the panel row for the fifteen 
lower frequency modes of the structure. The displacement, slope, 
shear, and bending moment distribution are calculated at each 
element node point and at NINT points interior to the element. 

All values are normalized to the maximum value and printed. 

SUBPROGRAMS 

REQUIRED: ELEM 

VARIABLES: NCO - Number of coordinates for the structural idealization 

NBAY - Number of panel bays of the structure 
NSUP - Number of elastic supports 

IBL - Logic number: IBL = 0, the left end of the structure is 

free or elastically supported; IBL = 1, the left end of the 
structure is clamped. 

NINT - Number of points interior to an element for which inter- 
polated values of displacement, slope shear and bending 
moment are to be calculated. 

NDATA - Four digit data case identification number 

RESTRICTIONS: For the declared size of the arrays: NCO 5* 38, NBAY £ 5, 

NSUP < 6, NINT < 5. 

ACCURACY: (See NROOT and EIGEN) 

SIZE: 0011 14 a 
o 

REFERENCE: Reference 1, pp. 9-11. 
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FLOWCHART: SUBPROGRAM SSBM (NC0, NBAY, NSUP, IBL, NINT, NDATA) 
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Appendix C 



STEP 10.10 

STEP 10.11 
STEP 10.12 

STEP 10.13 

STEP 10.14 
STEP 10.15 
STEP 10.16 

STEP 10.17 
STEP 10.18 


FL0W CHART: SUBPR0GRAM SSBM (NC0, NBAY, NSUP, IBL, NINT, NDATA) 
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FL0W CHART: SUBPR0GRAM SSBM (NC0, NBAY, NSUP, IBL, NINT, NDATA) 


97 





Appendix C 



STEP -10.26 


STEP 10.27 


STEP 10.28 


STEP 10.29 


STEP 10.30 


STEP 10.31 


STEP 10.32 


STEP 10.33 


STEP 10.34 


FL0W CHART: SUBPR0GRAM SSBM (NC0, NBAY, NSUP, IBL, N INT, NDATA) 
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: — : : ! ! -J Appendix C 

subroutine ssbmi nco . nb a y » nsup » i bl , n inj » ndata ) : . — - — 
DIMENSION EVL (38) , Eve (1156) rWB(5 ) »BLMt (5) 
dimension sug) »sc(6) »sr(6> »rl(6) »rc(6> »rr(<T) 
dim ension nel(5) >nCp(6) >si(4) »S2(4) »S3(4) | 

DIMENSION DUM1(293> *DUM3(5> 

dimension xi geo) >W(ibq ) >pwu80 ) *shr<i 80 ) >3nd(i8Q) : 
Common EVLrEVC,DUMl,WB»BLMT#BW,DUM2 
CoMMONTL,SL,SC>SR>RL>RC>rR>NEL>NCP 
Common DUM3rSl»S2,S3,Xl*W,Dw>SHR»BND 
integer C nt 

l0 = 6 , 

StP-,1 » /FLqAT (NINT ) . 

NlNT=NINT4l 



; nn=nco+i-i 

' WrITE(IQ>102) EVL(Nn) 

NFi-=NN*NCO-’i 

0M=6.28316*£VL(NN) 

qms=om*om 

NcUR=0 

■! ea-o 

1 Do 921 J= 1 »NBAY 

NlMT=NEL(J) 

■ NlMaX=NLMAX+NLMT 
• ' BMM=OMS*KB( J)'/386. 

; 1F(J-1) 9Q1>901»9Q7 

901 IF ( I8L-1 ) 907 » 902 » 902 

, 902 Ki=NCO*(NN-l)+l 

K2=K1+1 

Dl=Q. 

02=0. 

D3=EVC(K1) 

04=EVC(K2) 

CA-1 

call elem(J»i»si»s2,S3) 

ShRI=S2(1)*TL*D3-*-s2(3)>04 

BImUI =-S2 ( 2 ) ♦TL+D3-S2 ( 4 ) *04 

: Call ELEM( Jr 2 >Sl rS2>S3) 

ShRI=SHRI-OMS*S2 ( 1 ) *TL*D3-0|VS*S2 (3)*D4 

BnDI=BNPI+ 0MS*S2(2) *TL*D3-t-0MS*S2(4)*D4 

907 continue 

Xl (CNT+1)=0. 

LC=LA 

DO 920 K=i»NLMT 
iF(K-l) 903»903»90 i +. 

903 XS=X1 (CNT+1) 

Go TO 905 
’ • 904 XS=X1(CNT) 

■ 905 IfU-LA) 90a»914»90g 

’ T iSUBPROGRAM^ j 





Appendix C 908 _ 1 c=NCUR+2*K-2*LC-i " J 

■ — " — Ki=NCO*(NW-l) + IC 

K2-K1+1 
K3=Kl+2 
. ' RTT=Kl + 3 

ip'(Kl-iMFL) 909 » 910*910 
. 909 Of=EVC(Kl) 

(J2=EVC(K2) 

D3=EVCTK3) ' 

Oq-EVC ( K4 ) 

Go ¥0 911 
s 910 Dl=|:VC(Kl) 

U2=£V"C ( K2) ~ ; ■ : 

. 03=0. 

04 = 0. 

911 OG 912 KK=1»NSUP 
r F (NCP(KK)-IC) 912*913*912 

912 Continue 

J Go TO 914 

I 913 ShRP=(0MS*RL(KK)/386.-SL(kK) )*TL*D1+ 

? 2(0i v >S»RC(KK)/386.-SC(KK) )*D2+SHRT 

BnOP=(SC (KK)-OMS*RC(KK) /3b6. )*TL*D1+ 

< . 2TsU (KK ) -0"MS*Rf< ( KK ) /386 . ) *02+BNDI 

i j • Go TO 915 

914 5hHP=SHR7 
8nDP=BN0I 

. ' * 9i5 _ t'0NTTNuE “ 

I Op 920 JJ=l»NlNT ; 

Cn I =CIm1‘ + 1 1 ’ 

ZE=STP*FLOAT(JJ-l) 

■ XEZE + BCMTTOT “ ’ — — — - 

; I Xi ICNT ) =x+xs 

' ZESzZE*ZE 7 

ZEC=Zt*ZES 

WTCTTTT=TZ7>'ZEC - 3 . * ZUS+ i."n*T CTO 1 + X'=f ( Z E s - 2 . * z L + r ; ) * D 2 

i 22 eS*( 2.+ZE“3. )*TL*D3+X»ZE^(ZE-i. )*D4 

0 W ( C NTT=o7*ZE*TzE^ i . ) * TL*C1/QCMT ( J >^(37*ZES^f7*ZE+ TTT*!^ 
26 . *ZE* ( ZE-i • ) *TL*03/BL.MT (J)+ZE*(3.*ZEs-2. )*04 
5RR ( C N TT=5HRF F EMM * X * ( . 5 * ZEC^ZES^TTI^Tr+Ol + 
2 Bk.M*X*X*(,25*ZES-2.*ZE/3> + «5>*02- 

' 3BMt v i»X*ZE5* l . b*Zt-i • ) *TL*03+ " 1 

4BmM*X*X*ZE* ( .25*ZE”1 ./3« > *04 

i O|\;U(CNTr=5T^0P+SHRP*X+BMM*X*X*< •1 *ZEc-,25*zES+»5)*TL*D1 + 

' 2B Fj M*X*X*X*( .05*ZES“ZE/6. + l»/6. )*02” 

3SMM*X*x1FzeS* ( .1*Ze-»25) *TL*D3+ ~~ 

4BmM*X*X*X*ZE*( .05*ZE-1./12« ) *04 

ShRI=SHR(CNT') 

~0NDi=bNt)(CNT) 

920 Continue • 

" RcOR=TC+l 

921 continue 

' wmax=o. ' " ' 

Dy.iMAX=0 • 

. 5hRmaX=D . 

BnDMAX=0 . 

T L 1 MERCK A X * N I n T - ~ 

Do 924 IWrlrLlM 

' iFlAdb(WUW) )-WNAX> 9].6»91fa.»917 ‘ ? i 

SUBPROGRAM SSBM : CA R D IMAGE U STI N G 2/3 ; 

i 100 :. ; 



|AppendixC : 

917 'k M A X = A BS ( W C IwT) ' ' ■ ' 

9 16 i F( Aos(DW( itf) )-dwmAx) 9i8»9i8»9i9 , _J 

919 DwMAX=ABS(DW{IW) ) 

916 liF<ABS(SHR(lW) )"ShRMAX) 922>922>923 f ; 

923 SVi«MAX=ABS(SHR( IW) ) ; 

922 IF(A6S( BNP(IW) )-BnDMAX) 924 >924 >925 . • 

925 8iyDHAX=ABS(BND(IWj ) 

924 contin ue L / - 

D0\ 926 IM=1»LIM I 

W(1M)=W(IM)/WMAX . ' 

DwCIM)=DW(IM)/DWMaX : 

Shrt(lM)=S HR(lM)/SHRMAX 

BN®( IM)=BND(IM)/BnDMAX 

926 Cq#INUE 

LlWfl . 

DOi 928 0=1 » NBAY 
— Wi^lTE ( 10 #998) 

W[|ITE (I0»997) J ; 

LlWI=LIM+NlNT*NEL( U)-i 

Up 927 IRsLlMfLIMM ' ' 

write ( io » 999 ) xkirj »wdR> »dw<ir> »shr{IR) ,bnd<ir> 

q?7 Continue 

L1M=LIM+NINT*NEL(J) 

928 Continue . 

929 continue ~ ’ 

100 FORMAT ( 1H1 ) 

101 FoRMAT(25X»9HOATA CaSE»I4) 

102 FqRMAT(/»SX>HHFReQUEnCY =>E12.5»1X>3HHZ. ) 

997 FoRMAT(/»29X»3HBAr» 12) 

998 FqRMAT(/» 9X»1HX» lnX»lHW»8X»SHDW/DX»9Xt iHVrllXrlHM) 

999 FoRMAT(5X,F7.3f2X,4E12.5) 

return • 

end 
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SEGMENT LEVEL 



TWO-DIMENSIONAL PANEL ANALYSIS: 
PROGRAM CORE-OVERLAY STRUCTURE 
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PURPOSE: 


SUBPROGRAMS 

REQUIRED: 

INPUT DATA: 
VARIABLES: 


Appendix C 


PROGRAM PLTVIB (MAIN) 

(See Figure 2 for Flow Chart) 

Main program for computing the natural frequencies and normal 
mode shapes of a nine-bay orthogonally stiffened flat panel. The 
structure i$ assumed to have clamped edges along the edges and 
four orthogonal stiffeners dividing the uniform cover sheet into nine 
bays. The stiffeners are modeled with a finite element representa- 
tion of a thin-walled open-section beam . The stiffener element 
deformation fully conforms to the plate element deformation along 
the edges of the plate element. The rectangular plate bending ele- 
ment used is based upon a sixteen degree-of-freedom element with 
an interior element mode in the form of a clamped-clamped beam 
fundamental mode. 

RDNWRT, PLTSTF, PLTMAS, RIB, NONDIM, ASSYR, ASSYP, 

FILL, ZERO, ORDER, DELETE*, LOC*, EIGEN*, NROOT* 

See Appendix D for description - 

NCASE, A(l), B(l), EP, HP, PR, RHOP, ER(I), GR(I), RHO(I), 
SR2(I), C2(l), SR3(I), C3(l), AR(I), 

A22(l), A23(l), A33(l), SJ(I), RE2(I) 

RE3(I), GM(I) 

NDEL(D) - Constraint vector for an N degree-of-freedom system. 

If NDEL(R) = 1 the R coordinate is removed from the equations of 

th 

motion. If NDEL(R) = 0, the R coordinate is ignored in subpro- 
gram ORDER. 

SI (I), S2(l), S3(l) - The arrays represented the partitioned stiffness 
or mass matrices of the plate on rib element. 

SC(D) - The coupling term between the edge displacements and the 
generalized coordinate for the plate stiffness and mass matrices. 

SF(I),RF(I) - Structure stiffness and mass matrices 
(symmetric storage mode: LOC) 

EVL(I) - Array of eigenvalues 
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VARIABLES: 

(Continued) 

RESTRICTIONS: 

ACCURACY: 

SIZE: 006371 
REFERENCES: 

*Subroutine i 


EVC(I) - Array of eigenvectors 
R(l) - Dummy array 

M - Number of grid points in f-direction 
N - Number of grid points in y-direction 

For the declared size of the arrays: M S 4, N < 4 

Not Applicable 


Reference 1, pp. 11-16 

identical to that described for the one-dimensional analysis. 
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Appe ndi x C C 

SeGzERO ( SEGo t PLTV IB . LoC r ZERO ) ' ' "'P* 

segment (segi,rdnwrT) . j 

Segment (seG 2 »pltstF) ! 

segment (seG 3 »assyp) 
segment ( SEGt* , pLtmas ) 

SeGMENT<SEG5,RIB*aSSYR> 

segment* seg 6 # nond i m , order , f i ll » delete ) 

S egment <SE G 7rNR00 T» eigen) 

PR0GRAM4PLTvlB(INRUXr0UTPuT»TApE5=lNPuTrTAPE6=0UT(3UX)__l_ 

I . Common ER(4) fGR(4) »RH0(4)»SR2(4) »C2(4) rSR3<4) »C3<4> » j 

IAr.C.4) #X 22 (-4.)^.A23T 4L,.Ao3 l4T_».SJX4 L»-RE2X4.)^nRE3_(-4-)-»~GMx4J ; 

common R(27oD *a< 3 ) »b(3) ! 

1 XoMMON_NCASE»EP»HR»RR # RHOR»_DPL»TLfTK»XM-»NCO j 

• dimension ndel( 73) »si( 36 ) #S 2 ( 64 ) »s3(36> »sc(i7) » 

lSF(625)_fRF*625)rEvLX2s)-»EvC4.625L»ANAME(2) 

equivalence (R{326) »rf(1) ) * <R<952) #evl(1) ) * (R(978> »evc(D >; 

' data_lf.ile/ 3 HLGiOZ ; » 

IoATA=l 

, 100 ReAD15»105) NDATA . .. . 

105 FoRMAT(5X»l3) 

■ 110 Mr 4 ' . . " ■ . . : 

r~ n =4 ~ 

) NcP=4 

; MbAY=M-1 1 I 

NbAY=N-1 . i 

N G P =M I* I N I 

Nul=NCP*NGP 

NcOzNUl +MBA Y +.NBAX 1 

AnAmE(D=4HSEG1 ; 

An ame .(2 ) ^0 

Call segmentilfile* uaname) •' 

CaLl-Xdnwrxu.) . 

c Initialize constraint vector 

Do— 120 - -Is.l-f.64- 

IFU-20) 115* 115*111 

_1 11 I FJ.I -2.9-)_U.4.f-l 1.5 ill 2 _■ 

112 lFd-36) 115* 115* 113 

-1-1-3 —I F-( I “45 ) -1.1.4-»Xl5.f- U5 1 ’ 

114 NoEl(I)= 0 

_ GO TO. 120 4 

f 115 NoEL(I)= 1 

: 1 20 l_C o m T I NU.E. : 

TL= SQRT(A(2)*A(2)+B(2)*B(2) ) j 

Dp=Ee*HP*HP*BP/LU^^J4^Efit^EBJL) _J. 

call rdnwrt ( 2) 

Rh0P=RH0P*HP/.386 . 

Nul=NCP*NGP 

Do_.305_IQP=lt2_ ; 

ncOzNui+mbay*nbay 

C&LL ZERO(R»?.701 ) : ' 

IF(IOP-I) 200 *200 » 210 

C _-CoMRUJ-E_AND_ASSEM 6 LEXFREE*ErEE_P.LaIE_S-T-IEFNESS_MAt-RIX~— -- 

• 200. AnAME11>=4HSEG2 

• — - — — AnAmE-( 2J-=0 ■ — - . 

P ~ PROGRAM PLTVIB: CARD IMAGE LISTING 1/3! 
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App endix C 


' call segment (Lfile»i»‘aname> - " 

Aim A MELUJj: 4.HSEG3 : , 

AnAME(2)=0 

CaLI — :SEGMeN.T_(LEILE-#-?-#-ANAM£.) : 

Do 205 J=1»NBAY 

t Dy=B(J7 : - 

Do 205 1 = 1 > MBAY 

D.x=A..LUJ : l 

call pltstF(dx»dy,dp,pr,si>s2»s3»sC) 

■ C a LL-A SSYPJ R»Sl»_S2^3,-SC»I»M»si.»NfNCpj 

: 205 CoNTlNup: 

Go TO. .220 : : 

c compute and assemble free^free plate mass matrix 
_210 AnAME.U.L=4JHSEM 

<; A NAME < 2 ) = p . 

i . Call segment (LFileh*aname> 

A|mAME( 1 )=uHqEG3 : 

AimAME ( 2 ) =0 

, C. A L 1 S EG.MENXLL EJ.L.Ex2j- A NAME.) 

Do 215 J=1#NBAY 

; D.y=B-(-JJ_ - 

Do 215 1=1 f MBAY 

‘ Dx=A.(.IJ 

i Call PLTMAS(DXfDY,RH0p»Sl»S2,S3»SC> 

C. A L L__AS SY.E-(.R.»-SX*J5^.tS3.#-S.C.»Jt-fJ4.f_Jj-N.«LNC PJ. 


i 


! 215 continue 

-CoMpUTE „ AND„ ASSEMbLEL_FREE-£rEE_R .I.B_SXiEEMeSS_AND_M ASS_MAIr.I_CES- 
220 AnAmE(1)=4HS£G5 

AfjAME;(2)=0 : 

c a ll segment (LFile» i » aname) , 

! Do_27.0__I.8=l.f_4. : 

Go TO ( 225 » 235 » 250 » 255) »lR | 

' -22-5 U=2 : J 

Go TO 240 

—23 5 U=5 - — — 

240 A X = 1 • 0 

Do -2-4.5_I=1 jlMBAY 1 1 

Dr=A(I) 

c A L| R IBXI R»_IOP_iDr_»j51-,-S2»-S 3 ) : 

Call ASSYR(R»S1»S2»S3. I»M.r.U»NrNCP,AX) 

—245 CONTINUE— : 1 

Go TO 270 - 1 

-250 1-2 : ! 

Go TO 260 

_255 — 1=3 i 

260 Ax=0 • 0 

L D.O— 26.5_J=l_f_NBA Y ; : 

Dr=B(I) ' 

'■j CaLl_RI BXlR#..IGP-».DR^l.,.S2rS 3 ) __J 

call assyr(r,si»S2>S3, i»m» j,n»ncp,ax) 

—265— Continue - — — — — — — i 

270 . Continue i 

j: ; N 0 ND.I.ME N SION A L IZE_M ATrIC E S : 

•, T R = 0 .0 ; ; • _______ : - 

Program pltvib: card image usting 2/3: 
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* . . 

r 'I Do ?75 I=3>NU1 »NC p T — 

Il=I-l 

1 1 1 1=1 1 +J 1 1*11-1 1 )/2 : 

Il=I+(I*I-I)/2 

Tr=.TR+R ( II1)+R(II) ; 

i 275 continue 

A N A M E .(. U-=MH S E.G.6 

AjvjAME ( 2 > = 0 

: C a LL- -SEGMENTS LEI Le NAME) , 

lF<IOP-l\) 280 *280 » 285 


_2fl0 Tk=TR/-( F.LOA.U.MJ *ElO.AT-( NU 

CaLL NONiV>lM(R*M*N,TK*TL) 

!_ GO-TO_2.90 


285 

Tn;=TR/( float (M)*FlOAT(N) ) 
CflLL NONDIM(RpM»N,TM.tL) 

» 

j 

: 290 
C 

NqL=0 

Apply coordinate constraints 


1 CaLl ORDER(RfNDEL,NCP,M»N,NDL) 

NrOrNCO-NDL 


IF(IOP-I) 295*295,300 


295 

Call fillisf*R*nco*o> 

Go- TO. 305 . 


i 300 
— 305— 
C 

Call fill(Rf»r»nCo*o) 
Continue 


compute eigenvalues and eigenvectors 

AmAME ( 1 ) =4HsFGl . 


Af.jAME (2 ) =0 

CaLL SFGMFNT(LFILe»1pANAMf> 

• C A LL RONWRT (3) 

AmAME ( 1 )=4H9FG7 , 1 

AnAME<2)=0 i 

C A LL SEGMFNTfLFILF# ^ .aNAMf) < 

CaLl NROOT(NCO,SF,RF,eVL»EVC> 

i _ AmAME(1)=4HREG1 ! 

AkAME(2)=0 

! CaLI SEGMENT (LFILF* X »aNAME) 

Call rdnwrt (4) 

310 

IoAtA=IDATA+1 

Ip(NDATA-lDATA) 3l0»li0»H0 

Continue 

. . ... 


End 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 001026 g 
REFERENCES: 


SUBROUTINE RDNWRT (IT0) 

This subprogram contains all input/output statements for program 
PLTVIB. Input data definition and format is described in Appendix 
D. An example of output format is also included. 


None 

IT0 = 1, input data read 

= 2, edited input data is printed 
= 3, NCASE, TK, TL, TM are printed 
= 4, eigenvalues and eigenvectors are printed 

NCASE Four digit data identification number 

TK,TL,TM Nondimensionalizing constants for force (stiffness), 
length, and mass, respectively 

1 < IT0 ^ 4 

Not Applicable 


None 
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FLOW CHART: SUBROUTINE RDNWRT(ITO) 


109 




Appendix C 

subroutine rdnwrt (itoj 

Common ER(4) *GR{4) »RHo(4 ) #Sr2(4) »C2(4) » SR 3 ( 4 ) > 03 ( 4 ) * 

1Ar( 4) »A22(4) »A23(4) *A,33l4) *SU(4) * RE2 ( 4 ) , Re3 ( 4 ) ,GM<4) 

Common R(27od *a<3) ,b(3) 

common nc ase » ep » hp * pr , rhop * dp * tl ► tk * tm * nco 

Dimension Evl< 25) ,Eyc(b25 ) *rF(625) 

equivalence {R<326) *rf<d j » <R<952) »evi ( i) ) * {R(978 > »evc(D > 

In=5 

l0=6 

60 TO (100»200»400»50n) » ITO 

mo RfAD(lN* 600 ) NCASe 

RfAD(IN»605) A(l) f A(2) ,A(3> 

READ(IN»605) BCD *B(2) »B<3) 

READ(lN»6lO) EP * Hp 9 PR t RHOp 
Do 105 I=l»4 

REAO(lN»6lO) ER ( I ) » GR ( I ) » RHo ( I > 

RlAD(IN*6iO) SR2 ( I ) »C2(I)»SR3(I) *C3(I) 

Read UN f 610) ARd) *A22(I> ,A23(I) *A33 (i) 

Read <in» 61 .O) SJti) » rE2(I) ,Re3(i> *gmu) 
ins continue 
return 

200 WrITE( I 0»b20 ) NCAsE 
WrITE(I0»62S) 

WrITE(I0»630) A ( 1 ) »A(?) »A(3) 

WRITE ( 10*635) B(l) *B(p)»B(3) 

Wp ITE ( 10 * 640 ) 

WrITE( I 0r645) EP » PR » Hp » RHqP * DP 
WplT£(IO» b 50) 

DC 140 1=1 #4 

Go TO (130*130*135*135) *1 
130 WrITE(I0»655) 

WrITE(IO»660) I » Er ( I ) , GR ( I > * RHO ( I ) 

WrITE(I 0*665) SR2(I) »C2(I) »SR3(I) »C3{j) 

WrITE(IO» 670) AR(I) rA22(I) »A23(I) »A33(I) 

WrITE (10*675) SJ(D ,Re2(I) >RE3(I) #GM(i) 

GO TO 140 
135 WrITE< 10*680) 

WrITE(IO»660) I * Er ( I ) , OR ( I ) * RHO ( I ) 

WrITE ( 10*685) SR2(I) rC2(I) rSR3(I) >C3(i) 

Wr ITE ( 10*690) AR(I) ,Ar 2(I) »A23(I) r A33 < I ) 

WRITE (I0*fe95) SJ<i) »RE2CI)*RE3(I) »GM(i) 
mo Continue 
return 

400 Wr I TE( 10*620) NCAsE 
Wr ITE ( 10* 700 ) 

write ( 10 * 705 ) tk*tl*tm 
return 

500 Do 415 I=1*NC0 

IF(EVLID) 405*410*41(1 
405 EyL(I)= ABS(EVL(I) ) 

410 FrEq= 0.159155*SQrT(Tk*EVl(I)/TM) 

WrITE(I0»710) EVL ( I ) » FREQ 
WrITE(I0*715) 

Ix=NCO*(I-l)+l 
I Y-NCO* I 

Wr I TE ( 10* 720 ) (EVC(J) ,U=IX»IY) 

415 Continue 
600 FoRMAT(5X*l4) 

SUBPROGRAM RDNWRT: CARD IMAGE LISTING 1/2 
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! 605 FoRMAT(4X»3£12.5> 

_ 61 0 __. .FoRMAT <4X..4E12. 51 

C615 F 0 RMAT(4X#4I2) 

1 620 F 0 R M A J. (.1 Hi .» 2 fl X_f.-9.HD A J A_.C A SE.r 15 ) 1 

625 FoRMAT(/f 28Xf 19HFREE VIBRATION OF A,/,19X, 

l 137HNINE. BAY',. ORTHOGONALLY. STIFFENED _PANELf.// # 28X # 

21qHSTRUCTuR^L GEOMETRY*/) 

6 3 0_.F.O RM AT_( 15X» ^HA 1= » El2.5t2X»3HA2=».El2.5,2Xf3HA3= ».El2..5 * / ) 

i 635 F 0 RMAT(15Xf ^HBl=»El2.5»2X,3HB2=»El2.5,2Xf3HB3=»El2.5*//) 

_64.0 FoRMA _T_(_3 ftXJLl 6 H C 0 Y E 8. _S HE E I_D A T A *./.)_ 

\ 645 Format (5X»16HY0UNg»S M0DULUS=rE12,5»l5X»l6HP0lSS0N , S RATIo=# 

1 lEi2.5»/f5Xri6HTHlcKNEsS = f.El2.5 fj5Xrl6HWE IGHj/ VOLUME = », 

2Ei2.5r/»24Xfl7HBENDlNG RlGlDlTY=fEl2.5»//) 

■_,65o format ( 3ix h4Hsti ff fnfr oata*/) 

655 F 0 RMAT(24X»29HSTIFFENERS PARALLEL TO x-AXlS) 

_66Q FoRMAT,(.5X.fl3HSJ_IFpEN£R^a._Ll2jL2X.L4.HE =»Ej2_. 5»2XJ4.HG - » El2..J5j 

f 12x » 6HRH0 -f El2«5) 

L.665 FoRMATC 5X , 3 hSY= ».El2 . 5 # 2X ■» 4HC Y_= »_El2.. 5 , 2 Xf 4 HSZ=fEl 2 , 5 » 2X » 

j 16hCZ =»E12.5) 

1—670 Format ( 5X ,3hA.= ,Ei 2 # 5 #2X»4HIYY= » El2 .5 , 2X » 4HI YZ= »_Ei2 . 5 * 2X » 

i 16 hIZZ =fEl2.5) 

,—675 F.0RMATX5X ».3h J_.= ».E l2.5^2Xf-4HREY=.r El2 .5 ,2XfJ;HREZ=»_Ei2 . 5 r 2X» 

16HGAMMA=»El2.5f/> 

_68o . Format (24x»29HstifFeners parallel t.o_ . y-axiS) 

1 685 FoRMAT(5Xf3HSX=»Ei2.5f2X»4HCX =»Ei2.5»2Xf 4 H SZ =»Ei2.5»2Xr 

16hCZ =fEl2.5) 

690 FoRMAT(5X*3HA =r E a 2 . 5 . 2 X » 4HlXX=f El2.5,2Xr ^HjXZZf Ei2.5»2X* 

1.6HIZZ =r.El2-.J5J 

695 F 0 RMAT(5X,3HJ =r Ei2.5»2Xf4HREX=f El2.5,2Xr 4HREZ=»Ei2.5f 2X» 

A_6hQ AMM A-_f_E 1 2 »_5 »./ ) ] 

700 format ( / » 23x » 29hnond i mens i onal izing Constants ) 

—70 5 F 0 R M A T_(. 1 3 X 3HJ K = ». E A 2 • 5 » 2 X , 3hTL= f El 2 . 5 , 2X f 3 H JM= r El 2 • 5 ) : 

: 710 FORMAT (/»5X,llH£lGENVALUE=rEl2.5fl9Xfi0HFREQUENCY=fE12.5f4Hf HZ.) 

—7 15 FORMAT (33X» IAHEIGeNVE CTOR) 

720 FoRMAT(4X»El2.5»2X»El2.5»2XrE12.5»2X»El2.5»2X»EA2.5) 

Return ■ ■ • 

- — - flENOj - — -■ 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 002220g 
REFERENCES: 


SUBROUTINE PLTSTF (A,B,D, P, SI , S2, S3, SC) 

To compute the stiffness matrix of the modified sixteen degree-of- 
freedom plate element described by Bogner, Fox, and Schmidt with 
an interior mode in the form of clamped-clamped beam functions. 


None 

A - Dimension of the plate element in the x-direction 
B - Dimension of the plate element in the y-direction 
D - Bending rigidity of plate element (Reference 1) 

P - Poisson's ratio for the plate element 

S1,S2,S3,SC - Partitioned stiffness matrices for the plate element 
(see Flow Chart and Appendix B) 

A^O, B/0 

See References below 


o Reference 1, pp. 12-15 

9 Bogner, F. K.; Fox, R. L.; and Schmidt, L. R., Jr.: Confer- 
ence on Matrix Methods in Structural Mechanics, Wright- 
Patterson Air Force Base, Ohio. AFFDL-TR-66-80, 1966, 
pp. 391 - 443 

o Mason, V.: Rectangular Finite Elements for Analysis of Plate 
Vibrations. J. Sound Vib., Vol. 7, No. 3, May 1968, 
pp. 437 - 448. 
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FLOW CHART: 
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SUBROUTINE PLTSTF ( A / B / D,P,S1,S2 / S3 / SC) 
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SUBROUTINE PLTSTF ( A , B . D # P-» Sl » S2 tS3i.SC) ____ ' 

dimension sk2) »S2<2> ,s3(2> »sc<2> 

_AB=A*B _ ; 

F= B*B/(A*A) 

[ G= A*A/(,B*B) 

‘ R- F + G '= 

L _CO= .39647605 

Ci 1= .20771538 

C2l=_. 04^69616 : 

C3l=4. 64712757 

PhI=R+2.*C11*C11*(C31-2.)»(C31-2.) 

SMS CO*CO*PHl/(C2l*C2l) 

: : SM 1 z S MZ2 0 4 8. ,.-C 0 * C_i 1 * C i_i 

SM2=SM/40g6.-C0*CiltCi 1 

S lJl.).- 15.6.* R / 35 ,±2... fl 8+.» 0 6 2 5* S.M 

: i Sl(2>= ( (22«*F+78.*G)/35. + 1.2*( .2-*-P) + .0l5625*SM)*B 

• SlJ3> = ((4.*F+52.*G)/35. + .32+.00390625 *Sm)*B*B 

Si ( 4 )= -( (78.*F+22.*G)/35.-H.2»( .2+P)+.0l5625*SM)»A 

: Sl < 5 ) 1 . *R / 3 5. + .0 2+ 1 . 2*P+ . 003g0625*SM) *AB 

Si(6>= ( (52.*F+4,*G) /35.+.32+.0Q3906?5*SM>*A*A 

Sl_(-7-)_= (ll ._*. R /_3 5^-E.O 2+.2 * P.+ 8 • *SMl ) *Ar. 

i 1 Si(8>= ( (2.*F+22.*G/3.)/35.+2.*( .2+P)/15.+2.*SM1 )*aB*B 

; S i < 9 > = . - U 22 . *F/ 3 . 1 2 . * G ) / 3 5 . + 2 . * < . 2 4P)/1 5.±2.* S M 1)* A* A B - 

t- Si(lO)= ( 4 . ♦R/105, + 8 » /225 ,+SM2) *AB*AB 

L Si (. 1 1 ) z- ( 156. ♦F-54. *G ) /35 .-2. 88t . 0625*SM 

Si (l2)=-( (22.*F-27.*G)/35.+.24*{ .2+P)-.0l5625*SM)*B 

SlJ.l3)=_((7a.-*F-l3»»G4y-35. + .^4-.„0l5625*SM)*A 

Si ( 14 ) =- ( ( 1 1 . *F-6 . 5*G ) /35 . + . 1* ( . 2+P ) -A . *SMl ) *AB 

Si d 5»=_s Ld) : 

Si(l6)= Sld2) 

Sid 7 ).=- (_( 4 . *F-18 . *G ) /35 . + . 32- . 0 0390625*SM ) *B*B 

Si (18)= ( (ll.*F-6.5*G)/35.4.1»( .2+P)-,003g0625*SM)*AB 

S 1. d 9 > =- ( ( 2_. * F-13.*G/^.)/35 .+2 ,/7 5.-2 T *SMi ) *AB*B 

! Si (20 ) z Sl(2) 

Sl<21A=_Sl(3)_ i 

Sl(22)=-Sl(13) 

S ll2_3_ L=-S l < 1 8_> . _ 

| Si(24)= ( (26.*F-3,*G)/35.-.08-.0b390625*SM)*A*A 

51 (25)=-(<ll. *F/3.-l.S» G)/35 .-( .2-»P) /nQ«~P.*SMl ) *a*AB 

Sl(26)z-Sl(4) 

SlJ 27 .) l=-SU5 ) . 

Si ( 28) z Sl(6) 

SlJ 29)=-S 1.(14) ; 

Si (30)=-Sl (19) 

SlJ31)-=_Sl(25) 

Si ( 32) z ( (2.*F/3.-G)/35.-2./225.-SM2 )*AB*aB 

Si(33)=-S1(7) : , 

Si (34)=-Sl(8) 

SlJ 3 5 )z_S 1(9.) 

Sl(36)= SKlO) 

S 2 (l).z .(.54 «.*F-156 ».*G )./ 35 . “2 • 88+ » 0625*5M 1 

S2<2)= ( (l3.*F-7a.*G)/35.-.24-f.0l5625*SM)*B 

S2 (.3)_=_-( (27.*F-22.-*G)/35.-1.2*( .2+P) + .0l5625^SMl*A__ 

52 (4 ) = ( (6.5*F-1i.*G)/35.-.1*( .2+P) +«. *SN|1 > *AB 

LS2J5Az_=.,54..+R/35.+2.8fl+.0625*SM 

S2(6)= -( (l3.*F+27.^G)/35.-.24-.0l5625*SM)*B 

S2J-7J_z_=J_( 2? .♦F+l 3 .*G)y 35.. - . 24-.0l5625*SM)_*A_ 

S2<8)= <6*5*R/35.-,02-8.*SmD*AB 

’SUBPROGRAM PLTSTF: CARD IMAGE LISTING 1/3 


S2<9>= 1 
_S2<10)=- 
S 2 (ll)= 

_S2ll2)_=: 
S2U3)=' 
._52<l4) = 
Sp<l5\) = 
_S 2 C16)_=j 
5 2 <l7) = 
_S2il8.)=j 
52(19) = 
_Sp.(20.)=' 
S 2 ( 2 1)= 
_S2 (22)= 
S 2 ( 2 3)= 
_SpJ.24)J= 
S 2 <25)=' 
_.S 2 _( 26 )= 
S 2 (27)= 

...S 2 ( 2 a) = 

S 2 ( 2 9)= 
_S2(.3Q).=j 
S 2 ( 31 ) = 
_ S 2 ( 32 ) = 
S 2 (33)= 
._S2i34) = 
S 2 ( 35 ) = 
_5pJ 36) =: 
S 2 (37)= 
__S2(38) = 
Sp(39)= 
_S2(40)=' 

S 2 ( 41 ) = 
SpJ 42) = 


■S2(2) A P P endlx C 

•( (-3,*F-26.*G) /35*-+2./25. + . 003p0625*SM) *B*Q 

( (6.5*F-H.*G)/35.-.l*(.2-»-P) + .00390625*SM)*AB . 
■(( 1. 5*F-ll»*G/3«)/35»+_(.*2+P)/.-50.#t2 • *SMH*aB.*B.J 
■S2(6) 

( ( 3. *F+9 •*G)/35.+. Oar ..00390625*SM) *B*B ^ 

(6»5*R/35.-,02-»00390625*SM)*aB 

■1(1 . 5*F+l3.>G/6 .)/35 9 +1./150.-.2 . *sMl L*AB*B_ 

S2(3) 

;S2(JLi) i j 

( (l8.*F-4.*G)/35.-.32-«-.00390625*SM)*A*A 

-1(1 3. ♦F/3.-2.+G ) /35 .-2./I5 • +2 . *.SM i_U A* AB 

*S2 ( 7 ) 

_S2ll5l_ ; 


5p(43) 

S 2 (_44 )_ 

S 2 ( 45 ) : 

5 2( 46)= 52(10) 


( (9,*F+3« *G) /35, + , 08“*00390625*SM) *A*A 
^1(13 »*EZ6.»±.l . ft».GJV^5 ♦S.MlJ,* AMB- 

•S2(4) 

52(12) 

■S2 ( 20 ) 

-( (F -2,*G/3. )/35.+2./225.+SM2)*AB»aB 

S2 ( 8) 

-52 ( 16) 

•S 2 (24) 

“ ( R / 7 0.jrl_t^L45 0. -5M 2 1*A B*A B 

S2(5) 

52(6) 

-S2(7) 

rS2(8) : 

S2(l) 

52.(2) 

•S2(3) 

■5214) ■ 

52(13) 

_S2 ( 14 ) ■ • 

■S2 ( 15) 

■S2J 16.1 ; 

52(9) 


S 2 (47): 
S 2 (48)_ 
Sp ( 49) 
S 2 (50) 
S 2 ( 51 ) 
S 2 (52 ). 
Sp ( 53) 
.52(541 
S 2 ( 55 ) 
52156). 
S 2 ( 57 ) 
_S?,l5.fl.)_ 
S 2 ( 59 ) 
S 2 (60L 
S 2 ( 6l ) 
S 2 162 L 
52(63) 
_S2l64J 

s 3 (i)= 

53l2l=_ 


-S2(U) 
-S2(12)_ 
-S2(35) 
52(431 
52(23) 
_S2(.24l 
-S2(17) 
“52118 1 
S2 ( 19 ) 
52(201 
-52(29) 
-521301 
52(31) 
.521321 
“52(25) 
“52(261 
S2 ( 27 ) 
521281 


SKl) 

-5JL121 


S3(3)=._S1(3) 

• SUBPROGRAM PLTSTF: CARD IMAGE LISTING 2/3 
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S3<4>= Sl<4) 

Appendix C S3<5) = _tS1<5> 

■ S3<6)= Si ( 6) 

S 3 (-7J _=_=S 1 ( 7) 

\ S3<8)= Sl(8) 

S3<9) = --Sl(9) : - 

S3 (10)= Si (10) 

s 3d i ) =_s i d i ) 

S3d2)=-Sld2> 

S3d3-).=_si-d3)__ : 1 

S3d4\)=-Sld4> 

S3 ( 1 51=_S Id 5 ) 

S3d6^=-Sld6) 

; S3 d 7|=_S 1 (_l n 

S3da#=-SK18) 

S3d.9l)=_Sl.d9) 

S 3 ( 20 ) ="Sl (20 ) 

S3 (21 )_=_S.l ( 21 ) : 

S 3 ( 22 ) = Sl(22) 

S 3 ( 23 ) =-SU23 ) 

S 3 ( 24 ) = Si (24> 

S3.( 25 ).=Sl< 25 ) : 1 

S3«2fa)= Si ( 26) 

— S 3 (27 )=rSl( 27 ) 

S3(28)= Si (28) 

S 3 (29) =-Sl.( 29) 

S 3 ( 30 ) = Si ( 30 ) 

5.3. (.3 1 )_=.-Si_(,3 1 ) 

S 3 ( 32 ) = Si ( 32) 

S 3 ( 33.) j=rSl_< 33 ) 

S 3 ( 34 ) = si ( 34 ) 

S3l3S.).=^SlJ.35.) 

j 53 ( 36 )= Si ( 36) 

Sc ( 1 ) = - « 25 * S.M 

SC ( 2 ) = -.0625*SM*3 

SC-(.3)-= .0625*SM*A — — — 

; Sc(4)= -32»*SM1 *Ab 

SC (5)= ~.25»SM ' 

SC(6)= -.0625*SM*b 

Sc ( 7 ) --5C ( 3 )_ 

SC(8)=-SC(4) 

SC.(-9. )_" -»25*SM 

SC(10)=-SC(2) 

S.CJ 1 1..) j=_. . 0.6 2.5 * SM *A _ 

Sc(l2)= 32»*SM1 *Ab 

S cd 3 ).=-.» .2 5.+.SM 

SC(14)= • 0625*SM*B 

SC-(-l 5 )_=.-• 0625 * S M* A 

Sc(16)=-32.*SM1*Ab 

■ Scd7.)-=-;SM 

Do 25 1 = 1 » 64 

lF-(I-17J-5»5fl0 

5 Sc(I)= D*SC ( I ) /AB 

_10 I EXI-36 ) 1.5 , 15.f.20 

15 Si(I)= D*S1(I)/AB 

s 3 _d)-=_D-* s 3.(j )/a a : 

20 S2d)= D*S2(I)/AB 

_25 -C 0 N.T..I NUE. 

_ return 

i EN D ^ : 

SUBPROGRAM PLTSTF: CARD IMAGE LISTING 3/3 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 001520g 
REFERENCES: 


SUBROUTINE PLTMAS (A, B,R, Si , S2, S3, SC) 

To compute the consistent mass matrix of the modified sixteen 
degree-of-freedom plate element described by Bogner, Fox, and 
•^Schmidt with an interior mode in the form of clamped-clamped 
beam functions. 


None 

A - Dimension of the plate element in the x-direction 
B - Dimension of the plate element in the y-direction 
R - Mass per unit area of plate element 

SI, S2, S3, SC - Partitioned mass matrices for the plate element 
(see Flow Chart and Appendix B) 


None 

See References below 


See References for subroutine PLTSTF 


117 



Appendix C 


*SEE APPENDIX B 



FLOW CHART: SUBROUTINE PLTMAS (A,B,R,S1,S2,S3,SC) 
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SUBROUT.I NE RLTMAS ( A ,B » R»Sl »S2rS3*S.Cj Appendix C ' 

DIMENSION Sl(2) rS2(2) ,S3(2> »SC(2) 1 — 

AR=A*6 : 

Ci= 13./35. 

; C2= J11./210. 1 

C3= s 1./105, 

„.C4=.J_9 . /_70 , 

C5= 13./420. 

C Ll../ 1.4.Q_. 

Co= 0139647605 

C 1 1 = J, 2077 1538 

Cpl= 104469616 

L__ C3l=4 1),64-72757_ : 

Rl= Cd|M .0625*00-2. *CH*Cll) 

R 2=_O0 *X. 0 l5625.*0.Q-Cl.l*.(.C 2 1.+. .2 5* C_ll.) J 

R3= C0*( .00390625*Co-.5*Cll*C2l) 

R 4=_ C 0*J_..0 039 0.62 5*0.0 -C2 1 * C 2 1^.062 5.*CJL-lACXU 

R«S= Co*(CO/1024.-,25*C21*(C2l4.25*CH) ) 

86=__c o *J.C o / 4096.^.1 2 5 *c 2 l*C 2 l.) 

R7= C0*( .25*00-4. *Cll*Cll ) 

Ra=^C 0*J.,062 5*00-4. *-Cl-I.*02l) 

R 9 = Co*(. 015625*00-4. *021*021) 

Sldl= C 1 *C 1 +R 1 

Si(2)= (C1*C2+R2)*B 

; Sll3)= ( Cl*C3+R3j* 3*B 

S"i(4)= - (C1*C2+R2) *A 

; S i.(.5.)_=_=_(.C 2*0 2+R 3J_*AB 

Si(6)= (C1*C3+R3) *A*A 

‘ Sil7_).= ( C 2 * C 2±B_4 ).* A B 

Si(8)= (C2*C3+R5) *AB*B 

S i_(_9 )= ~(C 2*0 3±R 5J * A * AB 

Si (10)= (C3*C3+R6) *AB*AB 

Si tl 1J. =_ 0.1 *C4±R i_ 

Si(l2)= (C2*C4+R2) *B 

1 S i (13) = -1.0 1 *05 ± R 2.)_* A 

Si ( 14) = (C2*C5*R4) *AB 

■: Sid 5. ) =_s ij i ) 

Si(l6)= Sld2) 

Si d 7_)_=_J 0 3*c 4_+R3 )_*3.*B • 

Si (l’8)=-(C2*C5+R3)*AB 

Sill 9. )_= (. C_3 * 0.5 +B_5_)_*A B* B 

Si(20)= Sl(2) 

S l J 2 11= -SI (3 ) 

Si (22) =-Sl ( 13) 

Si_( 2 31 = -S l_(.l 8-)_ 

Si (24)=-(Cl*C6+R3)*A*A 

Si_( 2 51=_( C 2 * C 6-+.R5-)-*A* A B 

Si ( 26 ) =-Sl (4) 

S 1 .( 27.) =-S 1 (.5 ) 

Si (26)= Si (6) 

‘ S_i_(.29.)=-S 111 4 J 

Si ( 30 ) =-Sl (19) 

• S 1(31 )= .Si ( 25 )- r 

I. Si (32)=-(C3*C6+R6)*AB*AB 

Sj.133 ) =-S 1 .( 71 

Si ( 34 ) =-Sl ( 8) 

! Sl_(35J.=_3jJ.9.) ; 

Si (36)= Si (10) 

SUBPROGRAM PLTM AS7“cTrd1^ 




Appendix C S2 ( 1 ) = Cl*C4+Rl 

— Sp ( 2 ) = ( Cl +C5+R2 ) *B 

S2<3>= -(C2*C4+R2)*A 

S? ( 4 ) = ( C 2 * C 5 + R.4J _* A B ! 

S2(5)= C4*C4+R1 

_____ Sp (.6) - ( C4+C5+R2 )*B 

S2(7>=; (C4*C5+R2)*A 

S2<a)= \-(C5*C5-*~R4)*AR • 

S2<9)= 1-S2(2) 

SMlO )=r( Cl*C6+R3 )*n*» ; . 

S2(ll>= ( C2+C5+R3) *AB 

' S 2 (.12 )_=-_( C 2 * C 6 +R 5 ) * A B * R ' 

S2(13)=-S2(6) • 

' Sp ( 14 ) =- ( C4»C6+R3 )*R*b ■ . 

S2(15)=-(C5*C5+R3)*AB 

' S?(16 )= ( C5* C6+R5 ) * AB*B 

S2 (17)= S2(3) " ’ ’ ” " 

S 2 ( 18)_=“S2 (11 L ; : 

S2(19)= (C3*C4+R3)*A*A 

’ S2i 20 ) =- ( C 3*C5+R5 ) * A* AR 

S2(21)=-S2(7) 

S2~( 22 S2-(. 15.) 

S2(23)=-(C4*C6+R3)*A*A 

S ? ( 24 )_=_._( .C.5 *C6+ R5) *A *aB 

; S2(25)=-S2(4) 

_S.2_C.26. >_=__S2 (_1 2.) r T T 

S2(27)=-S2(20) 

L_ S 2 ( 28 ) =-<C 3.*C6±R.6_)_*AB* AB 1 

l S 2 (29)= 52(8) 

i_ S?_(.30 ) =-S2 (16.) : 

i S2(31)=-S2(24) 

| S2(.32)j=„_1C6*C6+E6)„*_AB*AB 

J S2 ( 33 ) = S2(5) 

S 2 ( 34 ) =_ 52 ( 6 ) •' 

' ^(Shu-Sad) 

S2.L36 ) =^S2i 8 ) ' 

S 2 (37)= S2(l) 

1 S2±38)=_S2J2) 

f S2(39)=-S2(3) 

Sp(40)=-S2 (_4j ; J 

■ S2(41)= S2(13) i 

_S 2 .( 4 2).=_S2( 1 4 ) 1 

S2(43)=-S2(15) 

S2-(-44).=-S2(16_) : 

S2<45)= 52<9) 

.S 2 J ..4 6 )j=_S 2 .( 10 .) 

S2(47)=-S2dl) . 

S2.(_48.)=.-S2.(12.) 

S 2 ( 49) z~S2 ( 21 ) 

' _Sp(.5Q).=-S2.( 22.) ! 

S2(51)= S2(23) 

S2l52).=__S2.(24) 

S2(53)=-S2(17) 

S 2154 )=-S 2(18 ) 1 

1 ; S 2 (55)= 52(19) 

59( 56)= S2(?0) ; 1 

S2(57)=-S2(29) 

( S_2_(_5.8JJ==S2_(.3.QJ • 

! S2(59)= S2 (31) 

_ 'SUBPROGRAM P L TMAS: CAR D IMAG E LISTING- 2/4 r — — ! 


S 2 <60)= S2<32> ' Appendix C 

.S?<61)=-S2<25> : : 

S2<62)=-S2(26) 

5 2 . (.63L=_S2(27-> 

S 2 (64)= S2(28> 

S3tl)=— SKI). 

S3<2>=: -SK 2 ) 

S3t3)=i-.S1(3). 

S3<4>= Si (4 ) 

S3<-5 )_=J-S.1(5J r 

S 3 ( 6 ) — \ Si ( 6 ) 

S 3 (.7.) =_rSl (.7 ) : 

S3(8)= 'Sl< 0 ) 

S3(9) = .-SU9) 

S3<10)= Si ( 10 ) 

S3-(l 1 )= T -SU 1 1 ) 

S 3 ( 12) S 1 ( 12 > 

S 3 ( 13) = . SI (13) 

S 3 ( 14 ) =-Sl ( 14 ) 

S 3 ( 1 5 ) = _S 1 ( 1 5 ) : 

S 3 ( 16 ) =-Sl ( 16 ) 

S3ll.7J=_SXU7J 

S 3 ( 18 ) =“Sl ( 18 ) 

S 3 ( 19 ) = _Sl ( 19 ) 

S3<20)=-Sl (20) 

S 3 ( 2 1 ) =— S 1 (. 2 1 .) 

S 3 ( 22 ) = Si ( 22 ) 

S 3 (23.L=-SU23) : ' 

S 3 (24)= Si ( 24 ) 

53. (.25.L=-Sl_( 25) 

S 3 (26) = Si (26) 

S 3(2 7 ) L==S 1.(.2.7J. . 

S 3 ( 28 ) = Si (28) 

S3 i 2 9 ).=rS l.( 2 9 1 

S3(30)= Si ( 30 ) ' 

S3( 3i)=-Sl(3l) 

S 3 (32)= Si (32) 

S3(33)=»Sl (33 ) 

S 3 (34)= Si (34) 

S.3(35)=~Sl (3 5) 

S3 <36)= Si ( 36) 

SC.(.l ) =_“R7 ! 

Sc( 2 )= -R8*B 

ScJ 3 )-=_R 8 *A • 

Sc ( 4 )= -R9*AB 

Sc.(.5)=~-R7- : 

Sc(6)= -R8*B 

Sc(-7-)=-SC.(,3) ; 

SC(8)=-SC(4) 

SC-(9)= -R7 

SC(10)=-SC(2) 

S C-( 1 1 ) R 8 * A ! . . 

Sc ( 12) = R9*AB 

Sc.(l3)=r.R7 : : ! : 

Sc ( 14 ) = R8*B 

SC-(-l 5 ).=.-R8* A , 

Sc(l6)=-R9*AB 

scii-7_i=_co*co : 

Do 25 1 = 1 » 64 ' 

SUBPROGRAM PLTMAS: CARD IMAGE L I ST ING 3/4 1 ’ , 

- - - * * * • 'irti 
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IF<I-17> 5 » 5» 10 

5 _..SC < I )i=_R*SC.a>*A0L 

In ' IF t I"v36) 15»15»20 

_15 Si (1)1= R*S1 ( I)*AB 

S3<I)= R*S3 ( I ) *AB 
.2.0 S? < IJ ^_R*S2UL*AB_ 

25 continue 


-RETURIC 
END 


y 

J - 
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SUBPROGRAM PLTMAS: CARD IMAGE LISTING . 4/4 


Appendix C 


SUBROUTINE ASSYP (A,S1,S2,S3,SC, l,M, J,N, NCP) 

PURPOSE: To assemble (add) the plate element stiffness or mass matrices 

(SI, S2, S3, SC) in the appropriate location in the system free-free 
stiffness or mass matrix (A). 

SUBPROGRAMS 

REQUIRED: LOC 

VARIABLES: A - System free-free stiffness or mass matrices 

SI ,S2,S3, SC - Partitioned element stiffness or mass matrix 
(see Subroutines PLTSTF and PLTMAS) 

I - Index for bay number in x-direction 

M - Number of grid points in x-direction 

J - Index for bay number in y-direction 

N - Number of grid points in y-direction 

NCP - Number of coordinates per grid point 

RESTRICTIONS: None 

ACCURACY: Not Applicable 

SIZE: 000270g 

REFERENCES: None 


123 



Appendix C 



FLOW CHART: SUBROUTINE ASSYP (A, SI , S2, S3, SC, l,M, J,N,NCP) 
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FLOWCHART: SUBROUTINE ASSYP (A,S1,S2,S3,SC, l,M, J,N, NCP) 
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Appendix C 


subroutine. j\ssy_e_< a rSirS2rS3».sc»i^M,j rhJ ,NCpJ_ 

dimension a< 2 ) #si( 2 ) »s 2 < 2 ) »£3<2> »sc<2) 

N C O =N C P*M»N.+.(M-1 ) * (N- i ) 

■ NlM=2*NCP 

_____ NgP=M+ { J-l ) + 1 

NoBAt=(M-l)*(J-l)+I 

ICN=MCP*(N6p_-l.)jtl : 

ICS=ICN+NCP*M 

D_0._20._K = l.rNLM : 

Do 20 L=1»NLM 

; C ALL _L0C J K »L».KLS#_NLM*_1) : 

Call loC(k»l»klg,nLm»o> 

i InK = IC.N±K-1 ■ ; 

InL=ICN+L-1 

; C a L L__L. OCJ I N KjlIML.lN K LlMC.0j_1_) 

! I$K=ICS+K-1 

ISL=ICS+L.-1 

fe call loc(isk»isl»ikl#nco»i> 

fc CaLL_LO.C(INK»_ISL_»_IKS*-NCO#_1.) 

^ IF(K-L) 15» 10 » 10 

i _io A(NkL)_=a_(nklj_+si_(.kls) : 

A(lKL)=A(IKL)+S3(KLs) 

-15 A ( IKS ) =A_UKS )-+S2 (kLG1 

20 Continue 

ICL=.NCP*M*N+N0BA1 : 

Do 70 K=1»17 

•j I FJ.K~4_) 2 5-»-25.>_3.0. 

,25 ICR=ICN+K-1 

__G0__TP_65. 

i 30 lF(K-8) 35r35*40 

_3.5_ I C.R= I C N+NC P+K-5 

Go TO 65 

•__4 Q I F_( K - 1 2 )_4_5 » 45 » 50 

: 45 lCR=lCS+K-9 

: Go__T_Q_.65. : 

50 IF (K-16) 55,55»60 

_55 I C R= I C S+NCP+K- 13 • 

Go TO 65 

6 0 I C R =1 C L 

I 65 Call locucr, icl,icc*nco*i> 

i A( ICC)=A(lCg+SC(K) 

i 7o continue 

! RETURN 

■ 1 L END ' 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 001175 g 
REFERENCES: 


Appendix C 

SUBROUTINE RIB (I, I0P, D, SI , S2, S3) 

To compute the stiffness and consistent mass matrices for a finite 
element representation of a thin-walled open-section beam. 


Zero 

I - Index denoting rib number for constants in common data 
block with PLTVIB (see Figure 3) 

I0P - I0P = 1, stiffness matrix is computed 

I0P = 2, consistent mass matrix is computed 

D - Length of beam element 

S1,S2,S3 - Partitioned stiffness on mass matrix (see Flow Chart) 
D 0 

See Reference 1, p. 8 


Reference 1, pp. 5 - 11 
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FLOW CHART: SUBROUTINE RIB (I, IOP, D, SI , S2, S3) 
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1 Subroutine rib( i »iop.rD#si*.S2.»S3) Appendix C 

dimension er(4) »gr( 4 ) ,rh0(4) »sr 2(4) »c?(4 ) ,sr 3(4) »c3(4) , I 

lAR(4)»A22T4).#A23„(4)-rA33(*4)-»-SvJ{4-L».RE2.(_4J_i:RE.3X4.)..#J5tL(j4j } 

Common er»gr,rho#sR^»c2»sr3»C3» \ 

_1 Ar f A22 » A23»A33»SJ4'RE2.RE3 ,-GM— — j — — 

dimension sk 2) rS2<2) »S3(2> ! 

- ;__CaUlZERO(S 1 ,36) i 

Call I ZERO (S2» 64) 

L C all iZER0(S3>36) 

If(I-3) 5» 10 » 10 

-5-1 A x =1. 0 : : T : 

Go TO \ 15 

— 1 0 A.x = _0 • 0 : 

15 Ay=i • ti|-AX 

■ D2 =d * d: : . : . 

D3=D*P2 

IrUO^-I) 20 > 20 » 25 

20 B22=ER(I)*A22<I)/D3 

_L1 B 2 3=ER(I)*A23d ) V-Q 3 

B33=ER(I)*A33(I)/D3 

R2=-ER(.I.URE2-(-IJ-/D3 * 

R3= ER(I)*RE3(I)/D3 

RX=-SR2(.I-)J* t B22rSR3< I-) *B23-R2- ___ 

j Ry= SR2(I)*B22-SR3<I)*B23+R2 

; Gl= ER( I )_*GM ( I ) /D3+GR il )-*s4( I ) /.( lQ • *D ) _ 

G2= ER(I)*GM(I)/D2+GR(I)*SP(I)/60. 

• G 3 = _£R-a.L*GM.a.)-Z.D+GR CL) *S jJ.I. )..* D /_3 Q_. 

G4 = ER(I)*GM(I)/D-Gr(i)*SJ(I)*D/60» 

1 T.x X= G.1 + SR 3 ( I )_♦ SR 3.( J )-♦ B 3 3 - 2.*-* S R 2_( I„). * SR 3 ( I ) * B 2 3 

; l+SR2(I)*SR2(I)*B22-2.*(SR3(I)*R3-SR2(i)*R2) 

. T 2 *= Q2-( SR 3JI)j* i R 3-SR 2- (IL*R 2) „* D 

Ti Y=Gi+SR3 (I ) ♦SRS ( I ) ♦B33-2 • *SR2 ( I ) *SR 3 ( I ) *B23 
' 1 + S R2 I.I. JJ?JSR2J J JjlB2.2±2 «-*-CSR 2J I ).*R2 ”S^3jLlJL*R 3-L*P / 

T2Y=G2+ ( SR3 ( I ) +R3-SR2 ( I ) *R2 ) *D 

Sld.) = 12.*B22 

j Si<2)= 6**(AY*B22*D-2»*AX*RY) 

Sil3_)_= * ( AY* B 2 2 * D2 + 3 . * A X * T 1 X ) 

Si (4 )= 6** ( 2 • *AY*RX-AX*b22*D) 

5 1(5 )= 6« »( AY *RX * Q±A X*ILY*.0J. . - 

( Si(6)= 4 (3 • *Ay*Tl Y+AX*B22*D2 ) 

S 1(7 )=„ 6»*(AYrjAx)* R 2_*-D ! 

1 Si(8)= 2.*(2.*AY*R2*D2+3»*AX*T2X> 

Si<Q)= — 2 >♦ ( 3« ♦Ay*T?y~2q *AX*R2*_Q2) 1 

Sid0)= 4 • *G3 

S2-<-l-)-= -S 1.(1) 

S2<2)= “Si (2) 

S2 < 3 ) = “SI ( 4 ) - 

S 2 <4>= “SI (7) 

S 2-( 5 )-= —6 • *CAY* B 2 2 *0+ 2 * * A X * R Y.) 

S2 (6) - 2»*(AY*B22*D2“6.*AX*T1X) 

! S 2 ( 7 ) = “6«* ( A X * R.Y L* D - A Y * R X * D.) 

S2<8)= 2»* ( AY*R2*D2-3**AX*T2X) 

S2X9 > = “6 • * ( 2 . *Ay*RX+AX*822*D ) 

Sj> do ) 2 “&•*( ay*rx^d-ax*Ry*d) 

L-I S2XI 1-)-= 2. *.( A X * B22 * D 2-6. * AY*TlY). __ 

S2d2)= 2**(AX*R2*D2+3**AY*T2Y> 

_S|J-13)=-_S1.(.7) j 

S2<14)= 2»*(AY*R2*P2+3»*AX*T2X) 
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Sp ( 15 ) = “2»* (3 .*Ay*T2Y“AX*R2*D2) 

S2 ( 16) =__2»*G4 

S3<1)= Sl(l) 

L. S3 < 2.) = “6** ( A_Y_i 822* D>.2 x*AX„*RYJ 

S3 ( 3,) - Si (3) 

S3 ( 4 ) = 6 •.* ( 2 •. * A Y * R X.+. A X *.B 2 2* D.) 

S3<5)= -SI (5) 

.... S3 <6 )=. Sl.( 61 _ 

S3(7)p “Si (7) 

S3 <-8lf 2 a*_( 2.tAALY*.R2*D2^3j-*AX*X2X ).. 

S3 1.9) s 2»* (3 .*Ay»T2Y-*t2.*Ax*R2*D2) 
S3 ( 1.0.) = 4_. *G3 

return 

25 R M=_ -RHO. ( I )*AR(I)*n/3ftA . 

E2= C2(X)-SR2d) 

1 E3-=— C3.( I.) -SR3 (.IJ : : 

1 R= (E2*E2+E3*E3)/D2 

Rp= ( A22CI )J-A33.a))/.(AR.a).*02) - 

Pl=l3./35. 

; P2-11 • /210 • 

' P3= 1./105. 

P4=_9.7L7Q.. 

P5=l3./420. 

u P6= 1./140. ; 

* S 1 < 1 > = Pl*RM 

Si.(.2 ) .= __(.p.l* A X* E 2±P 2 * A.Y_*DJj*.RM 

Si(3)= ( Pl*AX* ( R+RP) +P3+AY ) *RM*D2 

S;( 4 ).= -IP 2* A X .*. D +Pl *A.Y.* E2.L%B M 

Si (5)= -P2*E2*RM*0 

S i_( 6 > =-< P 3 * A X + P 1 * A Y *.( R+RP.) ).*RM*D2 

i Si(fl)= P2*AX*RM»D3*RP 

S 1 (9 )=-P2 ♦A Y-* R M * D 3 *RP.. : 

Si (10)=P3*RM*D*D3*RP 

SpJ..l.)= P4 + RM 

i Sp(2>= (P4*AX*E2+P5*AY*D)*RM 

S 2 J 3 )-=_j=l(.P.5* A X *Q±P4.*.AY *.E 2 )_* RM 

Sp(5)= (P4*AX*E2-P5*AY*D)*RM 

S ?(6 ) = ( P4* A X* ( R+RP) -P6*AY ) *RM * t P2 

Sp(7)= -P5*( AX-AY )*E2*RM*D 

; S2l8l= P5*A X*R M* D3 * rP ■ 

S2<9)= (P5*AX*D-P4*AY*£2)*RM 

L_ Sp(l0)=.-S2 (7J 

S2 ( 11 ) =- (P6*AX-P4*AY* (R+Rp) ) *RM*D2 

; Sp (.12 )..=__jt R5 * A Y* R M * D.3 *.R P 

. Sp(i4)= -P5*AX*RM*D3 *rP 

S? .(.15 )= P3.*A.Y*.RM.*.D3*RR : 

Sp ( 16) = -P6*RM*D*D3*Rp 

; S3dl= SjJ.l.)_ 

! S3 ( 2 ) = (Pl*AX*E2-Pp*AY*0) *RM 

: l_s 3I3 )= S 1 J.3.) 

! S3(4)= (P2»AX*D-P1*Ay*E2)*RM 

■! S3 .(.5 1=_-tS1.(.5 ) 

S3 ( 6 ) — Si (6) 

S 3 I 8 ) =~r.si.< 8 ) 

S 3 ( 9 ) = -Sl(9) 

S3lLOJ.=_Sllli)i 

_ return 

._ jENQ. . ...... 
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SUBROUTINE ASSYR (A, SI ,S2, S3, l,M, J,N, NCP,AX) 

PURPOSE: To assemble (add) the rib element stiffness or mass matrices (S1,S2, 

S3) in the appropriate location in the system free-free stiffness or 
mass matrix (A). 

SUBPROGRAMS 

REQUIRED: LOC 

VARIABLES: A - System free-free stiffness on mass matrices 

S1,S2,S3 - Partitioned element stiffness or mass matrix (see 
Subroutine Rib) 

I - Index for bay number in x-direction 

M - Number of grid points in x-direction 

J - Index for bay number in y-direction 

N - Number of grid points in y-direction 

NCP - Number of coordinates per grid point 

AX - Logic number: AX = 1 .0 for ribs parallel to x-axis; 

AX = 0.0 for ribs parallel to y-axis 

RESTRICTIONS: M < 4, N 5 4 

ACCURACY: Not Applicable 

SIZE: 000207 g 

REFERENCES: None 
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FLOW CHART: SUBROUTINE ASSYR (A / S1 / S2 / S3,I,M / J / N / NCP / AX) 
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subroutine .assyr ( a »si. s 2 rS3f i^M.J»N»NcP»Ax) 

dimension a ( 2 ) #si<2> »S 2 < 2 > *S3(2> 

A y = l.-AX 

NcO=NCP*M*N+(M-l) *(N- 1 ) 

ngP=m^( j-i.) + 1 ; 

ICN=NGf*(NGp-l)+l 

I CS=I0N+IF.I X( AX) *NCP+aEI XXAYJj* NCP_*M 

Do 20 K=1,NCP 

Do _ 2.0LL= 1 r.N.C P 

Call loc(k»l#kls»ncp,d 

C a L| — L OC ( K_?L»_K.LG.tNCP »oJ - 

lNK=ICN+K-l 

inl=icn+l-i 

Call loc(ink»inl»nkl#nco»i> 

IsKzICS+R-l 

I S L=ICS+L-I 

C A ll _loc„c i sk » I SL»lKL*NCO»lL 

Call loc<ink»isl#iks»nco»i) 

—I F.(K-L) . 15» 10 »10 

10 A(NkL)=A(NKL)+S1(kLS) 

A.(.lKL).=A.aKL-)-±S3-CKLSJ 

15 A(IkS)=A(IKS)+S2(KLG) 

20 continue • 

return 
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PURPOSE: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 000444g 
REFERENCES: 


SUBROUTINE NONDIM (A,M,N,TK,TL) 

To nondimenslonal ize the stiffness and mass matrices of the stiffened 
panel structure. 

A - Stiffness or mass matrix 
M - Number of grid points in the x-direction 
N - Number of grid points in the y-direction 

TK - Nondimensionalizing parameter for force (stiffness matrix) on 
mass (mass matrix). In each case the parameter is calculated 
as the average value of the stiffness or mass of the direct twist 
terms. See Program PLTVIB statement numbers 270 to 290. 

TL - Nondimensionalizing parameter for length. Taken as the 
diagonal length of bay 5 (see Figure 3) of the structure. 

TK / 0, TL / 0 

Not Applicable 


None 
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FLOW CHART: SUBROUTINE NONDIM (A,M,N,TK,TL) 
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FLOWCHART: SUBROUTINE NONDIM'(A,M, N,TK,TL) 
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FLOW CHART: SUBROUTINE NONDIM (A,M, N,TK,TL) 
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Appendix C ’ SuBROu.T I NE_ NONDIM ( A r^N » J-K jlT-L ) 

dimension A(2> 

; nc.p=_4 ' 

NCl=NCP*M*N 

NC0=NC I + 

N|_M=NCP*(M*N-M-2)+l j . 

Nb A ISID =N CP * (M +l ) ■ 

MbAY=M-1 

N rAY=N-1 

Do 100 IC=NCP»NCI,NCP 

II = I Ct-nCP+1 

; Do 35;I=1,NCP 

• IrPw=I1 + I-1 ___ 

Do 35 |J=I»NCP 

IcQL=tll+g-l 

ij=i+|g*g-g)/2 

C aLl li0CXlR0W* LC_0L.aR-.-NiC0_Kl) 

Go TO (5,l0,l5»10,l5,15»l5»20»20»25) »iJ 

_5 Cl=IL*TL/_TK 

: Go TO 30 

i_in Cl.=tl/.ik 

i Go TO 30 

_L5 C L=lvO 7TK 

1 Go TO 30 

2 0 _C L= 1.-0 / ML*TK ) 

Go TO 30 

_2 5 Cl= 1 . 0 /.. (T L*J L.*XK_) __ 

30 A<Ir)=CL*A<IR) 

_35 Continue 

IF(IC-NCI) 40.100.100 

L.4.0 i c 1= I.C.+.1 • : 

If<IC1-NLm) 50.50,45 

_ 33 . NrAND=NBAn D~NCP-H - 

50 Nu=IC+NBAND 

D.o_ .9 5_J.C.= l_c_l .NU f.N.cP . 

If(JC-NCI) 55»95»95 

55 D o 90 J=1,NC P 

IC0L=JC+J-1 

• PO 90 1=1, NCP 

IrOw=I1+I-1 

; c aLl_loc x.i R.o w-f-1-COL.ajR ,.n.C-0.,-1j__ 

! Ij=I + NCP*(nJ-1) 

Go _T 0_J 60,65. 6 5 , 7Q . $5 , ?Q ,J7_0,_75 ,65.70, 70» 75 » 

l7o.75.75.80) ,10 

_60 Cl=TL*TL/TK ■ • 

Go TO 85 

_65 C L=TL/IK 

Go TO 85 

_7.0 Cl=1,0/IK : 

Go TO 85 

-7-5 CL= 1-.-Q7.LT L*IKJ 

Go TO 85 

_8 0 C l= 1.. 0/ (-lL*XL*TK-) ' 

85 A(IR)=CL*A<IR> 

_90 CoNXlNUE 

'95 Continue 


Do 155 u=i»nbay 
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Do 155 I= 1 »MBAY 

IrAY=I+MBAY*(J-1) . 

— n 

•i 

IcOL=NCI+I 8AY 
IrS=(IC0L*lC0L-ICoL> /o 

IC=NCP*IBAY+NCP* ( j-i ) “ 
’ Dn 1 50 L=1 j 4 

1 1Q5 

Go TO (105»ll0»ll5»120) »L 
IrC=lC 


\ J.Q 

Go TO 125 
IrC=lC+NCP 

. 

—115— 

U-12.0— 

Go TO 125 
ICC=IC+NCP*M 


Go TO 125 
IrC=lC+NCP*(M+l) 


1 125 

i 

Do 150 K= 1 »NCP 
IrR=ICC-NCP+K 


IrC=Irr+ICS 

Go TO ( 130 » 135 » 1 3«S» 1 4n ) 

130 

Cl=TL*TL/TK 
Go TO 145 

• 

; 135 

Cl=TL/TK 
Go TO 1 45 


140 

_1.4.5_ 

C|_=1./TK 
A ( IRC ) =CL*A ( IRC ) 


150 

_155— 

Continue 

Continue 


l 

Cl=tl*tl/tk 

N«;=nCI+1 


DO 160 I=NSiNCO 
Il=I*(I+l)/2 • 

-160— 

A ( 1 1 ) =CL*A (II) 
Continue ■ 


return 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 

RESTRICTIONS: 

ACCURACY: 
SIZE: 00002 1 
REFERENCES: 


PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 

ACCURACY: 

SIZE: 0000 76 0 
8 

REFERENCES: 


SUBROUTINE ZERO (A, N) 

To set the first N elements of a single subscripted array to zero. 
None 

A - Array to be initialized 
N - Upper limit for elements of A set to zero 

N must be equal to or less than the dimensioned size of A in the 
calling program. 

Not Applicable 


None 

SUBROUTINE FILL (A, B, N, MS) 

To fill array A with elements of array B 

LOC 

A - Single subscripted array 

B - Single subscripted array (symmetric storage mode: LOC) 

N - Number of rows or columns in A or B 

MS - Logic number : MS 0 subroutine is bypassed; MS = 0, A is 
a general matrix (NXN) and B is a symmetric matrix; MS = 1, 
A and B are both symmetric matrices 


None 

Not Applicable 


See Subroutine LOC 
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FLOWCHART: SUBROUTINE FILL (A,B,N,MS) 
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-SUBROUT INE.ZEROXAj_NJL 
Dimension A(2) 

_D0-_5_I = Lf-N 

A{ I) =0.0 

-Continue 

return 

end : ■ 


ii SUBPROGRAM ZERO: CARD IMAGE LISTING; 


SUBROUTINE- FI LL.IA , B»N»MSi- 

DiMENSION A ( 2 ) »B(2) 

IF_(.MS)_25.*.5j_15 ! 

5 Do 10 I=1*N 

_____ D o -1 0- U = 1 » N_ _____ 

call loC(i»j»ijs»n»i) 

C A L| LOCJ IJGaNilO) 

A(IJG)=B(US) 

_io continue 

return 

_1 5 — N U P= N *XN±UY-2_ 

Do 20 1=1, NUP 

A ( IJ =BjLU - 

20 continue 

25 _8ElU-RN 


'SUBPROGRAM FILL: CARD IMAGE LISTING 
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PURPOSE: 

SUBPROGRAMS 

REQUIRED: 

VARIABLES: 


RESTRICTIONS: 
ACCURACY: 
SIZE: 000333 g 
REFERENCES: 


SUBROUTINE ORDER (A, NDEL, NCP,MI, Nl, NDL) 

To remove (set to zero) specified (constrained) coordinates in array 
A, reorder array A, and calculate the new size of array A. 

DELETE, LOC 

A(l) - Stiffness or mass matrix of structural idealization 

NDEL(I) - Array of logic numbers: see PLTVIB 

NCP - Number of coordinates at each grid point 

Ml - Number of grid points in x-direction 

NJ - Number of grid points in y-direction 

NDL - Number of coordinates removed by this subprogram 

NC(I) - Array of coordinate numbers for which NDEL(I) = 1 

NC(I) must be dimensioned the same as NDEL(I) in program PLTVIB 

Not Applicable 


None 
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FLOWCHART: SUBROUTINE ORDER (A, NDEL, NCP,MI, Nl, NDL) 
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FLOWCHART: SUBROUTINE ORDER (A, NDEL, NCP, Ml, Nl, NDL) 
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FLOW CHART: SUBROUTINE ORDER (A, NDEL, NCP,MI, NI,NDL) 
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FLOWCHART: SUBROUTINE ORDER (A, NDEL, NCP,MI, Nl, NDL) 
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SuB R OUJ.I N£_ORQj^J. AiJ^DELf NjC P jlMJjlN IjtliDLji 

' DIMENSION A(2) »NDEi-(2) ,NC(73) 

i NdL=0 

NgP=MI*NI 

: NCO=NCP*Nr,P-KMI«l)*tN l»l) 

i Do X I=l»NCO 

: Nc < I>). =0 

!~I continue 

; D.O_3Ll.C=l.rNG£ 1 

I=NCP*(IC-1) 

DO. 3 k=.l » N.C.P - - 

Il = l*iJ 

• I R i NDEL11 1 L)_3 jl3_l2 

2 NpL=NDL+l 

N C( NDL.)_=ii : 

; 3 continue 

I R_.( NOL)_loaOf-3l 

31 Do 9 K=1»NDL 

MzNCI K.) H 

Ml=NCO-K 

■ N-MLj-,1 

! Call delete(a»n#M) 

I OO. 7.-.I.=M.»-ML : 

1 Il=I+l 

■ Po._7_*L; i.»N 

Call loc<i»j,ij»n,1) 

: ,I.EJ.L=J.)__5-»JL»J5 

4 Jl=j+1 

_C. A L I L O.C..(Xl 

Go TO 6 * 

_5 CALL_LOC_(ILr-JjJDr»Ji#a-) 

6 A(IJ)=A(IJ1) 

7 C O.NX.INUE. 

' DO 9 I=l»NOL 

J Fj NC(I )-M) 9»9»6 

8 Nc(I)=NC(I)-l 

: 9 Con tin ue 

M-NCO-NDL+l 

D o 1 0 K=M, NCQ , ~ 

Call delete<a»nco,k> 

_xo .continue : . 

; RETURN 

I f END L - - — - 



i SuBROuxI NE_ DELETE 

I dimension A(2> 

1 : DQ..-1-KjOjlN 

Call LOC(K»J»KJ»N,l) 

L_ A ( K J) =0 *0 

j i Continue 

reTurn_ 

■ • . END] 
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APPENDIX D 


COMPUTER PROGRAM USER'S MANUAL 
INTRODUCTION 

This appendix defines the input data parameters, input data format and output 
data format for the one-dimensional panel computer program and the two-dimensional 
panel computer program described in Appendix C. An example problem is included for 
each computer program listing typical input and output data formats. 

ONE-DIMENSIONAL PANEL ARRAY 

The one-dimensional panel geometry is illustrated in figure D-l. The structural 
idealization considers a fundamental mode across the width of the structure (y-direction) 
so that the finite element model considers only parameter variations along the length of 
the structure (strip analogy). The finite element model of the one-dimensional panel 
array then reduces to a simulation of a spring -supported beam. The elastic supports are 
modeled using the finite element representation of a thin-walled open-section beam as 
described in reference 1. The lumped parameter model of the elastic supports is given 
by equations 20a and 20b of reference 1 . The definition of the geometric constants 
given in equations 20a and 20b of reference 1 are defined by figure D-2. The stiffener 
geometry for a stiffener parallel to the x-axis is presented in figure D-3 for completeness. 
The definition of the required input data is given below. The choice of structural ideal- 
ization is best illustrated by an example problem. 

ONE-DIMENSIONAL PANEL ARRAY: 

DEFINITION OF INPUT VARIABLES 

The modal analysis of one-dimensional panel arrays is preformed using program 
BMPROP and the associated subprograms (see page 63). Definition of the required in- 
put variables is as follows: 

NCASE Number of data cases to be processed 

NDATA A four digit data case identification number 

NBAY The number of panel bays of the structure (not greater than 5) 
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NSUP 

IOUT 


IBL 


IBR 

NINT 

BW 

PR 

NEL(I) 

EI(D 

WB(I) 

BL(I) 

NCP(I) 

SL(I) 

SC(I) 

SR(D 

RL(l) 


The number of elastic supports 

Data output option: IOUT = 0; print parameters NDATA, TK, TL, TM, 
frequencies and modal amplitudes at each grid point; IOUT > 0; print 
output for IOUT = 0 and the normalized modal displacement, shear, and 
bending moment distribution. 

Logic number for applying clamped constraints at the left-hand end of the 
structure: IBL = 0; left-hand end is free or elastically supported: IBL = 1 
the left-hand end of the structure is clamped (translation and slopes set to 
zero) 

Logic number for applying clamped constraints at the right-hand end of the 
structure (definition identical to IBL) 

Number of points interior to each element for which normalized displace- 
ment shear and bending moment distributions are to be calculated (equal to 
or less than 5) 

Width b of the panel, inches, (see figure D-1) 

Poisson's ratio for the cover sheet material 

th 

Number of elements in the I panel bay 

Bending rigidity of 1^ panel bay (see equation 21, reference 1, p. 11), 
lb. -in . 

th 

Weight per unit length of I panel bay (see equation 22, reference 1, 
p. 11), lb. /in. 
th 

Length of I panel bay (dimension a., figure D-1) 

Translation coordinate number for locating I elastic support in structure 
NCP(I) will be an odd number 

th 

Lumped spring constant, K , for the I elastic support (see equation 20a, 
reference 1 , p. 1 1 ) 

th 

Lumped spring constant, K z q, for the I elastic support (see equation 20a, 
reference 1 , p. 1 1 ) 

th 

Lumped spring constant, Kgg, for the I elastic support (see equation 20a, 
reference 1, p. 11) 

th 

Lumped mass constant, I* , for the I elastic support (see equation 20b, 
reference 1, p. 11) 
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th 

RC(I) Lumped mass constant, l*g, for the I elastic support (see equation 20b, 

reference 1, p. 11) 

th 

RR(I) Lumped mass constant, Igg, for the I elastic support (see equation 20b, 

reference 1, p. 11) 

Note: To add a lumped mass at a coordinate, input an elastic support with zero 

stiffness at that coordinate (i .e . , SL(I) = SC(I) = SR(I) = 0 at NCP(I)) 

The input data format is as follows: 

CARD 0 (ONE CARD PER DATA SET = NCASE DATA CASES) 


COL (FORMAT) 

1(15) 

NAME 

NCASE 


DATA CASE INPUT FORMAT 

CARD 1 (ONE CARD PER DATA CASE) 


COL (FORMAT) 

1(15) 

6(13) 

9(13) 

12(13) 

16(13) 

18(13) 

21(13) 

NAME 

NDATA 

NBAY 

NSUP 

IOUT 

IBL 

IBR 

NINT 


CARD 2 (ONE CARD PER DATA CASE) 


COL (FORMAT 

4(E 1 2.5) 

1 6(E 12.5) 


JCTURE BAY) 

NAME 

BW 

PR 


CARDS 3 through 3+ NBAY (ONE CARD FOR EACH STRL 

COL (FORMAT) 

1(13) 

4(E 12.5) 

1 6(E 12.5) 

28(E12 .5) 

NAME 

NEL(I) 

EI(D 

WB(I) 

BL(I) 


CARDS 4 + NBAY through 4 + NBAY + 2 NSUP (TWO CARDS FOR EACH ELASTIC 
SUPPORT) 


COL (FORMAT) 

1(13) 

5(E12.5) 

16(E12.5) 

28(E12.5) 

NAME 

NCP(I) 

SL(I) 

SC(I) 

SR(I) 


COL (FORMAT) 

5(E12.5) 

16(E12.5) 

28(E 1 2 . 5) 

NAME 

RL(I) 

RC(I) 

RR(I) 
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ONE-DIMENSIONAL PANEL ARRAY: EXAMPLE 

To illustrate a typical structural idealization consider a one-dimensional panel 
array consisting of three bays with two interior elastic supports and both ends clamped. 
Suppose that the span of each bay is 6.0 inches and that the width of the structure is 
20.0 inches. The cover sheet is 0.032 in. and the material is 7075-T6 aluminum alloy. 
The stiffener is taken as a zee section as illustrated in figure 13, reference 1, and the 
stiffener orientation is taken such that both stiffeners face in the same direction. This 
data case corresponds to specimen SPI-2-1 for the three-bay configuration described in 
figure 15 of reference 1. Assuming that the stiffener attach point is directly below the 
vertical web and that the mass of the skin directly adjacent to the stiffener flange is 
considered to act with the translational inertia term, I* , for the lumped mass represen- 
tation the stiffener data is 


= 630.02 

1* = 0.074299 

zz 

zz 

z9 - ,4 ' 271 

N * 
CD 

II 

O 

o 

99 ' ,18 - 30 

l* Q = 0.0074605 


Since the cover sheet is uniform, the lumped data are 

E 1(1) = 313.34 WB(I) = 0.03232 I = 1,2,3. 

Three* elements are used to model each bay. 

The above structural idealization is illustrated in figure D-4. The input data format is 
illustrated in figure D-5. The output data format is illustrated in figure D-6. In fig- 
ure D-6 the edited input data is printed, the eigenvalues and eigenvectors are printed, 
and the interpolated values of the displacement (W), slope (DW/DX), shear (V), and 
bending moment (M) are printed for the fifteen lower frequency modes (only the funda- 
mental mode is illustrated). The sequence for listing the eigenvector is dimensionless 
displacement and rotation for each grid point across the structure (see figure D-4). The 
experimental values for frequency, mode shape, and bending moment distribution 
(strain) are given in Table II and figure 26 of reference 1 . 

*Experience has shown that (NBAY-1) elements for each bay insures satisfactory 
frequency convergence of the lower NBAY modes. 
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NINE-BAY TWO-DIMENSIONAL PANEL ARRAY 


The two-dimensional panel array considered here is the nine-bay configuration 
illustrated in figure D-7. The structural idealization considers nine plate elements 
forming the cover sheet and orthogonal stiffeners. The plate elements are described in 
reference 1 and consider an interior fundamental clamped-clamped plate mode as a gen- 
eralized coordinate. The stiffener elements are taken as thin-walled open-section 
beams as described in reference 1. The stiffener geometry is illustrated in figures D-2 
and D-3. The element nomenclature is illustrated in figure 3. Definition of the input 
data is given below. 

NINE-BAY TWO-DIMENSIONAL PANEL ARRAY: 

DEFINITION OF INPUT VARIABLES 

The modal analysis of nine-bay two-dimensional panel arrays is performed using 
program PLTVIB and the associated subprograms (see page 103). Definition of the re- 
quired input variables (see Subroutine RDNWRT (1)) is as follows: 

NDATA Number of data cases to be processed 

NCASE A four digit data case identification number 

th 

A(|) Length of the I rib segment parallel to the x-axis (see figure D-7), inches 

B(l) Length of the l^ 1 rib segment parallel to the y-axis (see figure D-7), inches 

2 

EP Young's modulus of the cover sheet material, Ibf./in. 

HP Thickness of the cover sheet material, inches 


PR Poisson's ratio of the cover sheet material 

o 

RHOP Weight density of the cover sheet material, Ibf./in. 
ER(I) Young's modulus for the l^ 1 rib, Ibf./in? 

GR(I) Shear modulus for the l^ 1 rib, Ibf./in? 

RHO(l) Weight density for the 1^ rib, Ibf./in? 


SR2 (I ) S y (I = 1,2) and S x (I = 3,4), inches* 

C2(l) C (I = 1,2) and C (I = 3,4), inches* 

y x 

SR3(I) S z , inches* 

C3(l) C z , inches* 

*See figures D-2 and D-3 
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th 2 

AR(I) Cross sectional area of I rib, in. 

A22(l) I (I =1,2) and I (I = 3,4) of the l^ 1 rib, inf (equation 5a, ref. 1) 

yy xx 

A23(l) I (I = 1,2) and l^ z (I = 3,4) of the I* rib, inf (equation 5a, ref. 1) 

A33(l) l zz of the l^ 1 rib, inf (equation 5a, ref. 1) 

SJ(I) St. Venant's torsion constant, inf (equation 5a, ref. 1) 

RE2(I) R g y (I = 1,2) and R gx (I = 3,4), inf** (equation 5b, ref. 1) 

RE3(I) R> ez for the 1^ rib, inf** (equation 5b, ref. 1) 

GM(I) Warping constant for the l^" 1 rib, inf 

The input data format is as follows: 

CARD 0 (ONE CARD PER DATA SET) 


COL (FORMAT) 

6(13) 

NAME 

NDATA 


CARD 1 (ONE CARD PER DATA CASE) 


COL (FORMAT) 

6(14) 

NAME 

NCASE 


CARD 2 (ONE CARD PER DATA CASE) 


COL (FORMAT) 

5(E12.5) 

17(E12.5) 

29(E12.5) 


NAME 

A(l) 

A(2) 

A(3) 


CARD 3 (ONE CARD PER DATA CASE) 


COL (FORMAT) 

5(E12.5) 

1 7(E 1 2 .5) 

29(E 1 2.5) 


NAME 

B(1 ) 

B(2) 

B(3) 


CARD 4 (ONE CARD PER DATA CASE) 


COL (FORMAT) 

5(E12.5) 

17(E12.5) 

29(E12.5) 


NAME 

EP 

HP 

PR 

RHOP 

CARDS 5, 9, 13, 17 (ONE EACH PER DATA CASE) 

COL (FORMAT) 

5(E12.5) 

T7(E1 2.5) 

29(E 1 2.5) 


NAME 

ER(D 


mssm 



**See comment on page 6 of reference 1 . 


154 

























Appendix D 


CARDS 6, 10, 14, 18 (ONE EACH PER DATA CASE) 


COL (FORMAT) 

5(E12.5) 

17(E12.5) 

29(E1 2.5) 

41 (El 2 .5) 

NAME 

SR2(I) 

C2(l) 

SR3(I) 

C3(l) 


CARDS 7, 11, 15, 19 (ONE EACH PER DATA CASE) 


COL (FORMAT) 

5(E 1 2 . 5) 

17(E12 .5) 

29(E12.5) 

41(E12.5) 

NAME 

AR(I) 

A22(l) 

A23(l) 

A33(l) 


CARDS 8, 12, 16, 20 (ONE EACH PER DATA CASE) 


COL (FORMAT) 

5(E12.5) 

17(E12 .5) 

29(E12.5) 

41 (E12.5) 

NAME 

SJ(I) 

RE2(I) 

RE3(I) 

GM(I) 


NINE-BAY TWO-DIMENSIONAL PANEL ARRAY: EXAMPLE 


As an example consider a square nine-bay stiffened panel with identical stiffeners 
uniformly spaced in x- and y-directions . This example is selected to illustrate the na- 
ture of repeated roots and symmetry in the eigenvector (mode shape). Assume that each 
bay has dimensions (A(l) x B(l)) of 10.0 in. x 10.0 in., that the aluminum cover sheet 
is 0.032 in. thick, and that the stiffeners are described by the following data: 


ER(I) = 10.3 x 10 6 , Ibf./in? 
GR(I) = 3.9 x 10 6 , Ibf./in? 
RHO(I) = 0.101, Ibf./in? 
SR2(I) = 0.0, in. 

C2(l) =0.0, in. 

SR3(I) = 0.27191, in. 

C3(l) = 0.80525, in. 

AR(I) =0.120, in? 


A22(l)= 0.018389, inf 
A23(l)= 0.0, inf 
A33(l) = 0.022298, inf 
SJ(I) = 4.08 x 10“ 5 , inf 
RE2(I) = 0.0, in? 

RE3(I) = 0.0, in? 

GM(I)= 2.4445 x 10" 3 , in? 
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The stiffener cross-section shape and its attachment to the cover sheet is illustrated 
in figure D-8. Figure D-9 illustrates the input data format. Figure D-10 illustrates the 
output data format where the edited input data is printed and the eigenvalues and eigen- 
vectors are printed. The sequence for printing the eigenvectors is illustrated in figure 
D-l 1 . The mode shapes for the four lower frequency modes are given in figure D-12. 



FIGURE D-l . ONE-DIMENSIONAL PANEL ARRAY 
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FIGURE D-2. CROSS-SECTION GEOMETRY FOR RIBS PARALLEL TO Y-AXIS 



FIGURE D-3. CROSS-SECTION GEOMETRY FOR RIBS PARALLEL TO X-AXIS 
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a) Section through Centerline of Structure 



b) Coordinates for Initial (unconstrained) Structural Idealization 



I 1 3 5 7 9 11 13 15 


z 

c) Coordinates for Final (constrained) Structural Idealization 


FIGURE D-4. STRUCTURAL IDEALIZATION: EXAMPLE PROBLEM THREE BAY 
ONE-DIMENSIONAL PANEL 
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FREE VIBRATION of a ONE DIMENSIONAL PANEL ARRAY 

\ data 

CASE 12 

number of bays= 3 

number oe SUPPOrTS= 2 

panel WIDTHS 0.20000E+02 

POISSON’S RATl0= 0.32050E+00 


BAY NUMBER OF BENDING ' WEIGHT PER rAY 

NUMBER ELEMENTS RTGtDITy UNIT ApEA LfNGTH 


l 


1 

13 

e 

0«3l334E+03 

0.32320E-01 0.60000E+01 

2 

3 

0 • 31334E+03 

0.32320E-01 0.60000E+01 

• 3 

3 

0.31334E+03 

0.32320E-01 0.60QO0E+01 

SUPPORT NO. 1 



Input coordinates 7 

KZZ- 

0.63002E+03 

KZTH£TA= 0.14271E+03 

kTHEtA= 0.11fl30E+03 

IZZ- 

0 . 74299E-01 

iztheta= o.oooooe+oo 

ITHEtAs 0 . 74(Tj05E“02 

SUPPORT NO. 2 



Input coordinates i3 

. KZZ- 

0.63002E+03 

KZTHETA= 0.14271E+03 

KTHETAr 0.11b30E+03 

IZZ- 

0.74299E-01 

iztheta= o.oooooe+oo 

ITHEtAs 0.74605E-02 


FIGURE D-6. OUTPUT DA TA FORMAT: EDITED OUTPUT PAT/C EXAMPL E PRO BLEM 


free vibration of a one-dimensional panel’ array 
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FIGURE D-6. OUTPUT DATA FORMAT: NORMALIZED MODE, 
SHAPE AND STRESS RESULTANTS 
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FIGURE D-7. NINE BAY TWO-DIMENSIONAL STRUCTURE 
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FIGURE D-8. STIFFENER CROSS-SECTION SHAPE: EXAMPLE PROBLEM 
NINE-BAY TWO-DIMENSIONAL PANEL 
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Appendix D 


DATA CASE 1000 


free vibration of a 
NINE'BAy ortho g onally stiffenfd panel 


STRUCTURAL GEOMETRY 


a 

A l= Q .lQofo OE+ n ? A 2 = 0.10000E4.Q2 A 3 = 0 . ln Q OQE +02 
Bl= 0 « 10doQE+n2 R2= 0 » 1QQ00E+Q2 B3= 0 . 1 nOOQE+02 


COvEr sheet data 


• YOUnG*S MODuLuS= 0.l0.30oE+0fi POISSON’S RATlO= 0.32050E+00 

thickness - o,32 n onF-oi w Eight/VoLume = o.inioo F+oo 

BENDING RlGlDITYr 0.31346^+02 


■. stif fen er data 

sti ff e ner s parallel to x-Ax is 

STIFFENER no. 1 E = 0.10300E+08 G = 0.39000E+07 RHO = O.lolOOE+OO 

_SY= _0 . 0QQ00E+00 CY = O .OOOoOE+O O SZ = 0.27 l9l E+00 CZ = 0 ^80525E+00_J 

A = 0.12000E+00 IYy= 0.183a9E-6l IYZ= 0.00000'E+00 IZZ = O.22290E-Q1 ' 

__J = 0.40800E-04_ REY= 0 . QOOOOE+ QO REZ= O.OO OQOE +OO G A MmA= 0.2444 5E-02 


stiffeners paral l el tQ x -axis 


STIFFENER no. 2 

E = 

0.10300E+08 

G = 0.39000E+07 

RHO = 

0.10100E+00 

SY= O.OOOOOE+OO 

CY = 

O.OOOOOE+OO 

sz = 0.27I91E+00 

CZ = 

0.80525E+00 

A = 0 . 12000E+00 

iy y = 

0.163fl9E“0l 

IYZ= O.OOOOOE+OO 

IZZ = 

0.22298E-01 

J = 0.40800E-04 

REy= 

O.OOOOOE+OO 

REZ= O.OOOOOE+OO 

GAMmA= 

0.24445E-02 

STIFFENERS PARALLEL TO Y-Axis 

STIFFENER no. 3 

E = 

0.10300E+08 

G = 0.3g000E+07 

RHO = 

0.10100E+00 

SX= 0 , OOOOOE+OO 

cx = 

O.OOOOOE+OO' 

sz = 0.27191E+00 

CZ = 

0.80525E+00 

A = 0.12000E+00 

IXx= 

0.183fi9E-01 

IXZ= O.OOOOOE+OO 

IZZ = 

0.22298E-01 

J = 0.4060CE-04 

REx= 

O.COOOOE+OO 

REZ= O.OOOOOE+OO 

GAMmA= 

0.24445E-02 

stiffeners parallel to y-axis 

STIFFENER no. 4 

E = 

0.10300E+08 

G = 0.39000E+07 

RHO = 

0.10100E+00 

SX= O.OOOOOE+OO 

CX = 

O.OOOOOE+OO 

SZ = 0.27191E+00 

CZ = 

0.80525E+00 

A = 0 . 12000E+00 

IXx= 

0.18389E-01 

IXZ= O.OOOOOE+OO 

IZZ = 

0.22298E-01j; 

0.24445E-02V 

J = 0.40800E-04 

REx= 

O.OOOOOE+OO 

REZ= O.OOOOOE+OO 

GAMmA= 

L- ■ ■ ■ 
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Appendix D 


DATA case 1000 


\ nondimensionalizing constants 


. EI.Ge N V A L UE =_ _0_. 3 2 76 3 E+OX- 


_-0.».l 0 599E+00- 


eigenvector 

0_. JLQ.59 9E±00_ 


_-0-.6o316E-02 

0.10599E+00 0«'10599E+00 -0.47696E+00 

.-0...1 n 599E+-0 0 -0* 

-0.47696E+00 -0»24869E-02 0.62428E-08 

_0 • 37093E-09 Q..46L36F-^-09 Q> 

Si 

EIGENVALUES Q.24fe3F+n i 


■ m 


0.4 1)6 68E-0.2_ 
0.27893E+00 
_0 . 7791 IE-01. 

0.23973E+00 
-0 .41.871E-07 


/ \ eigenvector 

_o A 

0 • 779o7E-(1 1 0.50451E + 00 

0 *5 0.45 2 E*43 0 ^0__4_Q668E-Q2 

0.35S53E-02 0.98020E-02 

0-._3487_3.E-Q2 


_ErE Q UE N C_Yr_ 0-»-2 5 8 2 ft E + 0 4 r H Z . 


^a*4-7-696E+j)0 _o. .6 0 3 16E---0 2. 

0.60316E-02 -0.1Q599E+00 


0.24869E-02 


0.12935E-08 

:-oa. 


_EREaUENCy.-^o.22ll6E±OA»_HZ. 


-0.85587E-02 

_CU29630E+aI7 


0.27892E+00 

0.12675E+00. 


0.74823E-Q2 -0.34873E-02 

llQ_. 24.8 2.3E.7 0 2 rO*98ai9Ero2 rQ.35553E-.02_ 


.EIGENVALUES Q-.2~4.Q-23£±D_1. 


•0 . 85587E-02 -0 ._7_7.9l.lEnU- 

■0.12674E+00 . 0.29630E+00 

_0_. 2 9 63 .0 E + 0 Q -0 » 23 9.7-4£±a0_ 
•0.50451E+00 -0«74823E-02 

L0-.-5385OEr.Q2 0^980 20 EniZ. 


EIGENVECTOR 

0_.278.92E+00_ 

0.23972E+00 
...0*85 587 E “0.2 . 
-0.34873E-02 


JEI GfNVALUE= 0.2 31 87F+01 


_Q_. G-3.4 0 8E-Q 7 - Q.35356 E+fl n 


EIGENVECTOR 

-Q.353 56E+00 


0.35355E+00 -0.3535&E+00 -0.28779E-05 

_ 0_._3 5.3 5 5 E + 00 -0 .2 66^ 8F-na -Q.2 7 Q49E-Q7 

-0.5Q723E-05 0.4’394lE-07 -0 . 12487E~06 

.-0,4j365E-07 -0 . 19978E-0G 0 .55218 E-Q7 


EIGENVA LUES 0.89302E+00 


Q_._i_7.Z8 4 E - 0l_ 


EIGENVECTOR 
Q.35 l52 F+nn -0 . 35152E+00 


0,5Q452E+o0 Ol._40 668 E - 0 2_ 


-0.40668E-02 

-O.779O8E-0I. 


0.35553E-02 

Q^-3487ijLE-_g2_ 


-0.12674E+00 
.. 0.27893E+00._ 
“0.98019E-02 
0 »74 823£r-02_ 


FR EQUENC Yr Q ,2 l 72flE+Q4 . HZ . 


~Q.4B897 E~n5 -Q. lflQllF-06 
0.15584E-06 -0.35356E+00 

0.35355E+n0 0 .35 355F-t-00 

-O.I38I8E-06 0.24371E-06 

0 . 1643 3E-Q6 _ -Q . 1Q497E-07 


0.35152E+00 0.35152E+00 -0.38818E-01 

L 0 . 35152E+00 -0 . 388 1 9E -nl -0 .17764E-Q1 

0.3r818E-01 0.19067E-01 0.68695E“02 


FrEquENCY= 0 « 13484 E+ 04 1 HZ , 

0.388l8E~nl -0.1776 4E -01 
-0.17764E-01 -0.35152E+00 
-0. 3515 2E+0 0 0..-35 15 2 E + 0 .0 


0.19067E-Q1 0.68695E-02 


-0.5Q2 32E - 01 Q.66695E-02 0. 19067E-Q1 0.68695E-n2 0.19 067 E-01 


E IGfNvAL U Es 0 «fil502F+0Q 


frequency- o.i2882E+o4rHZ . 


0.-2 l3.90.Er01 r0-«.2 12.7.A£+(jO_ 

0.11434E+00 -0.44260E+O0 

.-0...44260E+.0 0 Q .550n lE-nl 

-0.11584E+00 -0.19860E-01 

-0.1nl39E-Q7 -Q.460A7E-n2.-_ 


ElGENVECTOR 

_Q_..47258E+0_Q_ 


_0-.ll565E+_oO. 


•0.55088E-01 -0.10171E-01 

l 0_.2.13_90E-0.1 -0.2l274 E+n0 


.0_t.l0.l7 1 £r0 1_ 
0.H434E+00 
„0..47758E±0.0_ 


0.12966E-01 -0.94438E-02 0.46087E-02 

_O...9_4_4.37Er_0.2 Q_dl9_8_59£r.ai_ 


r 
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EIGENVALUES 

0.^1502E+00 


FrEquENCY = n 

* 12882^+04 »HZ, 

; -0 ♦ 1q171E-01„ 

0..44260E+00 

eigenvector 

0.11433E+00 

- 0.550 78E-Q1 

0.21390E-01 

; -0.47758E+00 
: -0 .2i 274E+00 

-0.21274E+00 
0...1 15fl6E4fiQ_ 

-0.11586E+00 

0.10171E-01 

-0.21390E-61 

0.44260E+hO 

-0.47758E+00 
0 • 1 1 434F+00 

0.55060E-01 
i 0.35098E-07 

0 .94437E-02 
0«';12966E-01 

0.46087E-02 

0.19860E-01 

-0.l9859E-ol 

-0.46088E-02 

-0.12966E-01 

-0.94437E-02 

1 EIGENVALUE: 

0.75190E+00 

EIGENVECTOR 

FrEqUENCY= 0 

• 12373E+04 »HZ . 

1 

' 0.11632E-07 

■_L 0.35283E + 00. 
\ 0.35283E+00 
0.66447E-07 

0 .3|52a3E+oO 

r O « 3-5 2 3E + G.0_ 

0«ll360E-06 

- 0 .sii 2 PE.-Ga_ 

0.35283E+00 
0_. 44.4 39£_“06_ 
-0.68469E-08 
0.31864E-01 

0.14321E-06 

-0.35283E+00 

-0.74018E-O8 

0.42211E-08 ' 
_r0..._35283E+00_ 
“0.35283E+00 
“0 . 3l 864F— 01 

.5 0.45107E-08 

-0.31864E-01 

0.39202E-08 

0.3i864E-ol 

0.21881E-08 

j EIGENVALUES 

0.71075E+00 

EIGENVECTOR 

FrEqUENCYz 0 

. l2030E+04rHZ. 

•-0.14213E-04 
0 • 1g918E-02 

-0 • lb9i8E“G2 

0.169i8r-Q2 

0.1691SE-02 

0.49999E+00 

0.49999E+00 

0..14214E-Q4 

0.14214E-04 , 
-O.1G919F-02 ; 

-0.16918E-02 

0.4q999E+‘00 

0.49999E+00 

-0.12440E-O2 

-0.14214E-04 

-0.46570E-0B 

0. 16919E-02 
0.12440E-02 

-0.16918E-02 

-0.23489E-08 

-0.22390E-08 

0.12274E-08 

0.12440E-02 

0.33422E-0S 

-0.12440E-02 

eigenvalues 

0.37975E+00 

EIGENVECTOR 

FrEqUENCY = 0 

,87932E+03»HZ. 

0.87165E-04 

-0.74343E-03 

-0.11064E-04 

0.11077E-O2 

0.82126E-03 

-0.28351E+00 

-0.64778E+00 

0.3ftl49E-n4 

-0.38151E-04 

-0.74342F-03 

0.H077E-02 

0.64776E+00 

0.28352E+00 

0.15233E-82 

-0.87168E-04 

-0.91G32E-03 

-0.H038E-04 

-0.66667E-Q3 

0.82126E-03 
-0.23273E-02 ; 

0.89807E-07 

0.23273E-02 

0.66670E-03 

0.9l024E-n3 

-0.15232E-02 

eigenvalues 

0.37975E+00 

EIGENVECTOR 

FrEqUENCY= 0 

,87932E+03»HZ. 

-0.3«l52E-04 

0.82128E-03 

0.11077E-02 

0.11050E-04 

0.74341E-03 

-0.64776E+00 

0.28350E+00 

0.87167E-04 

-0.87171E-04 

0.82128E-03 

0.11093E-04 
-0 • 2p.353E+00 

Q.64778E+00 
—0 • 66665E-03 

0.38148E-04 

0.23272E-02 

0.11077E-q2 

-0.15232E-02 

0.74341E-03 
“0 »9l034E“03 

O.I 3 I 28 E-O 6 

0.91021E-03 

0.15233E-02 

-0.23274E-02 

0.66672E-03 

etgiwaluE= 

0.28821E+00 

eigenvector 

FrEqOFNCYz 0 

.76604E+03»HZ. 

* 

0.7 1 U24E-04 
0.41952E-03 

0.41956E-O3 

0.41954E-03 

-0.41951E-03 

-0.50000E+00 

0.49998E+00 

0.71026E-04 

0.71021E-04 . 
-0.41953E-03 

-0.4i953E-03 

0.4q999E+00 

-0.49997E+00. 

-0.12207E-n2 

0.71024E-04 

0.24255E-02 

-0.4i950E-o3 

-0.12207E-02 

0.41954E-03 

0.24254E-02 

-0.50403E-02 

0.24255E-02 

-0.12207E-02 

0.24254E-Q 2 

-0.12207E-02 , 

eigenvalues 

0.20602E+00 

eigenvector 

frequency^ o 

• 64767E+O3 » HZ . 

-0.23805E+00 

-0.2fil84E+00 

0.28184E+00 

-0.281a3E+n0 

-0.28183E+00 
0. J.7558E+00 

0.17558E+00 

0.23805E+O0 

0.23805E+00 . 
0.28183E+00 

0.2a183E+00 

0.17558E+00 

0-17558E+O0 
0 » 60366E-01 

-0.23805E+00 

0.26440E-07 

-0.28183E+Q0 

-0.60366E-O1 

0.28184E+00 

-0.10817F-07 

-0 • 13894E-07 

0 • 331q8E— 07 

-0.60366E-01 

-0.69448E-08 

0.6Q366E-01 
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EIGENVALUES 0.97947E-01 


0.12541E-01 

•0..14430E+0G 

•0.10926E+00 

•0.42720E-02 

•0.34003E-07 


0;423l4E+00 
.-0*A.0?26.Ei0^_ 
-o • iqir6e+oo 
_rLCL.24i5.9E- 02. 
0.p939lE-01 


ei genvector 

0.41250E+00 
0 .1Q18 6E+0 0 
■0.12541E-01 
_0 . 62 194E-Q 1 
■0/57600E-01 


Appendix D 

FrEquENC Y=~ 0.44 657F+03»HZ o 


0.42688E-02 

_O_.299O.2E+0O_ 

0.42314E+00 

0.576QQE-nl 

■0.82l94E-nl 


-0.2g902E+00 
JL.0_.lA4 3 0E+ 00. 

0.41250E+00 

JL0_._8_9A91E.-01. 

0.24158E-02 


0.89391E-01 0.24158E-02 

!_0_».2.4.i58EJJ3 2 _r_0_, .8.9 391 EAol- 


0.82194E-01 
J-Q ..5.760 0 E- 01- 


EIG'eNVALUEs 0.97947E-01 • - FrEquENC Y = 0 .44657ET037HZ^ 

i . eigenvector 

- -0 . 29902E+.0 0 -0.1'G926E + n0 0.1 4,430 E. 1 0.0 -0 ,,1.0 186 E + ojO - a.l 2 5.41 E-01_ 

; 0.41250E+00 -0.42314E+00 0.42724E~02 0.12541E-01 O.4125OE+0O 

, 0.1ol86E+00 0.5T5600E-ol 0.89391E-01 0.24158E-02 0»82l94E-01 

.r 0. 3 0 0 5 OE - 0 7 -0. 8 i 1.9.4. E - 01 - 0. .2 .4.158 EllO 2__^lQ_. .89 391 E^oi -n.576nop-Ql 

_i El GeN V ALUE=_0.. 2-903.9E.t01 FrIjQUENC1=_0-^2A316E-+_03.»-HZ.. 

t\ eigenvector 

__Q..l 9-0.9 7 E.f.0.0 0...2 8it2.7E.-M3 Q =JL.-284.22E±£LQ aJLaaQ2£ ±0JQ n . l QQ97F+nO 

0.2«427E+00 0.28427E+00 -0 . 10907E+00 0.19097E+00 -0.28427E+00 

0 . 2 h 4 2 7E.+ 0 0 - 0..1 0 9 0-7E .+ (W 0.19097E+00_ -Q.2a427E+n O n.2fl4P7F+0Q 

0.10907E+00 -0.35377E-01 -0.11373E+00 -0.35377E-01 -O.H373E+0O 

_r.O . 3 2l 4 0 E + 0 O.U. 1 37.3 E-+JiUl_=l(U3l 5 3_7_8E.- 01 -0 .ll373E±fl-Q tQ-.:35377E-01_ 


eigenvector 


lE_QUENC-Y-r:_a..2-43l6E-403»-HZ. 


EIGENVALUE S Q .60586F-0: 


*. 0 .« 63666 E-O.. 3 . 
0.4f,582E-02 
■ 0....4 6 5.8 O.E -0 2. 
-0.2il33E-01 


_-.0.465a3£-ri2 - 

0.465 ft 0E-02 

— o. 21138 E -01 - 

-0.47369E+00 ■ 


. .0. • 2 7 . 2-4 5.E+.0..0 -U .81 Q37r- n i 

..E.IGe.Nv A L.U E s _0_,.6 04.84 E--l0.2- 

-Q....1 f, 3 3 2 E- 0.7. Q_._7A_79 3.E - o2- 

0.74795E-02 -0.74795E-02 

_0. ._7.4.7 9.8E-.Q.2 rJ3-i_2 S6.7.5Er_£L5- 

-0.20951E-05 0.92277E-05 

“Q«4a667E-05 Q . 499a9F+q O 


eigenvector 

l 0-.A658-OE.~O2 -0-^21 l_3AE-nl_ 

0.21127E-01 -0.63664E-Q3 

..0.63665E-03 -Q-.A6.5_82E~oZ. 

•0.81052E-01 -0.47369E+00 


OllOLL^O-idiA 

-LLQ.-6 3-666E-Q-3_ 
-0.46583E-02 

0.4.65 81 E- 02. 

-0.81030E-01 


Freq uency^ n.n o97E+n hz . 

EIGENVECTOR 

0.74795E-Q2 - -0 .81798E-Q6 il.-3.7-34-4E-0.7_ 

0.65948E-05 0.88572E-08 -0.74794E-02 


•0.49988E+00 0.27817E-05 0.49989F+00 

Q.10172E-0 4 -0.49 98 9E+n0 Q.6fl 5l0F~05 


EIGENVALUES 

l0_._8 7 9 9 8 E - 03_ 
0.34193E-02 
_0_. lj.252E.-02_ 
0.4ol99E-01 
lO_.535l55E_-05_ 


0.60384E-Q2 

eigenvector 

0.56 074E-f, ? Q.45 847 E-Q2 

O.U253E-02 0.17870E+00 

17-867E+0Q 0_..87.999Ejl0_3_ 

-0.142?lE+00 -0.17723E+00 

_- ( 0 .11 2 1 lF +nn 0 . 63405E+00 


0.39097E-02 

__0-._56.075Ejv)_2_ 

-0.63406E+00 

Q.l772 4 E+nO 


EI GfNvAL U Es 0 . 60384F- 0? 

— 0_. 3g0 97E.-P 2 Q..U253E-02.. 

0.45847E-02 -0.56076E-02 

j=. 0_»_5.6 0.7 6 E -0-2. .0..A02l-0£-o.l_ 

0.178G9E+00 -0.63405E+n0 

-0_*.4iJ112EllO- 5 D- t l_7l2iE±o.Q_. 


__ FREQUENCYr ( 

EIGENVECTOR 

■ 0.34192E-0 2 -Q .l7867E+n 0_ 

•0.40227E-01 -0.87994E-03 

JDL* .39.0 97E ” 02 Q.U253E-n2 

• 0.H212E+00 0.14271E+00 


» 1 l08flE+0 3 »HZ 

-n.3qQ 97F-02 
0.34193E-02 
_Q-.A58A7E-0.2_ 
0.H211E+0O 
0.14270 E±0 0_ 

, 1108q E+Q3. HZ. 

__0_.8 7.993 E.-03_ ; 

0.45847E-02 . 
_-0-». .341.9.2 E-0.2_, 
-0.17724E+00 


' FIGURE D-10, OUTPUT DATA FORMAT: EIGENVALUES AND 
EIGENVECTORS, EXAMPLE PROBLEM , 
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'EIGENVALUE?' 


.Q.13580E-02 
0.33017E-02 
r 0_t_3j5-0 1 8 E - 0.2. 
O.2O068E+00 
„0 •.1.3.70 5E-05. 


0 •6Q356E- 02 ' Fr Eq uENCY = 0 . lY0 8ftE+03 *HZ . 

: EIGENVECTOR 

_r.b.»-3.3 0l7E-n.2 -.0..330.17 Ej-l 02_ Q.._2_0_0.6aE+G.CL: -0.1 3 580£- 0 2 

0.33017E-02 0.20068E+00 -0.l3580E-o2 -0.33017E-02 ' 

0>20.066E±HO_ 0j_ 1.3.53AE“02 Q..33Q17E-n 2 -n. 3^0l6F-n2 ' 

0‘.45793E+00 -0.20494E“05 -0 . 45793E+00 -0.10528E-05 

_ -_0 »,3 78 3 2 E- 0 6 -0.45794E+00 0 .34345E -n5 n . 46793F+00 


EIGENVALUE S 0. 59567FT 


L 0_».6 o _6 0 2 . E — 0 4 _ 

L ; . ” 0 • 3 1 6 8 1 E - 0 2 
■' \ — 0 «2aS69E— 02 
_< 0.857&3E-02 
1 0.66564E-05 


EIGENVECTOR 


0.» 8j5.9 a 2 E - Q 2 _ 

-0*28569E-£2 
0.^008E+n0 
. “0 59.5E- 02 . 

-0.36155E+00 


-0.47012E+00 

-0.60605E-04 

r0_.34862E+.00_ 

-0.16280E+00 


_=.(U3.E4i>2Er-O2_.-r.0 ...33269E.rO 2_ 
Q.33269E-Q2 -0.31682E-02 . 

0 . 85981E-02 0 • 84085E-O2 . 

0 . l628 lE+n0 0 .3&154E+00 

0.34860E+Q0 0.2g670E-02 


l EIGENVALUES 0.59S87E-02 


FrEquENCY= 0.1i015E+03»HZ, 


S-0.33269E-02 
4-. 0.84884E-02. 
T-0.85981E-02 
Jl-0 .4,7 0.11 £+00_ 
T -O.27019E-O5 


-0-28570E-Q2 
j-.Q» 8.59 B 2 E- a 2 
0 » 85674 E-n 2 

0.? E ±D JQ. 

0»3486lE+n0 


_E I GE N VE_C_TO R 

0.31681E-02 0.47010E+00 

■0.85575 E^02 0.6p579 E-n4 

0.33269E-02 -0.28569E-p2 


-0.29679E-02 0.36l54E+ 0 0 


-0.60569E-04 
__0*84884E-02_ 
0.31681E-02 
_^Qj_3.4j8.6 1 E±JO.O_ 
“0.16281E+00 


EIGENVALUES 0.57948E-02 


FrEquEMCY = 0 • l0862E+03» HZ » 


0.43013E-02 
_„0.54C55E-02 
V-0.5405GE-02 
Jlr. 0 • 4 3.0 2.1 E ±0.0. . 
-0.1n2Q4E+00 


0.54057E-02 
Q • 540 5 6F-D? 
0.43022E+00 
-0 • 5b5ft lF~n I 
0.24289E+00 


EIGENVECTOR 

■0 . 54055E”02 -0 .43021E+00 
0.43 Q20 E+00 0 .4 30 13E-P2 
0.43013E-02 -0.54055E-02 
0.24290E+ 00 -0 .56552E-P1 
■0.56565E-01 0.24290E+Q0 


0.43013E-02 

-0 .54057E-02 

0.54057E-02 

CL..2 4.2.9 0 Et.08. 

-0.56556E-01 


EIGENVALUES 0.55392E-02 


0.15788E-01 
0._2fs64lE-Ol_ 
■ 0 . 2f,64lE-0 1 
0 ._38l26Et00 
0.6Q574E+00 


0.26641E-01 
0 ♦ 266 41E-0 1 
-0 • 38127E+P0 
Q. 73416E-01 
0.764fl9E-01 


EIGENVECTOR 
■0.26641E-01 
• 0.38 1 26E+00 
0.15788E-01 
0. 76468E-01 
0.73419E-01 


FrEquENCY = 0.1062oE+03»HZ, 


0.38127E+Q0 

0. 15788E -Q1 

-0.26641E-01 

0.73411E-Q1 

0.76485E-01 


0.15788E-01 
-0»2 6641 E-01_ 
0.26641E-01 
0 .76 484 E-01 
0.73413E-01 


FIGURE D-10. OUTPUT DATA FORMAT: EIGENVALUES AND: 

EIGENVECTORS, E XAMPLE PROBLEM . i , 


i\; • , 


’ \ 1.. :.. 


Appendix 



Eigenvectors listed in the sequence:, 

||d,[ , )d 2 ( , ]d 3 | , |d 4 | , jw 0 ,/TL, W 02 /TL W^/TLl} 


FIGURE D-ll. EIGENVECTOR NOTATION: NINE-BAY 
TWO-DIMENSIONAL PANEL ARRAY 






APPENDIX E 


CHLADNI PATTERNS FOR STIFFENED PANEL ARRAYS 


As described in the main text, photographs were taken of each mode excited. For the 
one-dimensional arrays, only the structure modes with fundamentals across the bay width are 
given except as indicated. Under each photograph appears the following notation: 

f = N/M : R : Q 

where N is the experimental frequency for the mode. 

M is the computed value (if applicable) 

R is the speaker condition (see figure 15, reference 1) 

Q is the number of unsupported bays (Q = 5 or Q = 3) . 

Not all of the modes illustrated were predominant. 

The modal patterns for the nine-bay specimens are also given . All of the predominant 
modes excited consisted of coupled fundamental modes for the individual bays. Higher modes 
are also illustrated. The notation indicating the experimental frequency, the computed 
frequency, and the speaker phase condition is similar to that described above. 


1 77 
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f = 126/117:A:5 


f= 108/1 01 :A:3 


f = 127/151 :B:3 


CHLADNI PATTERNS SPECIMEN SPI-1 
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f= 103/99: A: 5 




f = 120/11 8: B:5 


f = 109/151 :A:3 



f= 11 3/-: A: 5 



f= 1 27/1 33 : B: 3 


CHLADNI PATTERNS SPECIMEN SPI-2-1 


179 
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CHLADNI PATTERNS SPECIMEN SPI-2-1D 


81 
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f = 98/-: A: 3 


f = 1 08/1 28: B: 3 




f= 1 34/1 47: B: 3 


f = 249/- :B: 3 


CHLADNI PATTERNS SPECIMEN SPI-2-1D 


181 






iti 
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i 


f = 57/70: A:5 f = 6l/-:A:5 



CHLADNI PATTERNS SPECIMEN SPI-2-2 


182 
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f = 55/-: A: 5 




f = 119/-:A:5 



f= 102/-. A: 3 




f = 96/1 02:B:5 



f = 64/70: A: 3 



f = 1 12/-: A: 3 


CHLADNI PATTERNS SPECIMEN SPI-2-2D 


183 
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CHLADNI PATTERNS SPECIMEN SPI-3-1 
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CHLADNI PATTERNS SPECIMEN SPI-3-1D 


85 
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f = 54/-: A: 5 




f = 62/-:A:3 



f = 74/-: A: 3 



f = 62/74: A: 5 



f = 111/-:A:5 



f = 10l/~:A:3 


CHLADNI PATTERNS SPECIMEN SPI-3-2 


186 
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f = 93/-:A:5 


f = 74/71 :A:5 


ill 

llli 


1 


warn 

*11 


Ml ■ 

■ I 
■ 1(111 


- 

ii 


— I 

^*F< I 

■■illililii* 


f = 99/1 02: B: 5 



f = 63/-: A: 5 




f= 107/1 28: B:5 


f = 125/-:A:5 





f = 65/-: A: 3 f=7l/65:C:3 

CHLADNI PATTERNS SPECIMEN SPI-3-2D 


187 
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CHLADNI PATTERNS FOR MACHINED 9 BAY SPECIMEN 


f = 88/82: A 


f = 1 88/1 97: A 


f = 270/-: A 


f= 148/1 35: A 


f = 94/1 12:A 


f = 175/1 26: A 


188 
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f = 279/-: A 



SgaggaBfflS88E8aa^^8»888g88« 

f = 139/-: B 



f = 264/-: B 



f = 346/-: C 



f = 697/-: C 



f = 287/-: B 




CHLADNI PATTERNS FOR MACHINED 9 BAY SPECIMEN 


189 


















ipllllllp 


till®; 


Appendix E 


f = 1 44/ 1 6 1 : B 


CHLADNI PATTERNS FOR 9 BAY SPECIMEN SPII-2 


. . . ; 

■ 

1 ^ V‘' * : ■ | 

I 

I /fel® ^ ? ; 

; 
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